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MAGNIFICENT . PROSPECTS 


The extraordinary 21st Congress of the CPSU presented to the Soviet people a soundly based plan of action 
for the development of the USSR in 1959-1965, in which period the foundations of a communist society will be 
laid. 


The main task the program sets {s to raise the output from all sections of the economy to an extent such 
that it will enable us to take a decisive step towards creatingthe material basis of communism, and towards en- 
suring that the USSR is victorious in its economic competition with the capitalist countries, 


Living standards must rise greatly if the country's economic resources are increased, if the economy shows 
progress in all sections, and if productivity increases continuously, 


The scale of industrial production envisayed by the seven-year plan is truly enormous, The socialist econo- 
my will take a giant stride towards communism in 1959-1965. Industrial production will rise in the seven-year 
period by about as much as it did in the previous twenty years. 


N. S. Khrushchey's report, and the Congress resolutions, emphasized especially the role that automation has 
in the national economy. 


Khrushchev declared,"The tasks set by the seven-year plan can be accomplished only if new techniques are 
used, production is fully mechanized and automation used, and specialization and cooperation are widely intro- 
duced throughout the economy." 


",.. Avtomation must be introduced at the same time as full-scale mechanization; the two must radically 
improve working conditions and lighten labor, and must improve productivity, Automation is of great social, as 
well as economic, value in a socialist society, It changes radically the nature of work, it raises the workers’ cultural 
and technical standing, it makes it possible to remove the distinction between mental and physical labor; man's 


function is then to control the machines and instruments, to adjust them, and to plan the processes and conditions 
that shall be used." 


Automation has become extraordinarily important in nuclear power, in the conquest of space, in the grandest 
discoveries in the physical and chemical sciences, and in the rapid and revolutionary development of technology. 
It reveals to us great prospects, Our present experience enables imagining the not-too-distant future: in factories 
the entire work will be done in the best way by sets of automatic machines; in mining there will be hardly any- 
one underground; computers will do the work of clerks, and will enable us to plan the economic lives of organiza- 


tions, of towns, and of entire countries, by handling vast amounts of information; and, finally, space-ships, con- 
trolled by automatic devices, will fly. 


Having noted some advances in automation, Khrushchey pointed out that it was still not widely used, 


The resolutions of the 21st Congress state that, in the next seven years, we must carry out a broad program 
of full-scale automation, and must create fully automatic factories, process plants, and production lines, To this 
end,we must produce highly specialized and up-to-date autématic devices in the next few years, It has also been 
decided that we shall set up more than fifty model factories that will use the latest automatic systems, These 
factories will test on industrial scales new technical ideas proposed by scientists and inventors. 


These tasks can be performed only if scientists and automation specialists devote all their energies to them. 
Modern automatic systems are so complex that long and careful research is needed to design them; new mathe- 
matical methods, computers, and automatic control theory, must be used, In automation, as in virtually no other 
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branch of technology, advance 1s impossible without scientific progress, Scientific advances must be still more 
widely used in practice, 


A primary task is to develop and to extend the scientific basis of automatic control, and to do extensive 
research in all branches of automatics and telemechanics; to make a scientific “breakthrough,” and to seek new 
major technical advances that will enable us to raise Soviet automation to a new level, 


Special attention must be given to combining automatic devices with computers, to providing highly effec- 
tive automatic control systems that will ensure that the processes will be run in the best possible way; and to com- 
bining these systems with physics (as in crucial problems in nuclear physics, solid-state physics, semiconductors, 
low-temperature physics) to produce new automatic devices that incorporate the latest scientific advances and 
discoveries, Designers and technologists have before them the major tasks of producing highly developed processes 
and the machines with which to run them automatically. 


A unified state nomenclature for automatic devices must be created; this demand presupposes that there 
will be available unified standard sets of transducers, secondary instruments, amplifiers, effectors, and of digital 
and analog computers, all with strictly specified input and output parameters, 


The research institutes, design organizations and pilot plants must be extended and strengthened in every 
way, and the departments of educational institutions must be drawn into scientific work on automation, The 
seven-year plan assigns great importance to improved training for specialists in automatics. 


Wondrous prospects have opened before our native land. Soviet scientists and experts who work on auto- 
mation will devote all their talents, strength and knowledge to carrying out the tasks that the 21st Congress has 
set them, 
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com- 
ors, . 
nd THE BEHAVIOR OF A VIBRATION SERVO IN AN ELECTRICAL 
AUTOPILOT 
me I. N. Krutova 
sital 
(Moscow) 
ery 
: An electrical vibration servo is considered, in which electromagnetic clutches 
control the speed of the servomotor, The behavior of the servo is studied in the 
uto- phase plane. 
has The equation that describes the servomechanism is considered in two forms, 


one of first order, and the other of second order, The relations betyv:een the param- 
eters which define the stability limits are derived. 


I, Introduction 


The servo shown schematically in Fig. 1 [1] is of a type that is often found with electrical servomechanisms, 
The servomotor is of low inertia because electromagnetic clutches contro] the output speed; they couple the out- 
put shaft to the continuously running motor. This servo ts of relay type, and is stabilized by vibrational lineariza- 
tion, which vibration is provided by an internal vibration circuit. 
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Fig, 1. 


Now the rate of change of the servo's output coordinate is much less than ts that of the vibration loop's out- 
put, so we can deal with the vibration loop as though it were a separate closed system; we can neglect transients 
in this loop when we consider the response of the servo as a whole, and need consider only the loop's static charac- 
teristics, 


The stability is one of the main problems that have always to be considered in taking any servo, Phase~ 
plane methods [2] are used here to decide on the stability of this essentially nonlinear servo, 
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Il, The Equations of Motion 


The input to the amplifier T,(Fig. 1) receives a signal (an ac voltage) from the sliders of the sensing units 
Ry, Rg, and Rg if the aircraft deviates from its set course, The feedbsck potentiometer R,, which is coupled to the 
servo's output, is connected in a bridge circuit with the sliders of the sensing units. 


The signal passes from T, to the other amplifying 


ld, j ores stages, which are fitted with their own feedback, and 
on, which form an internal vibration loop inthe servo. The 
oo. BoD contacts on the output relays Cry and Cry in this loop 
6 control the circuits of the working electromagnets W, 
ae ‘ . and Wg. The servomotor also has two braking magnets 
By and Bg, which operate when the is 
U, <0 on, and press the disks of the friction clutches F , and 
6% at; the brake rings, The output drum 9 is thereby 


Any signal at the servo's input causes one of the 
PIU) working magnets to operate, so shaft 10 moves and lifts 

, the clutch disk off the brake ring, by overcoming the 
braking magnet, The rotation is transmitted to the 


Fig. 3. 


tube 4 and cog 7, and thence via the differential 8 to the output drum, which drum is coupled to the aircraft's 
steering, The slider on the feedback potentiometer is also coupled to this drum. 


We shall make several assumptions in order to simplify the mathematical description of the processes, namely, 


that 


1, The first amplifying stage is linear and is free from lag. 


2. The vibration loop can be replaced by a lag- free sampled-data link, which link produces as output the 
function °(t, U,) of Fig. 2 from an input U,. The sign of the pulses in the output depends on the sign of the in- 
put, while the pulse and pause lengths (8 and the respectively) are the times for which the relay contacts are closed 
and opened, These times are determined by the static characteristic of the vibration loop; t, = f,(U,), and o . 
= f_(U) (Fig. 3). These parameters either are measured, or are calculated when the loop's response is studied. 
The °(t, U,) function starts with a pause, since each channel has an initial lag of th when it is switched on. 


3. The transients in the working magnets are neglected; we assume that the clutches cause torque changes 
that are the same as the changes in the output of a pulse generator, except that there are lags r in the operations 
(Ty ON operate and rz on release), 


4. The motor's speed is taken as constant and as independent of the loading. The motion of the servo fs 
then described by these equations: 


ty 
! 
: i 
| 
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114 


U; =k, AU,, 


(2) 
(t, (Ui1>0), 


and repeats periodically 


0 when) <t<tr, (4) 


U 
ee ( 1<0) 


= O° (t U;), 


Ma=% (3) (1) 
M, (zt) (8) 

AU, = Vex — . (10) 


Here AU, is the input to T;,U,, is the external perturbation, k, is the gain of the first stage, a is the rota- 
tion angle of the output shaft, My, Mp, and M7 are the driving, braking and load torques, respectively, Us is the 
supply voltage, R dand Rb are the resistances of the working and braking electromagnets, Cq and Ch are con- 
stants for these magnets, Us = ac is the feedback voltage, a is a constant, b it a feedback coefficient, and J is 
the moment of inertia of the moving parts coupled to the servomotor's output shaft, The set of equations can be 
simplified if the pulse generator unit is combined with the clutch unit, and their joint output is considered as a 
nonlinear function (t, Uy) which determines the magnet switching. Then the on and off popods, to and tp, are 
given by: to =to— Ty + Ty ty = units; f(t) = 

= = = = = 2 
Uf, €,5; the set then becomes, for f(t) = 0, of the form 


d*z 
ar = Mq— Sign 


mg= (t, z), (12) 
m+ (13) 
0 for D(t, 1, 
°<t<t 


1 for ty<t<tp+h 
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0 for 
for tp<t<tp+t, (z>0), (15) 
and again repeats 
to= hi (12), (16) 
tp = (yz). (17) 
ab am ax Jom ax 


u, \? 
We have put (11) in, assuming thatthe controlled part and the effector are both of high inertia, 
We consider below the cases J = 0 and J #0. 


III, Case J = 0 


The speed at the servo's output at the moment when working magnets are energized is + w»/n (where wy 
is the motor's speed, and n is the reduction factor of the gearing), and zero at the moment when they are de~ener- 
gized. Taking the load torque m; = 0, we transform (11)-(17) to the form 


T, — @(t,2), (18) 
where 
T, na 
0 for 
= 4 0), 
2) for tp<t<typ+t, 
and repeats periodically 
0 for 
O(,2)={_ for ty < ty + to (z>0), (19) 
and again repeats 
=f; (72), (20) 
tp = fa (12). (21) 


Let us see how the system moves in the x, t plane. The integral curves in this plane have equations of the 
form 


z= =) (t — te) + 20, (22) 


where Xp and ty are the initial values of the coordinates, 


The three, values of (t, x), which depend on t and on the sign of x, are 1,—1 and 0; there will be three 
sets of integral curves in the x, t plane, 


The motion is determined by ©(t, x); it is made of runup sections (ts, relay energized), and of braking ones 
(tp, relay de-energized). In the latter sections @(t, x) = 0, and x stays unchanged. The transient will be traced 
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19) 


(20) 
(21) 


of the 


(22) 


ing ones 
traced 


out in the x, t plane by fitting sections of curves, for 
‘ which purpose we take the final values of the coordinates 
I} in each runup or braking section to be the same as the 
initial values in the next section, 


Suppose that the deviation is x, at t = 0 (Fig. 4). 
Now x, is negative, so ®(t, x) will take the values 0 
and 1 in turn, The changes start with a delay equal to 
SAT N the initial lag in the system, so the system is braked 
during the first perfod,and x does not change, This 
period is represented by the vertical straight line 1 in 
the x,t plane, The end point in this period {s found as 
follows, We take the initial r, as zero for ty and con- 

I Za struct a static tp = fg(yx) curve, The point where 
Pry this curve intersects the straight line x = x, gives the 

, a point we seek, This static curve is restricted to some 
£ Mm bounding value xp of x; this value divides the region 
of oscillation ( |x| < xg) from the region where the 
relay is permanently energized (|x| > xg). We there- 
fore put tp = . to find the initial delay, and hence 
# “ the initial braking period, when |x| > xp; the static 
curve is produced into this region as a horizontal straight 


line. 


Fig. 4. 


The system switches over a point I; @(t, x) be- 
comes 1, and the motion transfers to another integral 
line, This line is straight, but inclined. Here x in- 
creases, since @(t, x) = 1. We find the section of line 2 (Fig. 4) that determines this runup section as follows. 

We take t, as zero for ts and construct a static to = f;(yx) curve, The ppint where this curve intersects trajectory 
2 gives the point where (t, x) becomes zero, The system is then again braked, this corresponds to part of verti- 
cal line 3 in the x, t plane, The end of this period is defined by the point where trajectory 3 intersects f(yx), 
which curve is constructed by taking the t at point I, namely tg, as zero for tp in the third section, 


This process is continued until the final point on the last runup period, which point corresponds to trajec- 
tory 10, is found. Now x at point X lies in the region |x|< €g; the system will therefore’ not switch over again, 
and x stays fixed, On the x decreasing side we have the static to = f,( yx) curve bounded by the line x = — €, 
so to find X we have to produce this curve. This we do as follows. We extend curve f;(yx) in the region 0 < 
<|x| < e,asa straight line, which touches the curve at x = €9; while for x's that lie on the other side of the 
t axis we take to = tyg and extend f,(yx) as a horizontal straight line, Here tgp is the to atx = 0. Here € and 
€, define the insensitive zone and the zone where f;(yx) and fz (yx) are not single-valued: 


y Urmax y 


Figure 4 also shows such a construction for an initial deviation xg > 0. Here the nonlinear function has 
values 0 and —1, and the runup period corresponds to the third set of integral curves; x is reduced in this period, 


A straight line suffices for the phase representation if J = 0. Let us take the x axis as our phase line, For 
convenience we consider the x axis in the x,to plane (Fig. 5), In this plane we draw the static to = f;( yx) 
curves for x > 0 and x < 0. Sections of the straight lines defined by (22) which lie between the x axis and the 
fi( yx) curve will correspond to separate sections of the runup in the x, to plane, The change Ax that corresponds 
to any one such section is given by the projection of that section on the x axis, The entire set of sections of the 
runup for a given deviation in x can be dealt with in x, to plane. Figure 5 shows the constructions for the motion 
of the system for initial deviations X, and xg. We get a set of projections Ax;, Axg, . . ., Ax, on the x axis, which 
Projections correspond to the energized periods, and points of rest B, D, . . ., etc., which correspond to braking 
periods, 
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The phase representation shows that the deviations 
are reduced up to the last energized period, no matter 
what the deviation in x. In the last period x may pass 
through zero and start to increase. 


We thus need consider only the last period in order 
to analyze the stability, 


Stability, Critical relation between parameters, 
Figure 6 shows part of the x, t, plane near the insensi- 
tive zone, Let us consider the motion in this region 
during the last energized period. 


Fig. 6. 


The motion may be such that its trajectory intersects the f;( yx) curve within the insensitive zone, so the 
system will come to rest; this depends on the relation between the parameters and on the initial x. The inter- 
section may fall outside the insensitive zone, though; we may get self-oscillation if Ax, is then larger than 2€,, 
The stability limit is given by Ax, = —2€,. Here the system hunts only if the initial x for the last period is €, 
or —€,, If this is so, the system will follow lines 1 and 2 (Fig. 6). We get the critical relation between the param- 
eters in the form given below by using the initial and final coordinates of these trajectories, and the equations 
thereof; 


(23) 


Effects of loading. We have assumed that wy is constant and is independent of the loading, so the latter 
can affect only r; and rz, i. e., can alter only f,( yx) and fg(yx). The changes are usually slight, and can be 
neglected in studying the servo's response, It thus suffices to consider zero load if J = 0. 


IV. J #0 


The only type of load we shall consider is that of the form —myz = kx, since the motion stays the same if 
my = 0 and if the only loading is due to dry friction, 


Here (11) takes the form 


ToS — m, sign — kz. 


The motion described by (24) and (12)-(17) can be represented in the x, y phase plane, 


| 
|__| 
| (24) 


ame if 


(24) 


The phase plane, We take a new variable t* >| Tr t, and replace (24) by 
0 


dy ° dz 


We eliminate t*, my and m,,, and integrate the equations, to get the phase trajectory equation: 


y?+ {z (t)}sign {20 _ (t)]sign y \ an 


The phase plane will consist of four sheets completely superimposable, in accordance with the values of 

$(t), ¥[(t)] and sign y. It will also have boundaries that are caused by the motor's maximum speed and by the 

minimum speed during braking. In sheet I, where ®(t) = +1, the boundary will be the straight line y =y,.,, = 

in sheet II, where —1, it will be the straight line y = = — and in sheets III 

and IV, where ®(t) = 0, it will be the straight line y = 0, In each sheet the phase trajectories are circles with their 

centers at x= 41/k,y =0 andx = 4 m,/k, y = 0. 


The image point in the x, y plane will follow a composite trajectory, since it is made up of energized and 
braking sections, This trajectory we shall trace by fitting together appropriate pieces of the trajectories, We shall 
avoid having to calculate the end-coordinates on each such section by using a graphical construction to derive the 
sections. The constructions may be made directly in the x, y plane. To this end we transform the static curves 
for the sampling link by replacing t, and t, by speed changes Ay, (runup) and Ay, (braking): 


(28) 


The static curves then become nonlinear functions: 
Ay ,=Fi(q2) for Ayg=—Fi (yz) for 
Ay, = Fa(12) for y<0, for y> 0. 
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We lay out the static curves in accordance with the signs of Ayg and Ay,» and delineate the areas where 
(t, x) takes its various forms, namely, where the relay is energized continuously (|x| > xg), where the oscilla- 
tory mode occurs (€_ < |x| < x,), and the insensitive zone, |x| < €,. 


To every period of runup and braking there will correspond that section of a phase trajectory that lies be- 
tween the x axis and the static curve appropriate to the given initial values of x and y. Now the trajectories on 
any given sheet are circles, so the sections are constructed by moving the circles" centers along a straight line 
parallel to the y axis by an amount equal to the initial deviation in y; the origin of the section is the point on 
the x axis which has the x yalue appropriate to the initial x for that section of the motion (Fig. 9). 


Figure 7 shows a phase plane in which the transients have been constructed for two initial deviations x, and 
Xg. Let us consider one of these, Let the deviation at the output be x, at the start, with y= 0, These values will 
persist for a time equal to the initial delay in the system, The relays then cut in, and, since the initial deviation 
is such that x < —xp, the motion will at first be continuous along trajectory I, until the boundary is encountered 
at Y = Yrnax? the motion will then continue at constant speed along y = ymax until x = —xp is reached. The 
relays release at point B, and the brakes are applied, We construct the braking trajectory from a point B', which 
has as coordinates those of B. The relays again cut in at C, which corresponds to C*, This next section we then 
draw from point C‘, for which the initial values are xc and yc. The motion will continue along trajectory 4 up 
to point D, and then along y = y;ax. Point E’, which corresponds to the release point E, is the point where the 
tangent to the trajectory at D" cuts the line K(yx), The other sections are constructed in the same way. 


The construction does not give an exact solution because the static curves relate to steady~state x's, whereas 
x changes during each part of the motion, The errors are small, though, because the Ax in any part is small, 


Figure 7 shows that the resulting phase trajectories are very complex, 
We shall analyze how the phase plane is split up by considering how the 
points on the x axis transform, The trajectories are the same for the half- 
planes y > 0 and y < 0, so we need deal with one only. 


Stability. Critical relation between parameters, Let us consider a 
point transform, and isolate a region nn on the x axis (Fig. 7). Here Ayp> 
> Ayy, and so all initial deviations x that lie in nn will have complete 
phase trajectories that consist of sections lying between the x axis and 
Fy(yx). Here F,(yx) is the line on which the jumps between sheets I and 
Ill occur, The deviation always diminishes during the motion along the 
phase trajectories in the nn region, except perhaps in the last period, The 
equilibrium state may be one of rest, or one of hunting. The system will 
hunt if the Ax, for the last period passes outside the insensitive zone, This 
hunting corresponds in the phase to a simple bounded loop composed of 
one runup section and of one braking section in each half-plane, The sta- 
bility limits are set by Ax, = 2€,. Figure 8 shows that in this case the 
image point falls on this loop only if the initial x for this last period of 
runup is €, or —€3. The relation between the parameters that corresponds 
to the stability boundary for this nn region we find from the coordinates 
of this loop, and from the phase trajectories’ equations, as: 


(a+) =[44Y aye, +P] + (30 


The trajectories are more complex if the initial deviations lie in the region x*>|x| > x,. The limit x* 
is caused by the fact that the first release occurs at point B if |x| > x* (Fig. 7), and so the final value (xx) is 
constant, since the same trajectory is always followed. 


The complete phase trajectories will include parts of those for the nn region,if the initial deviations fall 
near this region, and so the stability problem is as for the nn region, The trajectories in the upper half-plane 
may, as the deviations increase, change in such a way that they include no part of the nn trajectories, The x, 
may then fall within or outside the nn region, In the latter case it can happen that the image point may follow 


closed trajectories in the |x| > xp region, i, e,, there may occur closed 
composite loops that include several runup and braking sections, These 
types of hunting can be detected from the point transforms diagram, 


Here the point transform will define a correspondence function x, = 
= f(X9), where x, is the initial x. We determine x, graphically, as for the 
phase trajectoric: Figure 9 shows a graphical determination of x, for one 


Figure 10 shows the point transforms diagram in a general form, Closed 
z loops in the |x| > x, region can be detected from the position of x; = 
= f(x») relative to the straight line x, =x». Loops occur if the lines inter- 
Fig. 10. sect, We have no analytic expression for f(x»), and so we cannot examine 
the function further, or pursue the bifurcation relationships for the | x| > x, 
region. 
Thus our study of this servo shows that, if the amounts of inertia of the servo's moving parts are small and 
negligible, we can get simple hunting with one relay reversal, if certain relationships exist between the parameters, 
We can find the boundary between stability and oscillation in the parameter space of €4, tg) and T, from (23). 


The servo's behavior becomes more complex if the moments of inertia are not negligible, Simple hunting, 
as above, can occur under some conditions; complex forms, which involve several relay switching operations in 
each control circuit, can occur under others, We can find the boundary between stability and simple hunting in 
the parameter space from (30). The complex types of hunting can be detected by using the method, provided that 
the parameters have certain definite values. 
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THE STABILITY OF A NONLINEAR CONTROL SYSTEM WITH 
A NEUTRAL OBJECT 


L. I. Kupriyanova 
(Moscow) 


Lyapunov ‘'s direct method is used to solve for the stability and instability of 
the transient motion in a nonlinear control system that has a neutral object. 


It is shown that it is possible to estimate how variability in the coefficients 
affects the stability limits. 


It is also shown that it is best to consider nonlinear functions of a certain 
definite class when one considers the equation for the effector in such control 


systems, 


The stability of the transient motion must sometimes be examined for nonlinear automatic control systems, 


One type of transient motion is that caused by secular changes in the neutral controlled object's parameters, 
We then have to find the regulator parameters that will ensure stable transient motion, 


This paper deals with this problem, and with the stability conditions for such a system. 


1, Let us study the transient motion of a controlled system in which the perturbed motion is described by 
the following differential equations: 


$= 
o = + + Ard, 
= + / (9). 


The first equation here is that of the controlled object; the second, that of the regulator; and the third, 
that of the effector; » is the object's coordinate, ¢ the regulator's coordinate, o is the control signal, and n(t), 
Ao, Ay, Ag and p are coefficients, 


This control systern has the following features. 


1, The equation for the object has no terms for position and damping moments, 1. e., the object is neutral 
with respect to » [1]. 

2. Here n(t), which specifies the effect of the regulator on the object, is a function of time which is de~ 
fined for 0 < t < T (where T is some real positive number). 


3. The effector’s equation is nonlinear, The p€ term in this gives the effect of the load on the effector. 
The f(o) function represents the generalized force exerted by the effector, This function is bounded, is single~ 
valued, and is always (except perhaps at o = 0) continuous; it has the properties f(0)=0,of(o) > 0, when 
o # 0. Such a function falls in class A [1}. ? 
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(1) 


It is sometimes of value to consider f(o)'s which in addition have the properties 
d 


where 9(0)= f(0)—ho, and h is a certain number, 


Such functions form subclass A‘ of class A (Fig. 1). 
We have to determine Ag, Ay and A, in this problem. 


Il, Let us transform the original equations to the 
form we need, We use a new symbol r = A, n(t), and 
a new variable ¢ defined by 


C= (3) 


where p and q are as yet arbitrary constants, 


We also denote n(t) for short by n. Simple mani- 
pulations, and removai of €, give us 


= bum + 
= + + ngs, 
Pim + — Ro —rp(o), 


pr = — + [p — (Aint 2), 


(Ay t+ R= —L(Amt Ate +rh. 


We transform system (4) to a form from which the stability (or otherwise) of the system (1) will follow 
directly (1), 


We take new variables; 


2 
Deum =VVaule @=12, 4 Stet, m 
k=1 


where S is any real number, and F* is real and always of positive quadratic form in the variables n, and ng: 
2 


F* = >} 


Here the coefficients ®gg» %k» and ay are constant, and are subject to Sylvester's inequalities: 


(2) 
where 
items, 
meters, 
h Ne 1 q’ 2 A 
A 
by bn m=, 6) 
rd, 
ector. 
en 
| 


a, >0, A, 0. 


Some manipulations give us 


v= —wv 


2 2 
Ww = > 2 + 


8,r=1 s=1 


V 
By, = — | Aas Agr, 
\k=1 


a, B=1 


and A as, Agr Ags are the algebraic complements to the element that stands in the a, 6 row and the s, r column 
of the determinant Ag. 


The equations of the other variables are found by differentiating Sgr and ¢. These equations have as first 
integral 


2 
= 4, (10) 


8,r—1 


V 


2 


a,pB=1 


Now (10) defines an ellipsoid in the space of €;,¢, and ¢. The behavior of the image point M(¢,, C2, ¢) 
is described by the radius vector V, which has direction cosines ¢ 4, €,, amd ¢. 


The system is stable if, for any arbitrarily given number (which may be as small as we please), we can select 
a number ) such that, for any given initial conditions [¢(0), €(0), and o(0)] which lie in some bounded region 
P, and which satisfy the condition V(0) = A, we have the inequality V(t) < A fulfilled automatically for any 
t> 0. 


We thus have to find the regulator parameters Aj that will ensure that the system (1) will have a stable 
transient motion, no matter what the initial conditions may be in region P, or whether f(a) falls in class A or 
in subclass A‘, 


The new variables have an obvious geometrical meaning, so the stability is determined uniquely by the 
behavior of one function V only. Th: transform retains all the properties of f(o), so it is sufficient that W shall 
take positive values on surface (10) for the system to be stable. 


We need therefore consider only the abbreviated equation V =-WV. 
The condition that W must be positive can be put in terms of Sylvester's inequalities; 


where 
where 
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(8) 


r column 


is first 


(10) 


BB By 
By > 0, 0, Ba 0. (11) 


II], We now derive the conditions for W being positive. In this particular case we put agg = ag; = 0, and 
then from (9) we derive the expressions for By, Qj: 


By = — Bog = — bg, = By = — + (12) 
2V 
Va 


Sylvester's first inequality (the first stability condition) we put as 


By = +>0. 


This condition is not essential, and is easily complied with. The condition is 


By Bog — Byg* = by, — + > 0. (14) 


4 


Putting bg; = x and bg = y, we express the By, in terms of x and y: 


| P B 
By q’ 22 Yu, Bis 3 (1 ) 


Substituting for Byy, Bgg and Byg in (14), we get the bounds to the region the second condition defines: 
a a? 
a3, 2? + + + = 0. (16) 


The curve (a parabola) cuts the y axis at 
and touches the x axis at 


The parabola has its vertex upwards, The vertex coordinates are 


au 


Yo = 0. (19) 


= — 


The third condition, which we still have to consider, is 


2B 120102 — Boa — + R (Bu Bs — > 0. 


(9) 
| 
— 
Cao) 
y 
Pg an 
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The expressions for Q, and Q, in terms of the parameters are [see formulas (5) and (12)): 


Qs +=[+ (Ao — (Arn 


Now (5) gives us 
Ay = (ra + Plea) = — (2 + py)- 


Then Q, and Q, are, in terms of the generalized parameters x and y, 


R=—(2—y)—F +etrh. (23) 


Substituting into (21) for By, Qy and R from (16), (22) and (23), and collecting terms of the same degree 
in x, we get the equation that defines the boundaries of the stability region in accordance with the third condition; 


Bz? + Cz? + Dz+E=0, (24) 


Here B, C, D, and E are functions of the arbitrary constants p, q, 433, agg, and S, of the parameters Ag, p, 
h, fi and n of the system, and of the generalized coordinate y. 


Fig. 2. Fig. 3. 


We take p, n, i and h as known, and as given for some moment t [1]. We then select the arbitrary con- 
stants for fixed values of Ag by assigning different values to y, and find the roots x; of (24). We thus get the 
region in the x, y plane where (20) is complied with. 


The region where (14) and (20) arecomplied with corresponds to parameters for which the system (1) is 
stable. 
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1) 


(22) 


con- 


Example, We {llustrate the method by finding the stability region for the system (1), 


We define the boundaries by putting p=q=1. These values ensure that the first condition is complied with, 
and help us to construct the boundaries for the other two conditions, 


We have to choose ay, and ags in order to construct the parabola of (16) (second condition), Figure 2 shows 
the boundaries for the second condition for p = q = 1, for agg = 1, and for various values of a,;. The hatching is 


directed towards the stability side, The region that corresponds to the second condition widens as a,; increases, 
We take a,; = 4 for our future purposes, 


It is convenient to represent B, C, D, and E in the following form when we deal with the third condition: 


B= A} (a4+ au), 


C = Af (as + +a 
D = AS (a + ay + ars) + AB (ae + 


0.25 
E A} (a1 + + Gis) — AZ hn (1 + y) + (Gs + O18) — 
The relations that give the oy, in terms of p,n and h for the system (i), and of S and y, are 
s 
= (1 — — (8 — 12 + 4p%) y* + 8(1 — — Ay}, 


da = 7g 15 (1 — (—1 + —3e%) y + (3— Te) — 


S? S? 
@s = [(7-- 89 + e*)— 3y—y"], 


1 1 
ds = 7g [20(1 — 8(1 +e) y — By"], = 1(7 — 3p) — (26) 


i S? 
dy = [10(1 — y — (12 —8p) + 


n 
= 7g (10(1 — + (— 4+ + 


We have to choose S in constructing the curves that show the boundaries for the third condition, We choose 
a characteristic for the effector that falls in class A in order to find S, Let the derivative of n(t) in (1) be n = 0; 
Og, Clg, . . ., yg Will then be zero, The boundaries are then independent of n. 


We assign some value to Ag in the second equation of (1) in finding the boundartes, We take p as being 
known, 


We solve (24) to get curves in the x, y plane, which define the boundaries for the third condition, for different 
values of S (Fig. 3), In constructing these curves, we bear in mind the real roots x;, which to some extent de- 
limit the region where the second condition applies, or which lie near it, 


idition; 
h S? is? i 2S? 
= = ig (8 — 20), an = — ain =( =) (Oy — 
i \2 S? yes? 
) is 
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We find that the regions for the second and third conditions do not coincide. Moreover, we find that the 
regions for the third and second conditions are most alike if S = 0.1, We thus put S = 0.1 in the future, 


Figure 4 shows curves for S = 0.1 for values of p, namely, Pp, p= 0 and py > py the other values are the 
same as for Fig. 3, 


The regions shown in Figs, 3 and 4 indicate that the sufficient conditions cannot be satisifed if the effector 
characteristic is of class A. 


In the future we take f(c) to be a nonlinear function of subclass A’. 


We consider now how variability in n(t) affects the stability region. We take the slope of the chord h 
(Fig. 1) and assume that n(t) has a fixed value, say np, at some time t;. 


condition 


Fig. 6. 


We also suppose that the derivative of this coefficient at t, can be zero, or finite and of either sign. 


Figure 5 shows the stability regions for h = 0, h > 0 and h< 0, It becomes clear that it is best to simplify 
the problem by taking n = const. Then we have that cig, og, . . ., &yg = 0. We also have that the second and 
third conditions give boundaries that almost coincide, 


Figure 6 shows curves for i = 0, with n= np, m > mg and n< ng, It is clear that the stability region widens 
asn increases, 


Figure 7 shows the regions for i = 0 and n = ng, for h in various positions, The region also widens as the 
slope of the chord increases (hg > h,). 


These results give us the values of Ay and A, in the regulator equations, We take a working point of coor- 
dinates x and y on the plane and determine A, and A, with A, = const [see (5)]. 
SUMMARY 


A nonlinear system with variable coefficients can in some cases be dealt with as though it had constant 
coefficients that correspond to some fixed instants in time, 1. e., it is possible to consider the stability of the 
transient motion in terms of the system stability at each fixed time t in the range 0 to T. 
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STABILITY OF THE EQUILIBRIUM POSITIONS IN RELAY SYSTEMS 


D. V. Anosov 
(Moscow) 


The stability of the equilibrium positions in relay systems {is studied mathe- 
matically, The results are the same as Tsypkin [2] found by a method that lacked 


rigor. 


1. The Problem 


Tsypkin posed the problem of stability in the equilibrium positions in relay systems to L. S, Pontryagin and 
V. G, Boltyanskii, and indicated what he expected the result to be, His arguments which led to this result gave 
cause for doubt, though, Pontryagin and Boltyanskii first examined the problem fully [1], and showed that Tsyp- 
kin's prediction was correct, I have found a simpler solution, which is given here. 


The system we shall consider has two parts, one of which is linear (described by a system of linear differen- 
tial equations with constant coefficients), and the other nonlinear — a relay (Fig. 1). We deal with the exact 
mathematical character of the relay's nonlinearity below. 


Let Xg, «» Xp denote coordinates that give the system's state; we as- 
at sume that x, is the input to the relay. The response of the system is then described 


by 


= +hif (ts) (1) 


j=1 


where f is the relay characteristic. The equation for any system that has only one 
nonlinearity can always be brought to this form, 


We assume that the system (1) has an f(x) of sign x; type. The equilibrium position lies at the origin 
(0, 0, . . ., 0) im the phase space (x, Xg, . . ., Xp). The motion of the phase point (which lies initially on the 
hyperplane x, = 0) is not defined directly by system (1) if f(x,) = sign x, because the right halves in this system 
have discontinuities at x, = 0, We have therefore to state the law of motion in order to define the relay system 
with mathematical exactitude; this is done in section 3, 


2. The Formal Description of the Relay System, and the Stability Criterion 
Let p be the differential operator; we then describe the relay system 


n 
= >} Ay sign (i= 4,..., 2) 


j=1 


| 
— 
Fig. 1. 1) Linear section; 
| @ 
| 
in the form 
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(2) 


n 
pu — (= 
j=1 


Eliminating xX, .. ., Xp, we have 


K (p) 2; + ML (p) sign x, = 0, 
where M is a number, and K(p) and L(p) are polynomials *: 
K (p) = p* + a,p""1+..., L(p) = p™ 4+ 14+... 


(L(p) = 0 is not ruled out), 
We call (3) the formal description of system (2).** 
We shall prove the following stability criterion for equilibrium at (0, 0, . . ., 0) in system (2), 


Ifr =1, it is necessary that we shall have M > 0 and that L(p) shall have no roots to the right of the ima- 
ginary axis; it is also suffieient that L(p) shall have its roots to the left of the imaginary axis, 


If r = 2,the necessary condition is supplementedby condition a; = by, and the sufficient one by that a, > b;. 
If r <3 or L(p) = 0, we have always instability, 


Let us deal briefly with arguments that are in essence the same as those Tsypkin used to derive these cri- 
teria [2]. *** 


In (1) we replace f(x) by the functions 


—1 when 


= Nx when — + (4) 


| +1 when 2 


Then (1) will satisfy all the conditions of the existence and uniqueness theorem. The limit to the functions 
of (4) as N—> op is of the form sign x;. Hence, by the study of (2) we understand the study of (1) when N—> o, 
with f being taken as a sign function, 


The equilibrium position is the origin, near which (1) is linear, The characteristic polynomial for this 
linear system is easily found to be 


K (p)+- NML (p). (5) 


We can thus show that the stability or otherwise of the equilibrium position for (2) is defined by the positions 
of the roots of (5), if N is large. The stability conditions that result give at once the criteria stated above (except 
when L(p) = 0), 


Such a deduction would be invalid, though. There is no doubt that sign can be replaced by the functions 
of (4). However, (4) has the form f(x) = Nx, only in the band |x,| <1/N, and differs very greatly from that 


* K(p) is the characteristic polynomial of the matrix || aj||7' for the linear part, 

** Note that (3) of itself has no meaning, because sign x, is discontinuous apd cannot be differentiated, Hence, 
when we speak of a system described by (3), we mean that the system 1s that of (2), and that (3) is its formal de~- 
scription, 

*** Except when L(p) = 0 (a case of little interest), Here Tsypkin's argument indicates that (2) is stable if K(p) 
is stable, 
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form elsewhere; the stability can thus be deduced only for trajectories that lie entirely within that band, The 
band narrows as N—»> op, and becomes an (n— 1)-dimensional plane in phase space; its equation becomes X,=0, 
Thus, even if we can draw any conclusions about the trajectories for the limiting system (N-> oo), these conclu. 
sions will apply only to trajectories that lie entirely in the x, = 0 plane (if such trajectories indeed exist). The 

genera! behavior of the trajectories near the equilibrium position would remain unexamined, 


3. Mathematical Definition of the Relay System 
The equations that define the phase trajectories of (2) in the half-space x; > 0 are 


n 
4c” = Di (6) 
a \ j=1 


or in the half-space x, < 0 


ti = >} — hy (i= (1) 
j=1 


‘ 


But it is not yet clear how the phase point will behave on the plane 
xy = 0, One phase trajectory abc’ of (6), and one abc" of (7), will pass 
through each point b on the plane x, = 0. On abc’ we can consider only be* 


@ ce’ (which lies* in xy > 0), and on a"bc"— only a"b (which lies in x,< 0)(Fig, 2), 
7 4, The phase point has a definite direction of motion along any trajectory, 
ie” Three cases occur; which we have depends on the direction of motion along 
is a"b and be’, Figure 2a and b, relates to the first case, Fig. 2c to the second, 
and Fig, 2d to the third, 


The conclusions ** we reach by considering these cases at tae limit 
N-—> oo**® are as follows. 


In the first case the trajectory through a” has as limit when N-> @ 
the curve a“bc". Here the limiting trajectory meets x, = 0 either as x, rises 
(Fig. 2a) or falls (Fig. 2b); it coincides with the phase trajectory of (2) if 
0. 

In the second case we must have that k, < 0; and conversely, if ky < 0, 


there are in the x, = 0 plane points to which the second case applies (e.g., 
the origin), These points form a band 


Fig. 2, (8) 


j=2 


which the trajectories of (1) approach from both sides, We would naturally expect that, when the phase point has 
reached this band, it will have to move in accordance with some law in the x; = 0 plane until it reaches the edge 
of band (8). The behavior of the point is described by 


%=0, 
= any — (i= 2,3,...,n), (9) 
j=2 j=2 


* We limit ourselves to the case where a‘bc' and a"bc™ both intersect x; = 0, for simplicity, 

** f can be a function of a form moie general than those of (4). 

e** Here we mean that the limit is approached only for some area on the plane x, = 0, We can show, however, 
that any section of a phase trajectory that is finite (in time) approaches as limit some curve that fs a part of the 
relay system trajectory when N — oo, 
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(9) 


wwever, 
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and, if bb’ is some part of the phase prninion.- of (9) that starts at point b, the phase trajec- 
tory of (1) that par-es through a* or c* (or more exactly, some finite part t of this trajectory) 
will have as limit a"bb’ or c’bb*. Thus in the second case the limiting trajectory is the same 
as the phase trajectory of (2) as it approaches band (8), and is, while it is within band (8), the 
same as the phase trajectory of (9). 


The whole of band (8) is thus filled with the trajectories of (9). The mode of operation 
of the system to which (2) relates is then termed sliding, 


Fig. 3. Note that (9) can be derived formally from (1) if we equate xX, to zero, find f(x), and 
then substitute the value in the expressions for Xx, . . ., X,. The geometrical interpretation 
es the same result, We can derive the phase velocity vector y for a point b that lies in band (8) as follows, Let 
v; and vg be the phase velocity vectors at b that (6) and (7) respectively imply. We assume the vectors to emerge 
from b (Fig. 3). We join the ends of v, and v,; the point where this line meets the x, = 0 plane 1s, as is readily 
seen, the end of the vector v. 


The third case is of no interest, because b cannot be reached along the phase trajectories, This case be~ 
comes the second one if we change the sense of time (1. e.,putt=~—r), The results are of no value and are not 
given, 


As we have said, (2) can be taken as the limit of (1). We shall consider in just this way a dynamic system 
that has trajectories which coincide with those of (2) when x, # 0, and that has trajectories which either join up 
directly (case one) when x, = 0, or become those of (9) when the phase point lies in band (8) (sliding mode), A 
system so defined will be called a relay system, and will be described in the form of (2), In other words, it is 
supposed in (2) that, if the phase point reaches x, = 0, it will either pass through, or move in,sliding-mode trajec~ 
tories, 


4, Transforming the Relay System 


We are given a relay system of order n, for which the formal description has L(p)=* 0 and is of degree 
—r, Then we can use a linear nondegenerate transform of the type 


T (c,; =8,;, jf = 4; (10) 
j=1 


to transfer this system to the form 


if i<ir(for r= 1 these equations drop out), 
r+1 

is 4+ksignz, if i=r, 

j= 


(11) 


n 
>} aij 2; if i>r(for r=n these equations drop out), 
j=1 


(where y; has been replaced by xj). 


The formal description is not thereby changed. K(p) is the characteristic polynomial of || ay i, , and its 
invariance is obvious, though it is harder to show that L(p) and M are invariant, That the formal description ts 
invariant under the transform (10) is intuitively clear; the description is derived by eliminating all the xj except 
xy, while (10) does not affect x; The relay system of (11) has M = —k, while L (p) equals the determinant of 
i,j=r+1,r+2,...,m, Lastly, ifr = 2 the difference a,— by between the second coefficients 
in K(p) and L(p) equals agg in the description in the form of (11), but with the sign reversed, 


When L(p) = 0, a transform of (10) type can be applied to give 


The 
$X,=0, 
conclu- 
The ¥, 
b 

jane 
ass 
ily be* 
) (Fig, 2), 
jectory, 
along 
econd, 

n 
> oOo 
2) if 

ky < 0, : 
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when i<h, 


h 
> Ajj Xj when ij =h, 
y= (12) 


n 
aij + ky signa, for i> h. 


j=t 


We omit the proofs here. They are not difficult, but are lengthy. 


5. Proof of the Stability Criterion for r = 1 


It is easily seen that r = 1 occurs only if ky # 0. Then M ~ ky, and L(p) is the same as the characteristic 
polynomial of (9), no matter whether k,; < 0 or k; > 0 (i. e., whether or not there is a sliding mode), It is sim- 
plest to prove this by transforming to the form of (11) and using the invariance of the formal description under 
this transform. 


The definition of a relay system implies that all the phase trajectories are directed to the plane x, = 0 
near the origin when k, # 0 [more exactly, they are directed to the points in band (8)], or away from this plane, 
The origin is unstable in the latter case. Therefore, we must have that all the phase trajectories are directed to 
the plane x, = 0 near the origin if the system is to be stable; i.e,., we must have that k,< 0. 


If ky < 0, any phase point that lies near enough to the origin must eventually reach the plane x, = 0 and 
begin to move along the phase trajectories of (9). A necessary and sufficient condition for the equilibrium post- 
tion of systern (2) to be stable is that k, < 0, 1. e., that M > 0; it is also necessary for L(p), the characteristic 
polynomial of (9), not to have a root to the right of the imaginary axis, and it is sufficient if all the roots lie to 
the left of this axis. We get the criterion given in section 2. 


We now prove the criterion for r = 2; first of all we deal with a few minor topics essential to this proof in 
sections 6 and 7, The proof itself is given in section 8, 
6. Preliminary Discussion of the Phase Trajectories Near the Origin When r = 2 
Section 4 shows that when r = 2 the relay system is described by 


7, 


= + + + ksignz, (k~0), 


n 
2; (i =3,4,...,n). (13) 
j=1 
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(12) 


(13) 


First of all we may show that we must have that k < 0 (which means that M > 0); see section 4) if the equi- 
librium position (0, . . ., 0) of (13) is to be stable. To a first approximation we may take it that the k signx, term 
will play the main part in the expression for xg. Then x, and x, will behave near the origin almost as though they 
were solutions to 


(14) 


But if we take x, and xg to have positive initial values xy) and Xgq, (14) will show that, since x,(t) > 0, 
sign x,(t) = 1 and = k, then xq = Xgq + kt, xy = + Xggt + */gkt*, This means that x;(t) never becomes zero 

if k > 0; x (t) and x9(t) increase until they have got some finite distance from the origin; this itself shows that 
there is instability when k > 0, Owing to lack of space, we omit the rigorous proof that it is permissible to neg- 
lect the linear terms in the expression for x». 


Thus we need consider only the case k < 0 whenr = 2, Note that the system (13) is undefined if x, = x,=0, 
We shall define the relay system when x, = x, = 0 (i, e., we shall state the motion of the phase points for which 
X, = Xg = 0 at the start) in the next section. Here we deal with the phase trajectories that pass through points which 
lie outside the x, = xg = 0 subspace. 


These trajectories behave almost as do those of (14) (k < 0 now, remember), A phase point that is in the 
half-space x, > 0 and is near enough to the origin soon enters the plane x, = 0 and passes to the half-space x; < 0. 
Very soon again it reenters x; > 0, then x, < 0, and so on, if only it remains within the area near the origin for 
which we may neglect the linear terms in x2. A phase point x = (x4, Xg,. . ., Xp) that is at first on the hyperplane 
x; = 0, but has x2 # 0, soon enters either x, > 0 or x, < 0 (which depends on the sign of x3) and stays there for 
a certain time. 


We shall encounter these so-called derivatives of (13) repeatedly in the future, Let f(x), where x = (x,, 

., Xp), be some function of xy, . . ., Xp, Which function we suppose to be either smooth or of the form c |x| . 
Let x(t) =[x,(t), . . ., Xn(t)] denote the solution to (13) that has the initial value x(0) =x. This solution is de- 
fined if x3+ x3 # 0, which in the future we assume, The derivative function f(x) is, in respect to (13), defined 
at the point x as 


It is clear that it exists everywhere outside the plane x, = 0 and that it can be found from 


= 


df (=) 
i 


where the Vv, (x) denote the right halves of (13). > Ox v(x) is denoted by ¢(x); (x) is then defined every- 
i 


where outside the plane x;=0. In the plane x, = 0 itself we cannot define (x) and E soo), 

Let us now consider some solution x(t) = (x;(t), . . +» X(t) to system (13) and suppose that x4(0) + x9(0) # 0 
and that | x(0)| ts sufficiently _. Then x(t) will cut the x, = 0 plane only a finite number of times within a 
short interval 0 <t <ty, Hence — [x(t) Jand g[x(t)] are functions of time within this 
interval, and, except perhaps at a number of points, = 


x * Things would be otherwise if there existed an x(t) such that x,(t) = 0 for a certain interval. We would then 
be unable to say that < f[x(t)] was the same as [x (t)because y(x) would be undefined at x; = 0, although 


here it would be because v,(x) contains sign xy. 

Consider the following example to illustrate this, The system with x;= —2signx, and kg = sign x, is clearly 
unstable, But if we were to introduce Lyapunov's function W = |x,| + |x_] and to calculate its derivative in 
accordance with the usual rules, we would get W = —1, which would imply that the system is stable, In fact 
W [x(t)] = W[x(0)] — t when x(t) # 0, so W[x(t)] would be zero for a certain time, 
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We shall in the future repeatedly encounter equalities and inequalities that contain expressions which haye 
no meanings when x; = 0. We must deal with these as follows, 


Suppose that an equality (inequality) has no meaning at a point x = (0, Xg,..., Xp). We select a solution 
x(t) to (13) that passes through x, and insert x(t) in place of x in the expressions, Then, ift is small but # 0, 


the expressions are defined, and we can say that our equality (inequality) is meaning ful for all t small but # 0, 
This is why the phase trajectories near the origin play such a decisive part. 


1. Lyapunoy's E Function 
We introduce a Lyapunov E function of the following form: 


We explain our choice of this function by noting that 
4 


is the energy integral for = Xg, = + k| sign x, (xy constant), and E = agsxyx; the Last 
term has here been found by trial, 
d 


We calculate the derivative of E with respect to (13); a convenient short symbolism u[ tai 
an #0 


2 aes « 
— sign — Zs Sign 2; + sign 21). 


In a sufficiently small area V about the origin we have the following inequalities complied with* 
a|a|+ E<e|m|+ dz (a,b,c,d>0), 

when ag < 0 

—pi<i<—ak (a,8>0), 
when ag, > 0 
| ak<E<BE (a,B>0) 
and for any a2 

—aE<E<aE (a>0). (18) 


The following theorem is readily proved by using (18) (see appendix 2), 


Theorem, Let a point b approach in an arbitrary fashion some point b in the subspace xy = x,=0, Then 
near the origin origin (more exactly, in the region V where (15) and (18) apply) the relay trajectory of (13) which starts 
from point b will, within some finite time, approach the trajectory that starts from b for the system 


(i> 3) (19) 
* These inequalities are proved correct in appendix 2. 


(15) 


(16) 


(17) 


(18) 


(19) 


But in the region V there are no trajectories that start at b = (by, bg, . . ., bn) with b¥ + bf # 0 and enter 
the subspace x, = Xz = 0 within a finite time. 


Hence the system (19) as it were completes (in any case near the origin) the definition of the relay system 
(13) at points in the subspace x, = xg = 0, 1. e., at those points where (13) is not defined in advance, In this re- 
spect this subspace plays for the case r = 2 the role that the plane x, = 0 plays for r = 1, and (19) corresponds to 
a motion like that of the sliding mode, There is, however, an important difference between the sliding mode 
and that described by (19), Unlike the mode of (19), the sliding mode ts reached in a finite time. 


8, Proof of the Stability Criterion for r = 2 


The system (13) is stable if agg < 0 and system (19) has eigenvalues with negative parts (“general stability *). 
It is unstable if agg > 0 and (19) has even only one eigenvalue with a positive real part. 


Proof, The system is unstable if (19) has even only one eigenvalue with a positive real part because the 
motion in the xy = X_ = 0 subspace is described by (19). We now prove that it is unstable if ag, > 0. We assume 
the converse; then all the trajectories that start from a sufficiently small area about the origin do not leave V, 
But (17) shows that 


(x(t) > aE (x(t) (a>), 
E (x(t) > E(x (0))e*, 


so if E[x(0)] # 0, x(t) emerges from V. 
Finally, we prove that the system is stable if agg < 0 and (19) has the above “general stability.” Here (19) 
has a Lyapunov W function (see [3], for example), for which & When —aW. Its derivative with respect to(13) is 


[ae = Lae + + +B VWE, 


=3 


2 
which is less than zero if W > 5 E, 


We now consider the family of surfaces (Fig. 4): 
od am Wal 
W< re, E=e, W= Exe. (20) 


It is clear that, if c is small enough, these surfaces lie within V, and hence that (16) is correct, Thus the 
surfaces of (20) have no contact for the system (13), and the trajectories pass within them, which shows that sta- 
bility occurs, It is easy to check that the approach to the origin 1s exponential. 


These results, combined with the discussion of section 4, show that the criteria of section 2 are correct, 


9. Proof of Instability for r= 3 


We put the form as that of (11), The proof of instability for k > 0 is in general the same as for r = 2, We 
omit the detailed proof for the general case, and illustrate the proof that instability occurs for k < 0 on a very 
simple special case: 


== = Zp, tg=—ksignz, (k>Q. (21) 


Now (21) is symmetrical about the origin, We shall show that, if x;(t), xg(t), X(t) are solutions to (21) 
with initial values x4(0) = 0, x2(0) > 0, x3(0) > 0 [xg(0) and xs(0) can be as small as desired], <fter some time 
T we shall have x,(r) = 0, x9(T) < 0, X3(r) < 0 and |x, (r)| > 2x3(0). Then at some time r’ later we shall have 


hhave ff 
lution 
t £0, 
If 
= 
= 
Then 
ch starts 
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+ T°) = 0, + T") > 0, + 1") > and + 2 > 4x(0), and then after a time 
we shall have + 1") = 0, + 7° + < 0, + 7°47") < 0 and > 2 + 
> 8x (0), etc. This proves the instability. 


Now when t = 0, x; = 0, K; = xg(0) > 0, so ift is small enough x,(t) > 0. Let r be the least positive t for 
which x(t) is zero. (We shall show that r< a). When0<t< 1, x;(t) > 9, sign x,(t)=+1,k,= —k and X5(t) = 
= X3(0)— kt. Hence, xg(t) = xg(0) + x3(0)t— kt”, xy (t) = x9(0)t + x5(0)t?— kt®, The term 
is the main one in the expression for x, when t is increasing, so x,(t) in fact does become zero at some t =r, 

We have that 


(t) = (0) 4+ 5-25 (0)? — ke? =0, 


so T > = x3(0). Substituting into the expression for x3(t), we have 
L(t) = (0) — kt < (0) — 3x, (0) = — 22; (0), 
which was what we had to prove. 


10. Proof of Instability for L(p) = 0 


Here.the form becomes that of (12), We separate out a subsystem within the system: 


= 2p, 


= In, (22) 
h 


= > Xj. 
j=1 


We select some solution £3(t), . . ., &,(t) to (22) for which €,(0) > 0; then for some time interval 0 < 
=t = T we shall as before have €,(t)> 0. We divide all the €4(t) by the same factor to get a solution x;(t), 
- +» X(t) to (22) with the x; (0) as small as may be desired, and such that when 0 <t < r we have x,(t) > 0. 


n 
Let be such that ky #0, and > 0 isso small that >) and We 
j=1 


assume that the equilibrium position x; =... =x, = 0 is stable, i. e., that, ifthe |x,(0)|,.. ., |x_(0)| are 
small enough, we have that the solution x;(t), . . ., Xn(t) to (12) has [x,(t)|,...,|xp(t)| < € for t> 0. We 
choose a solution to (12) for which x,(t), . . ., x(t) are the same as our previous solution to (22), and for which 
Xh41(0) =... = Xn (0) = 0, Thus the |x;(0)] (j= 1,..., m) may rightly be taken as so small that all the | 


n 
< € (t> 0), But then when 0 <t <r we shall have x; = ky +> > ky — Day since 
j=1 j 
x(t) > 0 when 0<t<r. Hence, x,(r) >4/,kyr > €. The result is a contradiction; the system must there- 
fore be unstable, 


This work was done in a study group guided by L. S, Pontryagin, I wish té thank him, and V. G, Boltyan- 
skii, E, F, Mishchenko, and V, Ya, Lin, for advice and assistance, Pontryagin and Boltyanskii also placed at my 
disposal their manuscript works on relay systems; these I have used in compiling this paper. 


* Ifky < 0, we change the signs of the x1; we then have kj > 0. 


(22) 


which 


Appendix 1 


Passing to the Limit in (1) for the Relay System 


a, Passage of the phase point through the plane x,=0, The first case of section 3 occurs if x, has the 
same sign because of (6) and because of (7) near a point b that lies in the x, = 0 plane. We assume that x, > 0 
near b (Fig. 2a). In other words 


D 0 
j=1 


= @, 57; + (21) (23) 
j=1 


is positive near b no matter what f(x,) is, provided that it satisfles -1< f(xy) = +1, To comply with this con- 
dition the derivatives X», Xs, . . ., X, must, from (1), be bounded at every finite point in phase space, and must, 
in particular, be bounded near b, We also have that 


is bounded near b, because this 
i. the speed of the phase point in directions normal to 


The phase trajectories of (1) that pass near b must 
4 therefore have a ratio 


V +... 
0 


(24) 


Fig. 5. which is bounded for any f(x) for which | f(x,)] <1. 


We now compare the curve a"bc', which is made up from two parts of phase trajectories (Fig. 2) with the 
phase trajectories of (1), subject to the condition that f(x,) is the function of (4) and N-> @, We mark on a"bc* 


points p and q at which it meets the planes x, = -+ and x, = +> (Fig. 5). Now,when x; =< — 7 the func- 
tion (4) is just sign xy, so the phase trajectory of (1) that passes through a” is at first the same as a“bc" (from a” 
to p). The width of the band between planes x, = — ~ and x; = += tends to zero when N—> on, and so p and 


q move towards b as N—> om, When N is very large, p is so close to b that our deductions above about (23) and 
(24) become valid. From p onwards the phase trajectory of (1) is such that x, increases (since (23) 1s positive), 
The increase in directions normal to the x, axis does not exceed QAx,, where Q is a positive number that is larger 
than (24) fs, near b. This shows that, when N is larger, our phase trajectory of (1) must meet the plane x, = 1/N 
at some point r such that the distance between p and r along the x, axis is 2/N, while the distance between them 


in directions normal to the x, axis does not exceed Q-}, Hence r—> b when N-* a, Now (4) is sign x, when 
X, > = , 8o the latter part of this trajectory (part rs in Fig, 5) is (as is qc‘ also) a part of a phase trajectory for 
the system 


n 
a= 
j=1 


The theorem on the solution's being continuously related to the initial values indicates that qc’ and rs come 
indefinitely close together as N—> (i. e., a8 q andr approach b), The result given in section 3 is thereby 
reached, 
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Note, It is not essential for f(x) to be the function of (4), The above arguments all remain valid if 
f = f(x, N) is a continuous function that does not exceed unity in absolute value, that equals —1 when x; < 


<—a(N), and equals +1 when x; > a(N), where a(N)—> 0 as N—> o, Function (4) has a(N) = ~. 


b, Sliding mode. In the second case (of section 3) we have that, near a point b that lies in band (8), the 
n n 
inequalitites a,;2;—k;>0 and >) + #1 <0. Let M be a region in the x, = 0 plane which ts 
j=l 
bounded, closed and lies in the band. Then 


n 
— ky — and — ki + 
j=1 


n 
+ > (25) 
j=1 


have a positive minimum in region M, There is then 
a region U in M in the complete space of the variables 


X4, +» Xp) , aNd a number m > 0, such that both 
-e 0 . 2 expressions in (25) have values larger than m at all 
points in U. 


Fig. 6. Let ab be a point of a phase trajectory of (2), and 


let it be in the half-space x; < 0 near the point b, 
which lies within M, We choose a small positive number €. The following arguments show how small € must be, 


On ab we mark the point p at which the curve meets the plane x, = —e€ (Fig. 6). Now when N is large (i. e., 
N>= ) function (4) is sign x; when |x,| = €, and so the phase trajectory of (1) that passes through a,at first 


coincides with ab (from a to p). Point p approaches b as € —> 0. We assume that € is so small that p lies in- 
side U, 


Consider a surface f, that is defined by 


n 
+ = Ve, (26) 


j=-1 


where f is function (4); we suppose that « < m?, If at some point on F, we have |x,| >, then f(x, = #1; 
(26) then shows that one of the expressions in (25) takes the values + ve at this point (i, e., a value less than m), 
But this means that the point lies outside U. Hence, the part of F, that lies within U must lie entirely within the 
band |x,| < 


We now pursue the phase trajectory of (1), and start from p. At p,x,;= —€, and so f(x,) =—1; hence the 
expression on the left in (26) takes a value greater than m (see (25) , and the text that relates to it . In other words, 
at p we have x, > m, i. e., > Ve. Now,f(x,) is continuous, so x; > will be correct until the phase trajec~ 


tory meets F, [see (26)]. 


We take 2V€ to be so small that the phase point does not leave U as it moves away from p for a time €. 
At the end of this time it must have met F;. We suppose that this is not actually so; then at all times during 
this interval we must have x, > ve, and so x, must-increase by more than €2Ve =2¢€, Thus the trajectory 
reaches a point with x, > (and, as before, x, > Ve). But when x, >€ wehave f(x)=1, te., 

n 
= Dd ayzythn<—m- [see (25)], which result is contrary to x; > ve. 
j=1 


Thus the phase trajectory of (1) lies on ab in the section ap; the phase point passes through p and strikes 
F, at some point q in a time less than 2¥€. Now,the phase velocities are bounded within the region U, so the 
distance between p and q tends to zero as ¢ —> 0 (i, e., as N-* o). This means that q approaches b as N—> o, 


The arguments are analogous for the phase trajectory cb of system (2) which approaches b from the half- 
space x; > 0. Here we mark oncurve cb the point r at which the curve meets the plane x, = €; instead of F, 
we have a surface F, defined by 


+ = —V. 
j=1 


The part of F, that lies in U lies entirely inside the band |x,| < €. The phase trajectory of (1) that starts 
from c has as its cr section the cb curve. The phase point passes through r and then reaches a point s on F, in a 
time less than Point s approaches b as N—-> 


We now show that x, falls when the phase point passes through F,. Let At be a small positive increment 
of time. Now the phase velocities are bounded in U, so there is some positive number P such that | Axy|< PAt 
in U. Moreoever, on F, we have x, = Ve, so if At > 0 and is mall enough we have, near F;, that Ax, > 4/, veAt. 
Further, at points on F, within U we have U,—1< f(x)< 1 [this corresponds to the increasing section of function 


(4)], so if Ax, is small enough we have that Af(xy)> = Any. This latter is so because A f(x;) = NAx,, and N > 
> : . These results, taken with k, < 0, give us that 


Az = ay; Az, + kyAf (#1) < Af + < 
j=l 


<2 ant || WE att Pav. 
j=1 
Hence, if € is small enough, 


This argument shows that the phase point cannot meet F, a second time within the region U when it has 
passed through F;, (i, e., through the surface x; = Ve) once already, 


If, however, the phase point is found in the region U “to the right® of the surface F, (i. e., is found in the 
region x, < Ve), then it cannot intersect the surface F, in the region U. 


In exactly the same way we can show that the phase point cannot meet F, a second time within the region 
U when it has passed through it (i. e., through the surface x; = — Ve) once already, If the phase point lies to the 
left of F, in U (1. e., in the region x; >— Ve), it cannot meet F, within U. 


Hence, if the phase point lies between F, and Fz within U (1. e., in the region — ve < &, < ve), it must 
remain between them till it has left U. Returning to the phase trajectory that passes through a (and correspond- 
ingly, c), we find that, from q (and correspondingly, s) onwards, the trajectory lies between F, and Fy (i, e., in 
the region — Ve < x, < Ve) until it leaves U. In other words, we have, from q (and correspondingly, s) onwards, 
that the following relations are complied with 


Inj<e, Inl<Ve, (27) 


until the trajectory leaves U. 
Now (27) implies that 


n 
j=2 


where 11 e., @ is a quantity of the order of Ve, Hence, if (27) is complied with, 
we have 
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n k “n 
+0, (i = 2,3,...,n), (28) 
j~-2 j=2 


: where © is a quantity of the order of Ve. It is clear from (27) and (28) that every finite section of the phase 
trajectory that starts at q (or s) and lies in U approaclies the part bb" of the phase trajectory for (9) that starts 
atb ase —> 0 (i. as N—> o). 


y Note, It is not essential for f(x,) to be function (4), The above argu- 
<a ments all remain valid if f = f(x,, N) is a continuous function which satisfies 
these conditions: 1) the function has three parts (Fig. 7), namely | (the straight 
2 0 line y = —1), Ill (the straight line y = 1), and the middle section II, which 
meets the ordinate axis, but not necessarily at the origin; and 2) in section I] 
a f must rise monotonically and satisfy the condition | Af | = K(N)| Ax;|, 
where K(N)—> as N— o@ [K(N) = N for function (4)]. 


For example, if y(x,) is a function that satisfies the first of these con- 
---+-! ditions, and has a positive derivative in section II, then f(x,, N) = 9(Nx)) 
satisfies both conditions, Further, any function that satisfies these conditions 
Fig. 7. also satisfles the condition formulated in the note to point a. 


Appendix 2 


Proof of the Assertion of Section 7 
Proof of (15)-(18), It is clear that we have 


4 1; 


i 1 i 


In the region of the origin 


1. 1 i 


Hence (15) follows if a=-¢-|k|, 


Further, let ag, # 0. In the region of the origin 
| 2aes 
= | k| — an | 21 | — sign — agers Sign + Sign 
differs from| not more than */,|k|.* But 


4 
= > + | 21 


which proves that 


1 


* x, is expressed linearly in terms of x;, Xz, . . ., Xp, and so can be made as small as desired near the origin. 
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Now, (15) implies that any expression of the form A|x,| + Bx§ can, for some region near zero, be estimated 


—ak <A|x|+ (a>0). 


We apply this to (29) and get (16) and (17), and (18) also if ag #0. If ag =0, then E = aggxsks, and in 
some region near the origin [E |< A|x,|; whence follows (18), 


Proof of theorem, Let X(t) =[0, 0, X(t), . . .» X,(t)] denote a solution to (19) with initial values (0, 0, 
By, « - -» Dy) = b, which lies in V when 0 <t <ty and let x(t) =[x,(t), . . ., Xn (t)] denote a solution to (13) 
with initial values x(0) = b = (by, . . ., by), where b? + bg# 0. We have to show that for any € > 0 we can find 
a 5 > 0 such that when |b—b]<6 we have | x(t)—X(t)| <5 when 0 <t <t, 


Let it be given that € > 0. We can assume € so small that X(t) (0 = t = t) differs from the limit of region 
V by not less than €. If |b— bj 1s sufficiently small, then b € V and, consequently, we have tg > 0, such that 
x(t) €V for 0 <t< tg. Butsince x(t) V, we observe the inequality 


d 


E [x (t)]) < E(b) 
So, ifwe take |b—b| small enough,we can make E[x(t)], 4, e., x4(t) and xg(t) as small as we please 


n 

when 0 <t < min (tg, ty), Let 5; > 9 be such that when 0 <t the solution to >) + 
j=-3 

(i= 3), where | f;(t)| < and | €(0)—B | < 5,, differs from x(t) by not more than Denoting + 
+ 4X9 (t) by f;(t), we have that, when |b —B| is small enough, x(t) differs from X(t) by less than */, € when 
0 <t < min (ty, t,). We show finally, that t, = t. If this were not so, we would have t, =< t;, tz = min(t;, t,), 
so for all t(0 <t < t,) we would have |x(t)—x(t)| < €, and hence | x(t,)— = */,€. But from x(t,)CV 
(which denotes that x (t,) does not lie in V) we have | x(t.) — X(t,)|>€. A contradiction is reached, so the first 
part of the theorem is proved, The second part is proved similarly, by using the second part of (18) (namely, 
E >— aE). 
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STATISTICAL EVALUATION OF THE ACCURACY OF AUTOMATIC 
CONTROL OF THE GRINDING PROCESS IN A BALL MILL 


V. Sh. Bereza 


(Leningrad) 


It is shown that the powder deviates from its set grain size in a random 
fashion, The random errors in the grain size are shown to follow a normal distri- 
bution law when automatic control is used, and a law with a sharp peak when 
manual control is used, It is shown that the root-mean-square error can be used 
as an index to the accuracy of grinding control for a ball mill. 


The scatter in the grain size is one of the main indices to the quality of the output from a ball mill, Elec- 
troacoustic regulators are now in common use in many cement works [1]; these regulators lighten the work of 
the mill operators, increase productivity by 10-15%, reduce the electricity consumption, and improve the quality 
of the product, because the supply of material to the grinder is more accurately controlled. 


We first consider the effects that occur in ball-mil! grinding (e. g., of cement), and that affect the fineness 
of the powder, in order to evaluate in a scientific fashion the accuracy of automatic grinding control, and to com- 
pare such control with manual control. 


We assume that the grain size of the output powder is measured at intervals, and that the method of measue 
ment, and all external conditions that might affect the results, remain the same. In general, successive measure- 
ments differ. The size of the moving sections between the baffles changes, the size, humidity, temperature and 
hardness of the clinker vary, and so the grain size of the product fluctuates. The result is usually ascribed to mally 
small perturbing factors, each of which makes its contribution to the total error found in any one measurement, 


In essence it is impossible to predict the grain size of the cement powder exactly*; the individual sam~- 
ples show deviations from the mean typical of a random sample. But if we transfer our attention from individual 
samples to the total set of tests, we find that things are different, and that an important effect occurs; although 
the individual results are in error, the mean results from a series of adequate length show a reasonable constamcy, 
This feature enables us to use statistical methods to evaluate the accuracy of automatic controls applied to grind- 
ing. 

A systematic summary of the hour-by-hour grain sizes would give us a set of statistical data from which to 
judge the efficacy of the control. 


Mathematical statistics can give us reliable deductions from statistical data; the accuracy of the control 
can be gauged, 


“The errors in the grain size are divided into systematic and random, The systematic include, for example, 
failure to follow one of the rules for grain-size measurements, or failure to allow for changes in the regulator 


* The existing State Standard, GOST 970-41, specifies that the grain size shall be determined from the residues 
in two sieves, namely, No. 0200 (gap between wires 200y ) and No, 0085 (gap 85 4). The weight of the residue » 
in No, 0085 (in%) is the main index to the grain size as measured by the sleve method, 
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settings, for lack of material in the hopper, etc. Errors which show a regular trend are also systematic, and not 
random; they will not be dealt with here, The random errors have many causes, which cannot be accounted for 
precisely, and which act differently in every case, They cannot be eliminated, and allowance made only for 
the average,by using the laws governing random errors. 


Table 1 gives the hourly data for the errors in the grain size; these data have been taken from the logbook 
of the Central Laboratory at the Vorovskii Cement Works in Leningrad; they relate to the periods May-July 1954 
(manual control) and September-December 1954 (automatic control) and to mills Nos, 1 and 2, each for 1000 hrs, 


TABLE 1 


Number of errors 
Limits to - manual control automatic control 
erLOrs as % mill Mi mill M2 il Mi il M2 
of residue in a m 
a No, 0085 
steve + |— |total | + total| + + total 
0O—0.5 |1411137] 278 |124| 4114 238 218 183 401 220 194 | 414 
0.5—1.0 |103/142] 245 1132] 119 251 451 157 308 131 140 | 271 
1,.0—1.5 65] 88} 153 88 85 173 99 87 186 90 106 | 196 
1.5—2.0. | 60] 58] 118 63 67 130 41 34 75 43 42 85 
2.0-—2.5 41} 33 74 39 33 72 9 13 22 15 10 25 
2.5—3.0 37] 17 54 27 15 42 2 6 8 2 4 6 
3.0—4.0 20) 13 33 35 20 55 0 0 0 3 0 3 
4.0—5.0 18} 13 31 9 4 13 0 0 0 0 0 0 
5.0—6.0 6 3 9 7 6 13 0 0 0 0 0 0 
6.0—7.0 0; 3 3 1 7 8 0 0 0 0 0 0 
7.0—8.0 0; O 0 0 2 2 0 0 ‘0 0 0 0 
8.0—9.0 0 1 4 0 2 2 0 0 0 0 0 0 
9.0—10.0 0; O 0 0 0 0 0 0 0 0 0 0 
10.0 0 1 1 0 1 | 0 0 0 0 0 0 
491 | 509| 1000 |525| 475 | 1000 | 520 | 480 | 1000| 504 | 496 | 1000 


The data show that the grain~size errors are random during routine production, . 
The errors have been calculated as the deviations of the actual values from the arithmetic mean, 


Systematic errors caused by wear to the grinders have been reduced by calculating the arithnietic means 
over 25-hr periods. The means have been calculated for identical dial settings when the automatic control has 
been subject to adjustment, 


The frequency of any random error is the number of times that error is found; Table 1 therefore shows that: 


1) small errors are more frequent than large ones, i, e,, the frequency of an error falls off as the error’s 
magnitude increases; 


2) the frequency is independent of the sign of the error, i, e,, errors of equal size but opposite direction are 
equally frequent, 


These two features are exactly those of random errors, 


Now the random errors are the results of many independent causes, so we may well assume that these errors 
will show nearly normal distribution (the ideal form) if sufficient data are available. 


The normal (Gaussian) law is 


where x is a random error, of frequency y, and o is the root-mean~square error. 
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Table 2 gives the o that result from working up the data, 


TABLE 2 
I Actual Root-Mean~Square Errors for 1000 Hours Running 
Type of control Period of observation Mill Actual o 
Manual control May-July 1954 No, 1 1,82 
No, 2 1.92 
Automatic control September- December No, 1 0.97 
1954 No, 2 1,03 


We shall show that these random errors follow a normal law, The sets of 1000 hourly observations for man- 
ual and automatic control were worked up. 


Figures 1 and 2 were used to compare the actual error distributions with those predicted by the Gaussian 
law. The actual and Gaussian curves coincide to a first approximation. 


1 L i L 

“7-6 -S-$6-3 -2-7 01234567 -8-7-6-§-4-3-2-1 01 234567 
Fig. 1, Histograms of the (1) empirical and (2) theoretical distributions for manually 
controlled production; a) for mill No. 1, b) for mill No, 2, The abscissas are the errors 
as percentages of the residue in a No. 0085 sieve. 


Kolmogorov's test [2] was used to determine more rigorously how closely the two curves agree, 


Let F,)(x) denote the summed frequencies for N of the empirical data, and F(x) the integral function of 
the normal distribution. 


We assume that N is large enough (1000), and that F(x) is continuous, Then the probability that the largest 
deviation of Fyy(x) from F(x) will exceed some set value } /N given approximately by 


P(X) = P{| Pr (2) — F (2) lmax> 


P(A) has been tabulated as a function of d [3]. 


We compare our Fy (x) with F(x) to find the largest deviation Dof the one from the other, We multiply 
this largest D by VN and get the number ) = DVN, from which we can read off P(A) from the tables [3]. 


Improbable events are impossible in practice, so we may consider the discrepancies between F(x) and 
F(x) systematic, and not random, if the P()) for an observed ) is small (in practice, not larger than 0,05), Con- 
versely, if P(X) is not small, we can consider the discrepancies between Fyy(x) and F(x) accidental, so the actual 
and ideal distributions agree closely. 
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TABLE 3 


Fig, 2, Histograms of the (1) empirical and (2) theoretical distributions for auto- 
matically controlled production: a) for mill No. 1, b) for mill No, 2. 
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Limits to the 
x integrals 


Empirical 


Summed | Summed 
empirical | frequency 


frequencies Fy (x) 


Integral distri- | 
bution F (x) 


y (x) —F(x)| 


+ 
on 


LL 


act 
Ty 


Mill No, 1 
0 0.000 
6 0,006 
19 0.019 
53 0.053 
140 0.140 
297 0.297 
480 0.480 
698 0.698 
849 0.849 
948 0.948 
989 0.989 
998 0.998 
1000 1.000 
1000 1.000 
Mill No, 2 
0 0.000 
0 0.000 
4 0.004 
14 0.014 
56 0.056 
162 0.162 
302 0.302 
496 0.496 
716 0.716 
847 0.847 
937 0.937 
980 0.980 
995 0.995 
997 0.997 
1000 1.000 
1000 1.000 
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34 
87 
457 
—0.5 183 
218 ; 
99 
+1.5+ 41 
+2.04+2.5 9 
+2.54+3.0 2 
4+3.0+3.5 0 
—4, 0 
0 
4 
10 
42 | 
106 1 
140 
194. 
220 
134 | 
90 
44 43 
+2, 15 
+2. 2 
4+3,.0+3.5 3 
+3.54+4.0 0 
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TABLE 4 


D=\F y(x)—F (x)! 


Integral dis- 
F (x) 


tribution 


Summed 
frequency 


Mill No, 1 


Summed 
empirical 
frequenci 


Empirical 


Limits to the 
x integrals 


1 

0 

1 

0 

1 

2 

2 

7 

6 

4 

9 

7 
33 
58 
88 
142 
137 
141 
103 
65 
60 
37 
12 
8 
11 
7 

4 
2 


Mill No, 2 


1 0 0.0000 0.9010 
1 0 0.0000 
2 0 0.0000 
2 0 0.0000 
—7 3 0 0.0000 
5 0 0.0000 Omp050 
6 0.0010 
8 0 0.0025 omp055 
—5 45 0 0.0070 
—4 21 0.0135 08075 
—4 25 0 0.0270 20 
34 0.0490 
— 51 0.0840 330 
— 84 ) 0.1350 510 
142 0.2040 320 
—{ 230 0.2900 
. —{ 372 0.3915 200 
509 0.5000 
650 0.6085 10 
+ 753 0.7100 430 
+1 5 818 0.7960 220 — 
+1 878 0.8650 — 130 
+ 5 919 0.9160 30 
+3 0 956 0.9510 50 
+ 5 968 0.9730 50 
+ 0 976 0.9865 100 
+ 5 987 0.9930 
+ 0 994 0.9975 0040 
+ 5 998 0.9990 0010 
+ 0 41000 4.0000 0000 
1 1 0000 0.0010 
0 1 OmO01 0 00 10 
2 3 00 30 
4 bi 0 
5 0m005 0 0 050 
5 13 00010 120 vari: 
4 17 00020 4150 size 
4 18 135 
2 20 0020 00095 105 trol | 
2 22 0m0180 40 The 
7 29 0m029 060345 55 To 
13 42 0042 0m0580 160 
: 15 57 1385 
33 90 0m1475 0575 
67 157 57 160 590 
85 242 1585 
419 364 0m3970 360 
414 475 0m5000 250 to sp 
124 599 06030 0 error 
132 731 31 0m6995 315 
88 819 0m7840 350 
63 882 OMS82 08525 295 the 
39 921 921 0m9045 165 c 
27 948 099420 peak 
24 972 0m972 65 
11 983 983 09820 10 
6 989 45 
148 


TABLE 4 (continued) 


Summed Summed Integral dis- 
Limits to the | Empirical | empirical: | frequency | 
x integrals y tribution D=\F py (x)—F (x) 
frequencies Fry (x) F (x) 
+4,5+5,0 3 992 0,992 0,9955 0,0035 
+5,0+5,5 6 998 0,998 0, 9980 0,0C00 
+5,5+6,0 2 1000 1,000 0, 9990 0,0010 
+6,0+6,5 0 1000 1,000 0, 9995 0,0005 
+6,5+7,0 0 1000 1,000 41,0000 0,0000 


Table 3 gives the data that show that the errors for automatic control obey a normal law. 
Table 4 gives the data that show that the errors for manual control do not obey a normal law, 
Table 5 gives the results from applying Kolmogorov's test to the results for automatic and for manual con- 


trol. 
TABLE 5 
Correspondence between em- 
Mill er A= DVN P(A) pirical and normal distri- 
butions 
Automatic control 
Mill No, 1 0.020 0.632 0.82 Present 
Mill No. 2 0.029 0.916 0.36 og 
Manual control 
Mill No, 1 0.062 1.960 0.0011 Absent 
Mill No, 2 0.059 1,860 0.0023 


It is clear from Table 5 that automatic control gives a normal law, whereas manual control does not. 
Pearson's x? test [4] gives a similar result. 


Thus only automatic control gives a distribution that agrees exactly with the theoretical one, Manual 
grinding gives a peak sharper than the Gaussian function predicts; conditions are unstable, and the dispersion is 
variable, Peaked curves are obtained as the dispersion increases [5], which shows that the scatter in the grain 
size can be reduced by making the grinding process more stable (e, g., by using automatic control), Manual con- 
trol gives a distribution that differs from the normal one, Let us give a quantitative estimate of the difference, 
The manual-control curve shows no appreciable skewness, so we consider only the sharp peak and broad wings, 

To this end we use a dimensionless index (the “excess*), which in [6] defined as 


Ey =" —3, 


to specify the sharpness of the peak; here, p, is the fourth-order central moment, and o is the root-mean-square 
error, 


The normal distribution gives zero excess, A finite value indicates a deviation from the normal law. If 
the curve is nearly normal, a positive excess implies that the probability density curve has a higher and sharper 
peak at the mode (the most frequent value) than does a normal curve with the same root-mean-square error. 


The fourth-order central moment is given in terms of the mathematical expectation M for x by 


a = M (2 — 
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TABLE 6 
Deviation of the Size and Sign of the Excess for the Curve Given by Manual Grinding Control 
x xi P(x,) xt P(x,) | x; x! P (xj) xt 
Mill No, 1 | Mill No, 2 
—10.0 {10000 0,001 10,00 —9.5 8100 0.001 8.100 
— 9.5 8100 0.000 0.00 —8.75 |5850 0.000 0.000 
— 8.5 5220 0,001 5.22 —8.25 | 4630 0.002 9.250 
— 7.5 3160 0.000 0.00 —7.75 | 3600 0.001 3.600 
— 6.75 | 2080 0.001 2.08 —7.25 | 2760 0.001 2.760 
— 6,25 1520 0,002 3.04 —6.75 | 2080 0.003 6.230 
— 5.75 1100 0.001 1.10 —6.25 | 1520 0.005 7.600 
— 5,25 760 0.002 1,52 —5.75 | 1100 0.004 4.400 
— 4,75 510 0.007 3.56 —5.25 760 0.001 0.760 
— 4,25 327 0.006 1,96 —4.75 510 0.002 1.020 
— 3.75 198 0.004 0.79 —4,25 327 0.002 0.654 
— 3.25 112 0.009 1.04 —3.75 198 0.007 1.380 
-- 2.75 57 0.017 0,97 —3.25 112 0.013 1.450 
— 2,25 25.6 0.033 0.84 —2.75 57 0.015 0.855 
— 1.75 9.4 0.058 0.54 —2.25 25.6 0.033 0.845 
1.25 2.4 0.088 0.215 —1.75 9.4 0.067 0.626 
0.75 0.3 0.142 0.045 —1.25 2.4 0.085 0.208 
0.25 0.004) 0.137 0.001 —0.75 0.3 0.119 -0.038 
0.25 0.004) 0,144 0.001 —0.25 0.004) 0.114 0.004 
0.75 0.3 0.103 0.033 0,25 0,004} 0.124 0.005 
1.25 2.4 0.065 0.159 0.75 0.3 0.132 0.042 
1.75 9.4 0.060 0.056 1.25 2.4 0.088 0.214 
2.25 25.6 0.041 1.050 1.75 9.4 0.063 0.890 
2.75 57 0.037 2.110 2.25 25.6 0.039 1.000 
3.25 112 0.012 1.340 2.75 57 0.027 1.540 
3.75 198 0.008 1.580 3.25 112 0.024 2.690 
4.25 | 327 0.001 3.600 3.75 198 0.014 2.180 
4.75 410 0.007 3.570 4.25 327 0.006 1.960 
5.25 760 0.004 3.040 4.75 510 0.003 1.530 
5.75 1100 0.002 2.200 5.25 760 0.006 4.560 
5.75 | 1100 0.002 2.200 
6.25 | 1520 0.000 0.000 
6.75 | 2080 0.000 0.000 
7.25 | 3600 0.000 0.000 
Here we have a centered distribution, 1. e., Mx = 0, so the central moment degenerates to jy = Mx‘ = v4, 
which is given by 
n 
w= (x), 
i=—n 
where x; is the mean value of the error within a given interval, and P(x;) is the error frequency for that interval, 


Tables 6 and 7 give v, for manual and automatic control, 
Table 8 gives the magnitude and sign of the excess for the two forms of control. 


It is clear that the excess for the manual curve is 3-5 times larger than that for the automatic curve, and 
(an important point) is of the opposite sign, 


This result is decisive in evaluating the accuracy of grinding. If the two forms of control gave error curves 
that were very similar (as regards their excess values), it would be enough to use the relative number of errors if the 
falling between set limits (say, : 1% of the residue in a No, 0085 sieve) to evaluate the accuracy. tt is c 


The excess values differ too greatly for this to be possible, though, It would therefore be best to use some 
measure which would summarize in a single figure all the random errors found during an observation period, The 
root-mean-square error g is suitable, Tchebycheff's test [7] shows that dispersion about the mean value {s small 
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TABLE 7 
Derivation of the Size and Sign of the Excess for the Curve Given by Automatic Grinding 
Control 
| xt P(x;) x} P(x,) Pix) | xt P(e) 
Mill No, 1 Mill No, 2 
—3.25 | 112.0 0.000 0.000 —3.75 | 198.0 0.000 0.000 
—2.75 57.0 0.006 0.342 —3.25 | 112.0 0.000 0.000 
—2.25 25.6 0.013 0.332 —2.75 57.0 0.004 0.228 
—1.75 9.36 0.034 0.318 —2.25 25.6 0.010 0.256 
—1.25 2.44 0.087 0.212 —1,75 9.36 0.042 0.393 
—0.75 0.316 0.157 0.049 —1.25 2.44 0.106 0,258 
—0.25 0.0039} 0.183 0.001 —0.75 0.316 0.140 0,044 
0.25 0.0039; 0.218 0.001 —0.25 0.0039} 0.194 0.001 
0.75 0.316 0.151 0.048 0.25 0.0039} 0.220 0,001 
1.25 2.44 0.099 0.242 0.75 0,316 0.134 0.041 
41.75 9.36 0.044 0.384 1.25 2.44 0.090 0.220 
2.25 25.6 0.009 0.230 1.75 9,36 0.043 0.402 
2.75 57.0 0.002 0.114 2.25 25.6 0.015 0,384 
3.25 | 112.0 0.000 0.000 2.75 57.0 0.002 0.4114 
3.25 | 112.0 0.003 0.336 
3.75 | 198.0 0.000 0,000 
WN 
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Fig. 3, Histograms of the errors for (1) manual and (2) automatic grinding control; 
a) mill No, 1, b) mill No, 2. 


if the root-mean-square error is small, Figure 3 shows histograms for manual and automatic control, from which 
it is clear that o, < Cm? where subscript a relates to automatic control, and m to manual, 
Automatic control thus gives a scatter in grain sizes that is only 94 /0,, of the value for manual control, 


The small o values indicate that small random errors predominate, and that large errors are very few; the 
two together are the main sign that the control works accurately, 
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TABLE 8 These data for 1000-hr working periods were worked 
up as 1000-hr series, and as 10 series of 100 hrs each, 


Cement Mill Table 9 gives the resulting root-mean~square errors, nole 
When we compare the o values of Tables 2 and 
Automatic control 9, we see that 100-hr periods are sufficient to evaluate 
the control accuracy approximately, This fact was used eo 


= res rs A in drawing up the working rules ‘or routine use of auto- 
™ matic control applied to the feeds of clinker, gypsum 
and other materials; the period used for checking on 
24 hr working was 100 hrs, The accuracy of the control 
process was evaluated from the o values given by 100-hr 
observations; the data also indicated to some extent 


Manual 
52. 


TABLE 9 
Values for Manual and Automatic Control of Grinding (Results from 1000-hr Working Periods) 


o from 1000-hour working 
Mill 


3 4 


Manual control 
| $8 1.82 | 1.92 | 1.65 | 1.78 
4.58 | 1.54 | 2.09 | 1.92 | 1.40 


0.94 1.04 
6. 1.11 146 
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TABLE 10 how well the electroacoustic regulator had been set 


up. Table 10 gives the numerical critera used for 


aan a aetieel Very Good | Satisfactory evaluating the accuracy of the automatic control, 


accuracy good 


Tables 2 and 9 show that the automatic control 
gives an accuracy about twice as good as manual con- 
trol does, 


Experience with 250 electroacoustic regulators fitted to ball mills in 60 cement works has shown that manual 
control usually gives a scatter in grain size that exceeds the permitted limits, whereas automatic control in nearly 
all cases keeps the grain size within the set limits. The scatter (as indicated by the root-mean-square error) given 
by manual control is 1.8 to 2.5 times as large as is that given by automatic control, 


Our criterion for the control accuracy of a ball mill used for grinding cement (the root-mean-square error 
a) is clearly suitable for evaluating control accuracies of other types of ball mills (e. g., drum and tube ones), 


Root-mean-square |<0.85 | <1.00 =1.15 
error 0 
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UNIFIED TWIN=-DRIVE SERVOSYSTEMS* 


P. F. Klubnikin 
(Moscow) 


The basic aspects of the theory of unified twin-drive servosystems are ex~- 
pounded, Transfer functions and error coefficients of unified servosystems are de- 
rived, and the results of their experimental investigations are given. 


Today many automatic systems of :egulation and control,as well as servosystems,are having more and more 
demanding requirements placed on them with respect to their band-width, transient times, maximum accelera~ 
tion and other dynamic indexes. 


Servosystems endowed with a high dynamic accuracy are necessary for the control of high speed industrial 
processes as well as for different types of control units and computers, It is known that the capabilities of a system, 
insofar as its dynamic accuracy is concerned, are determined, all other things considered, by the limiting values 
of the coordinates and their derivatives characterizing the behavior of the elements comprising the system, The 
restrictions imposed on the velocities and accelerations of the system's drive are of primary importance, These 


restrictions prevent one from increasing the dynamic accuracy and frequency band-width of the system. However, 


even those limiting values of the coordinates which the drive provides cannot be fully utilized in all systems; com- 
plex computers are necessary for this [1]. This can reduce the reliability of the system. 


This paper is devoted to the description of unified twin-drive servosystems developed by the author, the 
practical application of which confirmed their high dynamic qualities, These unified servosystems are based on 
the concept of parallel, combination open-cvcle closed-cycle systems, Such a system was described by Moor [2}, 
This idea was developed in the direction of the design of systerns possessing high dynamic accuracy and a wide 


band-w idth. 


It should be noted that with unified servosystems, in contrast to Moor's combination systems, considerably 
greater limiting values of the coordinates can be obtained at the output of the system, which 1s very important 
for improving its dynamic properties. 


Basic Theory of Unified Twin~-Drive Servosystems 


&) Block diagrams, Unified servosystems can, in general, consist of n loops,**each of which comprises 2 
drive and associated components, All the drives are connected to a common load through a differential reduc- 
tion gear. 


Here will be considered in very general terms only the case having the greatest practical value, that whet 
n = 2, i, e., when two drives are used in the system. 


In this case there are three possible types of unified servosystems, the block diagrams of which are show 
in Figs, 1-3. 


* Delivered at the Seminar on Automatic Control at the IAT AN SSSR, May 8, 1957. 


** In practice n can be as high as four. 
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For convenience let us introduce the following terminology. 


The system, the block diagram of which is shown in Fig, 1, 
we shall call a speed~speed type of system, In this system the 
first loop, including the drive with transfer function W,(p), is a 
speed servosystem.* The second loop, including the drive with 
transfer function We(p) is closed through the feedback loop by the 
total speed of both drives, 


The system, the block diagram of which {s shown in Fig, 2, 
we shall call a speed-angle type of system. Here the first loop 
is the same as in the speed-speed system: the second loop is more 
complex — it contains the feedback loop for the total angle of 
both drives, i, e., the systern‘s output angle. 


A speed~angle type of system can be applied very success- 
fully to the transformation of voltages from an integrator into 
mechanical displacements in the analog computation of automatic 


control systems, 


The system, the block diagrams of which are shown in Fig. 3, 
we shall call an angle-angle type of system, It can be construc- 
ted in two ways, shown in Fig. 3 as a and b. The difference be- 
tween the two alternatives can be easily determined from an exam~ 
ination of their block diagrams, 


The angle-angle type of system is a position-type servo- 
system, i, e,, that type of system which is most widely used in 
industry, 

b) Dynamics of drives connected to the load through a dif- 
ferential, Since there is present in the unified servosystem (USS) 


a differential reduction gear, let us examine the dynamic equa~ 


tions of drives connected to the load through a differential.** 


Assuming that the reduction gear does not contain any sel f- 
braking couples and neglecting the friction moments of the rotat- 
ing members of the reduction gear and differential, we can write 
the following set of equations describing the motion of the drives 
and the device acting as the load: 


= = M,— Mris, 


(1) 
Ji =r = + = 


Here 6,4 and 6, are the angles of rotation of the drives’ 
output shafts; 6° is the angle of rotation of the shaft connected 
to the load; 1, and {, are the ratios of the reduction gears from 
the drives to the differential (i < 1); 1, is the gear ratio of the 


differential; M, and Mg are the torques on the shafts of the drives; M_ is the reaction cog on the shaft of the 
differential connected to the reduction gears of the drives; J, and J, are the moments of inertia of the drives, 

with the reduced moment of inertia of the reduction gear taken into account, J, is the moment of inertia of the 
load, My 1s the load torque. 


* Any system in which the output speed is proportional to the input voltage we shall call a speed servosystem, 
** These problems are treated in greater detail in [3]. 
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In the simplest case which we shall consider first: 


M, ="M, = M,= M,—8 ©, 


where 6 is the rigidity of the mechanical characteristics of the drives. 
In this case from (1) we obtain 


(J + 2h G2) 4 =2M aig. 2) 
The solution of Eq. (2) for 70 Wilt have the form: 
where 
J +25 


In this way the starting time~constant of two drives with a differential in which ty < 1 is smaller than that 
of one drive.* The maximum acceleration of the load, determined from Eq. (2): 


2M, oigi 


is greater than in the case of one drive. 


c) Transfer functions of a USS and their dynamic accuracy. On the basis of Eq (1) and the block diagram 
for a speed=speed system shown in Fig. 1, we can write the system of equations 


p, = W, (p)e, — (p) 
p% = — (p) 
0, + 9, = 0’, (3) 
= — k,,p,, 
= k, pd — 


where W;(p) and W,(p) are the input signal transfer functions of the first and second actuator with the reduction 
gear and amplifying unit taken into account; Wj(p) and W3(p) are the load torque transfer functions of the first 
and second actuators; k--p@ is the input signal to the system; r, and Tr, are the coefficients characterizing the 
magnitude of the inertial load on the system's output (when Mz = 0). 


From Eq. (3) we obtain the system's transfer function 


0’ kW, (p) + k,Ws (Pp) [4 + ky, Wi (p)} 
Wilt + + (p) + pW, (p) [1 + ey, Wilp)] 
* Taking into account the inertial load. 


(4) 
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On the basis of (4) the expression for the error will be 


= 6 0’ = 
[1 Wi (p)} [1 + tapW, (p)] + (p) — kW: (p) 
[1 + ky, Wa + (p)] + tupW, (p) + topW, (p) [1 + ky Wi (p)) 


(5) 


From Expression (5) it is apparent that the system's error can be made equal to zero by including in the 
first loop a correction unit having a transfer function equal to 


1+ (p) + tapW, (p) 
cor (p) — ky, — ky, (P)) 
Such a method of error cancelling is not expedient in this case because of the considerable increase in the 


complexity of the system. However, this method points out the possibility of satisfying the conditions of invari- 
ance [4] in the system. 


Making use of (5) let us obtain the first coefficients of the error [5]: 


_ 1+ (0) —k,Wi(0) 


(0) [4 + ky, (0)) + ta, (0) [1 + ky, Ws — (0) 


C1 = 


{i+ kW, (0) — k,W, (0)} {k, Ww) (0) [1 Wi 
{1 {i+ k,W 2 (0)? 


(0) + (0) [1 + OD), 


where we have denoted 
d 


From the expressions for the error coefficients C, and C, it can be seen that by selecting the system's param- 
eters ky and ky so that the condition 


1 + ky,W, (0) = k,W, (0), (6) 
is satisfied, we obtain 
Co 0, 
k, WY (0) 
4 11 W; (0) (0) 


It can be shown that by proper selection of the transfer ratios of reduction gears 1, and 1, from the condition 
for obtaining maximum acceleration (see above),each of the quantities +,W4(0) and rgW$(0) cannot be greater 
than twice the value of the actuator's electromechanical time constant, Hence, the quantity C, in accordance 
with (7) is very small, 
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In this way, when condition (6) is satisfied a speed-speed-type system theoretically possesses first-order 


astatism; in practice, however, the accuracy of its operation corresponds to that of a system having second-order 


astatism. 


What has been said is confirmed by the transfer function of the equivalent open-loop system which, when 


(6) is satisfied, has the form: 
(p) + ky Wa(p) 
P + 1Q, (p) + w| 


where Y,,(p) is the transfer function of the closed-loop system, 


(p) 
Q1(p) = 1+ k, Wi (p ’ (p) = T+k, Wi ip)’ 
_ Ni(p) _, PN, (p) +14 
= pM, (p) + 1 


(8) 


For 2 system of the angle~speed type based on the block diagram shown in Fig. 2, we have the following 


expressions for and 
For this type of system we obtain the transfer function in the form 
hp PQi (P) + kpQz (P) 


0 P+ (p) + + taP*Q,(P) 


where 
(p) 
Q:(P) = Q2(P) = ° 
We determine the system's error: 


P {1 — (p) + (p) + PQ, (P)I 


— 7 
P+ kpQa (Pp) + (p) + (Pp) 


(9) 


(10) 


To evaluate the system's accuracy we shall adopt the same method as before and make use of the error 


coefficients, 
On the basis of (10) we have 
1 — k (0) 
Ce = 0, Cy = (0) 
(0) 
(ON [1 + kp Qe” 
(OP 
where we have denoted 
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From expression (11) it is apparent that if condition 


i= (0), 


or 
1+ (0) = kW, (0), (12) 


is satisfied, which is easily achieved by selecting the corresponding values of ky; andkp, the error coefficient 
C, becomes equal to zero while the value of C, is determined by the expression 


2Q (0) (0) (0) 


By considering the expressions for Qf) 0), Q2(0), Q}(0), Q3(0), as well as what has been said above, it is 
not hard to see that the values of Cy are also small. 


Let us note that, as can be seen from expressions (11), C, can be made equal to zero by the proper selection 


of the system's coefficients, However, in this case C, will no longer be equal to zero, although, it is true, it will 
still remain sufficiently small, 


Thus, the theoretical angle-speed type of system cannot possess third-order astatism without the use of 
additional correction loops, 


In practice, experiments show (see further on) that such a system is close to a system with third-order asta- 
tism. 


For version a of the angle-angle type of system (Fig. 3) we have 
e, = (p) — ky,p8,, & = hp — — k,,pb,, (14) 


where Ry(p)p is the transfer function of the differentiating unit. 
On the basis of (14) we obtain the system's transfer function 


kppRs (P) Quip) + kpQu(P) 
Qa (P) + (P) + (P) 


(15) 


Making use of (15) we can write down an expression for the error in this version of an angle-angle type of 
system 


_ PIL + (P) + (P)— (P) 
p+ kpQs (p) + (p) + tap*Q, (P) 
It can be seen from expression (16) that R,(p) can be selected in such a way that the error will be equal to 


zero, However, the achievement of this brings about appreciable complications because the expression for R,(p) 
in this case becomes very complicated, Determining the system's error coefficients we obtain 


(16) 


1 — Rs (0) Qi (0) 
kpQs (0) 


10; (0) + (0) 
(0) 


C, = 0, = 


C, = 2. 
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(0) (0) + ABRs (0) (0) (0) — 1 — (0) (0) Qa (0) — 
[Qu (0)]* 


— (0) QM (0) Qs (0) — tp (0), 


Ri 0) (P)], _, 


From these expressions it is easy to obtain the conditions for C; and C, to equal zero. Thus C; = 0 when 


or 1 = ky Ri (0) (0) 
1 + ky,W (0) = kpR, (0) W; (0); (17) 
Cz = 0 when equality (17) and condition 


(0) (0) + hy Rx (0) Qi” (0) = (0) + (0). (18) 


are satisfied, 


Analyzing (17) and (18),one can verify that the first expression is satisfied, irrespective of the nature of 
Ry(p), by changing the system's coefficients, For example, when R,(p) = 1 we obtain condition (12), In the case 
when 


R, (p) Tp+ 
where T = RC, T' = mT (i, e,, when using an RC network for differentiation), expression (17) takes on the form: 
1 + ky W, (0) =k (0). 


This last condition can be satisfied by selecting ky, or T'. The fulfillment of the second equality (18) is 
connected with the choice of the nature of the transfer function R,(p), which in this case must contain terms in 
p in the numerator, For example, when R,(p) = Tp + 1, equality (18) is satisfied when 
or 
0) (0) (0) 


hp Qs (0) ~ 


Hence, version a of the angle-angle type of system (Fig. 3) can always be made to possess astatism of the 
second order by variation of its coefficients. Third-order astatism can be achieved in this system by the selection 
of the transfer function Ry(p). in Fi; 


fici 
The transfer function of the b version of the angle-angle type of system (Fig. 3), obtained in an analogous +s 
manner, has the form: 


here | 

9: kppQua (p) + (P) simp! 
r—, (21 

(P) + (P) + aP*Q, (P) 


where 
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where 
| | 


(P) = k,,Wy (p) R(p)" Qu (P) = (Pp) Rip) 
On the basis of (19) the system's error is equal to 


+ + (P) — Qs 
P + (p) + (P) + (P) 
The expression for the error coefficients become: r 


1 — kp Qis (0) 
C=0 


Qs (0) (0) — (0)] + [1 — (0)] [1 + (0)] 
1Q2(0)]* 


Qis (0) ’ 


(p) (p)— Wi (p) [1 + R™ 
(0) i+ (0) RY (0) 
12 


k2 W, (0) [R(O)]* 


From these expressions it is apparent that C, = 0 under the condition 


1 = or R() - (20) 


Coefficient C, can also be equal to zero, if, beside condition (20), the following conditions are satisfied: 


k Q1x(0) = 


kytaW ,(0)W = 0. 21 
1(0)R@(0) + T+ Ws (0) ( ) 


By examining conditions (20) and (21),we can verify the fact that the same situation exists for the version 
in Fig. 3b as for the version in Fig, 3a, Condition (20) can be satisfied in all cases by changing the system's coef- 
ficients, while condition (21) can be satisfied only by the selection of the transfer function R(p). 


However, in comparison to the version in Fig. 3a, the selection of the transfer function R(p) is complicated 
here by its effect on the stability of the first loop and the system as a whole, For example, ifR (p) = Tp +1, which 
simply corresponds to the existence of a speed- and angle- feedback connection in the first loop, then condition 
(21) cannot be satisfied because the derivative R(0) is positive. 


Obviously, the derivative Ro) can be negative, for example, in the case 


¢ 
where 

or 
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The condition (21) is satisfied when 


4 (0) 
P—Ti= + 


The values of T, and T are selected from the system's stability condition, 
The above results of investigation of USS are valid for the linear regions of operation of its elements, 


However, it can be shown that the high dynamic accuracy of the USS is preserved even when the first loop 
operates in a nonlinear, saturated region (limiting of coordinates, velocities, accelerations, etc.), i. e., precisely 
under those conditions when the capabilities of the system's drives are utilized to the limit. 


This Last situation makes it possible to obtain a wide band of signal frequencies which are passed by the 
system without any error in the presence of large amplitude input signals, 


This will be demonstrated further on by means of an example of a concrete system. 


d) Transient processes in USS, One of the features of the transient process in USS ts its oscillatory nature, 
This follows directly from the existence of astatism in the system, 1. e., the zero values of the first error coef- 
ficients, 


Another characteristic of the transient process is its short time of duration as compared to an ordinary system, 
This is explained by the possibility of obtaining a much higher acceleration atthe output of the system than ts 
possible in an ordinary one, which was pointed out earlier on in this paper. 


€) Stability of a USS, The investigation of a USS consists of two stages. The first stage consists of the 
determination of the stability of each of the system's loops separately, This is done using conventional methods 
(5). The second stage consists of the determination of the stability of the system as a whole, Here one makes use 
of the USS transfer functions listed above, namely (4), (9), (15), and (19), as well as the transfer functions of the 
equivalent open system. Transfer function (8),given 


¥ above,is an example of such a transfer function. 

; (is x. | In investigating the stabilities of USS by means 
to, 4 of the transfer function of the equivalent open-loop 

as a 4 system, it is necessary to take into account the condi- 
4 tions under which the first error coefficients (C, and 

a < | C, depending on the type of system) are equal to zero, 

as is done in (8), 

50}- = We have examined the basic theory of twin-drive 
7 USS, Now, let us examine an example of an experi- 
~ p: mental investigation of a USS in which reversible, 
electromagnetic, induction clutches [6] are used as 

drives, 


Results of Experimental Investigation of USS 


Let us list the basic characteristics of the components of a USS, Two drives with clutches of equal power 
were used, The drives were connected to the inertial load through a planetary, differential, reduction gear, Angles 
and speeds of rotation were monitored by means of potentiometers and direct-current tachogenerators, 


The voltage transfer function of a clutch drive at the input of an electronic amplifier has the form: 


ap+i 23 
(Typ + 1) (Tip + 1) (Top + 1) 


/ 


| 
R(p) = when T> 
con 
can 
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| 
curv 
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where K, is the total amplification factor, Tj, is the chectveusaiiileaiiih time constant,T, and T, are the time 
constants of the loop containing the coils conttaliiog the clutch, for increasing and decreasing current,a = */,(T, + 
+ Ty»). 


The load-torque transfer function can be represented by the expression 


k’ 


where k* is the load-torque amplification factor. 


Substituting (22) and (23) into the expressions obtained above,it is easy to obtain transfer functions, condi- 


tions for the errorcoefficients to equal zero, and other expressions necessary for designing a USS with a given 
type of drive. 


We list now typical parameters of a USS, the results of an experimental investigation of which are given 
below: Kep,Wy(0) = ky kyiy 1000, kp Wy(0) = kpkyly 1000, kp,W(0) = kT 100 (for systems of the angle- 
speed and angle-angle type), 71W3(0) = 0.05-0.15. 


a) A system of the speed-speed type. The solid lines in Fig. 4 are the amplitude-frequency and phase- 
frequency characteristics of a speed-speed type of USS (Ag*/Ag) is the amplitude-frequency and ¢ is the phase- 
frequency characteristic), For comparison similar characteristics of a conventional system are shown as broken 
curves, 
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As an illustration, Fig. 5 shows the oscillograms 
of a sinusoidal input signal processed by a USS (below) 
and an ordinary system (above), Here it can be seen that 
the USS processes the given signal having a frequency 
of 6 cycles practically without an error, while the 
ordinary system already introduces a phase shift of 
about 10° at a frequency not much greater than 5 cycles, 


Figure 6 shows the transient process in a USS when 
a fixed signal is turned on and off (1 and 2 are the 
speeds at the Output of the first and second drives), 


Figure 7 shows oscillograms which illustrate well 
the operation of a USS and its loops with sinusoidal 
input signals of different frequencies, 


It can be seen from Fig. 7 that,for a low frequency 
input signal, the second loop (curve 2) works at a slow 
rate and only eliminates occasional, small errors of the 
first loop. As the frequency is increased,the second 
loop begins to work more intensively, At high frequen- 
cies, when the first drive operates in a nonlinear, satura- 
ted mode, the second loop works most intensively with 
fuli use being made of the power, speed, and accelera- 
tion of the drive, 


It is important to note that in this case there is 
practically no error at the output of the USS (Fig. 7). 
This confirms the assumption which is given at the end 
of point "c" of the first section, which was based on a 
theoretical investigation of a USS, 


b) Angle~speed system. The amplitude-frequency 
and phase frequency characteristics of an angle-speed 
type of system are shown in Fig, 8. The characteristics 
of an ordinary servosystem are shown on the same draw- 
ing by means of broken curves, 


The simplest form of an ordinary servosystem was 
considered, with no correcting unit whatsoever to im- 
prove its characteristics, By using such units one could 
improve its dynamic accuracy and increase its frequency 
pass-band, but even in this case it would fall far short 
of the USS, 


Fig. 10, 
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Fig. 13. 
Figure 9 are oscillograms showing the effect an angle-speed type of system has on a sinusoidal input signal, 
. thereby illustrating its great accuracy, 


Figure 10 shows oscillograms of the transient process when a signal corresponding to a constant speed is fed 
to the input of a USS, These oscillograms clearly point out the presence of astatism in a USS. 


c) Angle-angle system, Figure 11 shows the amplitude- frequency and phase- frequency characteristics of 
an angle-angle type of system (version in Fig. 3a) where Ry(p) = T'/(Tp + 1). The same characteristics for the 
version in Fig, 3b are shown in Fig. 12, where R(p) = m(Tp + 1), 


It can be seen from the curves that USS of this type, even in the case when the functions R(p) and R,(p) are 
of the simplest form, possess characteristics which are approximately twice as good as far as frequency bandwidth 
is concerned, as those of ordinary systems, The characteristics of an angle-angle USS can be improved at the ex~ 
pense of a slight increase in the complexity of the functions R(p) and R,(p), as can be concluded from the above 
presentation, 
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As an illustration of the dynamic accuracy of a USS, Fig, 13 shows the oscillograms of a sinusoidal signal 
processed by a system of the type shown in Fig. 3a. Figure 14 shows the same oscillograms for the version in Fig, 3, 


It can be seen from Figs, 13 and 14 that,up to certain frequencies ,the error in both types of systems is close 
to zero, However, the type in Fig. 3a follows the given signal without error up to higher frequencies, espectally 
when R,(p) = 1. 


Figure 15 shows the oscillogram of the transient process in a USS for a given constant speed (upper oscillo- 
gram). For comparison, the same oscillogram 1s also shown (below) for an ordinary system. The oscillograms 
demonstrate the high dynamic accuracy of a USS and the presence of astatism in it. 


SUMMARY 


Theoretical and experimental investigations of USS have shown that USS enable one to obtain high dynamic 
accuracy and a wide frequency bandwidth, By a simple selection of parameters, it is simple to obtain in a USS 
astatism of up to and including the third order. USS are appreciably superior to ordinary servosystems with respect 
to dynamic accuracy, having frequency bandwidths typically twice as wide, The structure of USS is simple and 
their design.does not pres. “t any particular difficulties. 


The shortcoming of USS rests in the necessity of having two drives and a differential reduction gear in them, 
USS are a logical choice in all cases of application where it is necessary to have very high dynamic accuracy 
and a wide frequency bandwidth accompanied by a small error, 
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DIGITAL INTEGRATOR DEVICE FOR PROGRAMMING SECON D-ORDER 
CURVES 


A. A. Voronov and G. N. Sokolov 
(Leningrad) 


The feasibility of applying digital analog computers to the programming 
of second-order curves is investigated, The circuit of a programming device 
applied to a milling machine for the working of profiles consisting of sections 
of straight lines and circles is presented. 


Programs for machines with programmed control can be developed using either universal calculating ma~ 
chines or specialized devices, 


Universal calculating machines are usually used for the preliminary recording of a program in an intermediate 
memory unit (in the form of a series of pulses) with the subsequent feed-out of data directly to the device control- 
ling the machine, In this case programming is possible either from the equation of the profile's curve of from 
points of support with a given law governing interpolation between them, The complexity of the programmed 
curves is determined only by the capabilities of a given computing machine, A shortcoming of programming by 
means of universal calculating machines {s the very large volume of the intermediate memory unit, 


Of all the specialized devices, simple, linear interpolators have found the widest application. In this case 
the pro file of the programmed curve is broken up into pieces ofstraight lines by means of a series of support points, 


The more complex the nature of the profile and the higher the requirements on the programming accuracy, the 
greater is the number of points. 


Available data on specialized units for the programming of second- and higher-order curves, operating on 
the same principle as universal calculating machines, is very general, One can assume that the possiblities of 
such devices are limited since they can interpolate only along a definite group of curves for the reproduction of 
which a given specialized unit is designed, 


Digital analog computers can be used very conveniently as programming devices, In this case the shape of 
the interpolated curve will be determined by the number of integrating units in the digital analog computer and 
the circuitry interconnecting them, 


Below are described circuits and principles developed at the Electromechanical Institute of AN SSSR, 


Programming of Second-Order Curves by Means of Digital Analog Computers [1) 


In many practical cases the trajectory of the center of a cutting tool can be approximated with any degree 
of accuracy by sections of second-order curves, 


The general form of the equation of a second-order curve is: 


Az* + 2Bry + Cy? + 2Dz + 2Ey + F =0. (1) 


j 
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Differentiating the equation with respect to x and determining dy/dx we obtain the equation of the phase 
trajectory: 


Az + By+D 
dx + Cy +E (2) 


dt 


_w(At+By+D), = —o(Br+Cy +2), (3) 


where w is an arbitrary quantity characterizing the angular velocity of a point moving along the programmed 
curve, 


Equations (3) determine the operating characteristic 
of a programming device consisting of integrators, A 
block schematic of the programming arrangement is 
shown in Fig. 1, 


For the reproduction of circles the circuit is con- 
siderably simpler, The equation of the phase trajec- 
tory in this case has the form 


t— Zz, 


y—y,’ 


(4) 


where x, and y,, are the coordinates of the center of 
the circle. 


The block schematic of the arrangement for reproducing circles is 
4, shown in Fig. 2. By introducing extra controls in this block schematic in 
J z = the form of switches in the circuits connecting the integrators, one is also 
able to reproduce: a) straight lines (when switches K, and Kg are open); 
[—*¥ — b) parabolas (when one of the keys K, or Kg is open); c) equilateral hyper- 
% (iad | bolas (when both keys are closed but the sign of the feedback loop {s reversed, 
Ye Se and K; is closed), 


Fig. 2. The table lists some possible methods of using the described device 
ig. to develop a curve between two given points O and A in the first quadrant 
of the coordinate plane. 


Columns 1 and 2 of the table list the signs of the feedback connections between the integrators for the x 
and y coordinates, while colurnns 3 and 4 list the signs of the initial values introduced into the integrators, A 
zero denotes the absence of a connection, 


Operating Principles of a Digital Integrator 


Figure 3 shows the main block schematic of a digital integrator, A digital integrator consists of an integra 
tor register R— a scale-of-two reversible counter, a frequency divider FD (using a binary counting circuit), and 
voltage pulse circuits P. 


Let the following quantity be entered in the integrator‘s registor 


n 
z= 2'z;, (5) 
i=0 


where x; is either zero or unity. 


or 
| y 
4 
7 
110 


Some Types of Curves Produced by the Programming Devices 


Sign of 
feedback 


Initial values of 
quantities 


| 


Shape of 
curves 


8 4 


Straight line 


Every gate P is open if the register's cell connected 
to it is in a condition corresponding to unity, The fre- 
quency divider supplies to the gates P pulses with a re- 
petition rate of F/2, F/2®,..., 


The number of y pulses appearing at the output of 
the integrator in an arbitrarily chosen interval of time 
t, which is a multiple of F~*, is not, generally speak- 
ing, proportional to x. Only over a time r = 2"*!p~?, 
corresponding to the transition of the FD from a zero 
condition again to a zero condition, 


n 
Mat =z. (6) 
img 


9 Ps curve 
y 
0 
= + + + 0 >. | Circle 
\ 
Circle 
z 
| 
/ j 
0 z 
+ 0 + + } \ Parabola 
g 
+ + + Hyperbola 
; 

+ + + + Hyperbola 

Fig. 3. 
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However, the maximum error A, measured along 
a perpendicular to the straight line z = xt/r for any 
discrete value of t/, is approximately equal to only 
one y step, Figure 4 shows the integrator's construction 
of a straight line when x = 1010 ... 0. 


Hence, when x = const, the relative error can be 
made, for a given t/r, as small as required by decreas 
ing the size of every step which can be achieved by 
the corresponding selection of the value of x and the 

number of rows n in the register. Thus, for example, 
Fig. 4. if the scales of the discrete values of Ay and t/r along 
the coordinate axes are equal, then a straight line with 
the tangent of the angle of inclination equal to */, can be constructed with x = 5 and n= 2, x = 10 and n= 3, ete, 
The value of the maximum error A will be correspondingly less than 1/5, 4/4, etc. If x varies during the integra- 
tor’s operation then the expression 


F n 
y >= 
2 
has to be written ip the form 


F 1 F k=N 
= > (332!) At = >) At, 
k=0 


where At is the time interval between the pulses with a frequency F, The last expression can be represented, by 
analogy to continuous integration, as 


t 
F 
0 


When x = x, sin wt, the error will be as large as the ratio w/ ar Plotting shows that when w < our 
the error A usually does not exceed two-four steps, When w « oar the error decreases appreciably and as 
w tends to zero, the maximum error tends to a value approximately equal to one step, 

In this way, the error in digital integration can be reduced first, by decreasing the size of the step, and 
second, by decreasing the rate of change of the integrand, 


Programming Arrangement for Reproducing a Straight Line and a Circle 


An example of a programming arrangement based on digital integrators with successive carriers (binary 
multipliers) is shown in block diagram form in Fig. 5. The apparatus is designed for the reproduction of straight 
lines and circles by the method of programmed control of a two-coordinate milling machine. 


Before the programming equipment is put into operation, the integrator registers Ry and R, and the displace 
ment-measuring devices DM-1 and DM-2 are set at zero, (The circuits for setting the displacement meters at 
zero and introducing the initial conditions into them are not shown, They are analogous to similar circuits for 


the integrator registers.) 


For this, the starting pulse S_- puts trigger A in such a condition that the gates P, which it controls are opetied. 
After an interval of time determined by the delayed multivibrator ZM, a negative pulse, formed by the NE clr 
cuit, passes through the gates P, to the trigger units of registers R, and R, and the displacement meters, thus putting 
them in a zero condition, 
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Once the zeros are set and after an interval of time related to the firing of ZMg, trigger A closes gate P; 
whiic trigger B opens gates P', P, and Py, 


A pulse displaced in time and formed by the delayed multivibrator ZMsg carries out the following simultan- 
eous operations through a hole in the perforated card PC. 


1, The initial values x, and yg are transferred from the card, through gates P,, to the integrator registers 
Ry and Re. 


2, The total displacements along the x and y axes, Ax and Ay, recorded on the perforated card by a supple 
mentary code, are introduced into the displacement meters DM-1 and DM-2, 


3, The nature of the programmed curve is established (straight line or circle), For this the curve form 
trigger C, controlled through opening 36 or 37 in the perforated card, opens or closes gate P35. 


4, When a circle is programmed the initial state of the sign trigger S.Tr. is set up through openings 1 or 
2 of the punched card and gates P", 


5. A link is established for the transition of +x or —x and +y or —y pulses from the programmer to the 
milling machine's reversible counter [2]. (Openings 46 or 47, 48 or 49 in the punched card, gates P, and triggen 
L or K). 


Subsequently trigger B cuts off gates P,, Pz, and P,;, Frequency F is connected to the frequency divider FD 
by means of triggers D, M, and N, gate P, and the gate controlled by them P,. Gates P, connect the x and y out- 
puts of the programming device to the integrators. © 


The use of delayed multivibrators ZM;, ZMg, ZMg, and ZM, in the circuit is connected with the necessity 
of establishing small time intervals between the control pulses. The time delay is basically determined by the 
resolution of the triggers r. The value of r lies between 5 and 10 wsec. Thus the delay created by the indica- 
ted multivibrators must be somewhat greater than r. In this way, the over-all time for preparing the 
programming apparatus for operation is of the order of 40-60 psec. 


Programming of Straight line segment. To achieve linear interpolation,gates P,; open the feedback loop, 


The angle of inclination of the straight line is a = arctan 2, The values of xq and yo, entered in registers Ry 


and Rg, and the frequency F must be such as to ensure the um rate of displacement of the instrument, The 
coordinates of the end of the straight line segment are given by the values of Ax and Ay which are introduced 
into the displacement meters. 


Programming of circles, For the programming of circles, gates Ps close the feedback loop, The radius of 


the circle is R=) x3 + y?. 


Because the origin of coordinates is shifted every 
time to the starting point of the programmed circle, 
the position of the center of the circle is given by the 
coordinates x and y, and the initial state of the sign 
trigger, The direction of the instrument's motion Is 
determined by the +x or —x and +y or ~y channels, 


4 The lar frequency is w = F /2"**, In order 
to ensure linear velocity V = wR the pulse 
repetition frequency F must have a particular value. 
If it is necessary to program circles with both large 
and very small radii, then it is necessary to provide 
for the possibility of changing the frequency F over the required range, The value of the frequency F must be 
either produced automatically or recorded earlier on the punched card, However, it is possible, without changing 
the maximum value of the frequency P, to redesign the circuit so as to provide for the possibility of decreasing 
the number of working trigger units in the registers and the frequency divider for the programming of circles with 
small radii, 


ber of 


the 
case 
servi 
pulse 
(gate 
The 
of Z] 
zero, 
one 
inclu 
150 r 
as, f 
consi: 
DM-% 
vibra: 
5 cry ste 
a proj 
deterr 
114 


The error can be reduced, as it was mentioned above, either by the reduction of the size of one step or by 
the decrease of the rate of change of the integrands, In the first case it is necessary to increase the number of 
columns in the registers and the frequency divider while leaving the absolute value of R unchanged, In the second 
case, it is necessary to increase the number of columns in the registers only, At the same time the first m cells 

in the registers are not connected to the Pcircuits. In this case w = F(2"**2™yFig, 6). 


In order not to complicate the feedback unit, the number ofpilses fed from the integrator to the discrete 
servosystem of the milling machine must be decreased, By the simultaneous choice of size of the step and de- 
crease in w it is possible to reduce the integrator's error to less than one step on the milling cmaine, 


The circuit operates until the time when the displacement meters DM-1 and DM-2 overflow. Overflow 
pulse 50 or 38 (overflow does not, as a rule, occur simultaneously) disconnects the corresponding output channel 
(gate P,) by means of triggers N or M, and disconnects frequency F (gates Ps, Py) at the completion of a process, 
The system is ready to process a new profile, The starting pulse can be initiated either automatically by means 
of ZMs; or manually by key T. The punched card must be transferred when a given profile is to be processed. 


In the programming of circles it is necessary to reverse the register counters when the quantities x or y reach 
zero, and to switch over the +x or —x and +y or ~y channels when the quantities x or y attain their maximum 
values, Since the maximum value of x or y is not exactly fixed in time the reversing is carried out only when 
one of the coordinates passes through zero, The last is achieved by means of voltage coincidence circuits [3], 
including the resistors R, connected to the register triggers, the NE circuit and the former ®. 


The circuit of the programming apparatus was designed for the reproduction of circles with radii of up to 
150 mm, and straight lines with Ax and Ay up to 600 mm, with a step size of 0,01 mm, 


The circuit is sufficiently simple in comparison to the circuits of the other second-order interpolators, such 
as, for example, the squarer, etc. 


The circuit of the programming apparatus, which was built at the Electromechanical Institute of AN SSSR, 
consists of six units. In each unit are assembled separately registers R, and Rg, displacement meters DM~1 and 
DM-2, frequency divider FD and the control circuit, including the auxiliary triggers A, B, etc., the delayed multi- 
vibrators and the P circuits. An over-all view of the register unit is shown in Fig. 7. The circuit contains 323 
crystal diodes and 325 type P, junction transistors. 


APPENDIX 


As a result of an experimental verification of the accuracy with which circles are reproduced by means of 


a programming device based on digital integrators with successive carriers, the following maximum errors A were 
determined: 


1) without auxiliary units (m = 0) A = 27 units; 
2) with m = 1, A = 10 units; 
3) with m = 3, A = 3 units, 


Subsequent increases in the number of auxiliary cells results in a relatively slow decrease in the error, 


Errors were determined for a large number of circles with various values of x, and yy) by counting the num- 
ber of output pulses in the x and y channels, 


Fig. 7. 
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Radii of circles were chosen in the range of 12 to 8191 pulses. The quality of the surface of a profile (circle) 
worked on a milling machine controlled by the given programming apparatus is considerably better than that ob- 
tained when the machine {s controlled by programming apparatus with linear interpolation. These checks were 
made on an LR 6441B machine, exhibited at the All-Union Industrial Exhibition in 1958, 
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CHOICE OF AN OPTIMUM AMPLIFICATION FACTOR FOR A 
SELF-ALIGNING CONTROL SYSTEM 


Perel'man 


(Moscow) 


A method is given for determining the optimum (from the point of view 
of reducing the mean=square control error) amplification factor in a control sys~ 
tem with a self-aligning program, when a cyclic disturbance at its input is de- 
scribed statistically, 


A control system designed on the principle of a self-aligning program is described in [1]. There it is also 
shown how this principle can be applied to the solution of the problem of controlling the thickness of the finished 
product in hot~-rolling, strip mills, 


This paper is devoted to the determination of an optimum amplification factor for this system for the case 
when a cyclic disturbance f(t), described statistically, appears at its input. Just as in [1j, the problem of thick- 
ness control in a hot-rolled sheet is used as an illustration, 


The investigation is restricted to the case when the system under investigation (Fig. 1) possesses a low inertia 
regulator 1 and a controlled member 2, This implies the situation when transients arising in the regulator and the 
member as a result of their being excited with a unit 
step function, vanish in a time less than T, where T is 
the length of the interval of integration [1], Despite 
such a limitation, the problem under consideration does 
not lose its practical importance, Thus, for example, 

a system for thickness control in hot-rolling can be 
looked upon as a low inertia one, because satisfactory 


(a7 results can be obtained by dividing the strip along the 


at) whole of its length (with the exception of the "tail" 

section, which is about 10-15 m long) into 25-30 m 

Fig. 1 intervals [2, 3] which corresponds, if one considers 
—~ existing rolling speeds, to intervals of integration 4-5 


sec long, At the same time transients in the drives of the pressure devices can die down in tenths of a second, 


Determination of the Mean-Square Deviation of the Output Coordinate 


Let us suppose that a disturbance f(t) was applied to the member at a time t = 0 coinciding with the be- 
ginning of the course of work of some cycle of f(t), Denoting this cycle by the number 1, we renumber in order 
all the cycles following it. Denoting the times at which the work cycles begin by tos, tog, etc., we can replace 
the function f(t) by a series of functions 


h(t —ta) = A(t), 
(t — toe ) = fa (t*), 


“fn (t —ton) = fn(t*), 
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where each function f ; (t*) coincides with the function f(t) between the limits of the jth cycle and at all other 
times is identically equal to zero, 


In an analogous manner, a change in the output coordinate x(t), which is dependent on the application of 
the disturbance f(t), can be represented as a series of functions x,(t* ), xg(t*), xs(t*), .. .. X,(t*). From the 
definition of a cyclic function [1] the condition 


(t*) = (t*). (1) inf 


is imposed on the behavior of f ;(t*). 


Condition (1) can be satisfied for a limited number of cycles (Nj) of the disturbance. We shall call the to- 
tality of these cycles a series, When a given series of cycles ends, a second series of N, cycles can begin, then 
a third ofN, cycles, etc, In each series condition (1) must be satisfied, but in the transition from the last cycle 
of one series to the first cycle of the next series, the function f(t *) can change its form in a jump, tf.e., condi- 
tion (1) can be violated, 


The numbering of the cycles of each series will be carried out independently of the other series, 1. e., the 
first cycle of a given series will be numbered 1, etc. A set of strips of hot-rolled steel sheet can be used as an 
example of a serles, In this case, as it is shown in [1, 2], within the limits of the given series f ,(t*) can be re- 


presented by 
fn(t*) = f, (t*) + fon(t*), (2) 


where f,(t*) is the regular constant (for a given series) component function, f op (t*) is the random function of 
the arguments t* and n, determining the deviation of f,(t*) from f;(t*). 


As a consequence of the definitionof a cyclic function, the duration of the intervals of integration can be 
chosen in such a way that the continuous functions f j(t*) and xj (t*) can be replaced with sufficient accuracy 
by the step functions f, (i) and x (1), where i is the current number of the interval of integration, while f ; (1) 
and x; (i) are equal to mean values of f jct* ) and xy (t* ) within the limits of the ith interval of integration, 
Here, {= 1, 2,...,M, where M is the number of intervals (the number of integrating links) in a given system, 


By introducing into the analysis the functions f (i) and xj (1) and under the condition that the control system 
is a low inertia one, the control processes on intervals having different numbers, for example, on the ith and com 
(1 + 1)th, turn out to be practically independent of each other. ; func 


Taking into account the assumptions noted above, as well as the fact that the duration of the interval of 
integration T is always chosen to be less than the delay time r in the member, we can represent the system under 
consideration (Fig. 1) as a set of M not interconnected systems, each of which (Fig, 2) maintains control over a 

definite (ith) interval of integration (switch Kj closed when (i— 1)T = 


case 

F =< t*< 1T), At the time, the transfer functions of the regulator, mem= omg 
ber, and integrator in the system of Fig. 2 are represented by three real 

i-th numbers Kreg, Kme, and Kip. 
K, : Let us examine the control process taking place in the system of 


Fig. 2 chosen for some fixed interval of integration 1. Since i is chosen 


Fig. 2. to be fixed, this symbol will be dropped henceforth, 1. e., f ; (i) and 
8 x; (1) will be denoted as fj and xj, From an examination of Fig. 2, where 
one can deduce a recurrent relationship for determining the deviation of the quantity being stabilized: and x 


Xn = Kme(fn~ Kreg n), (3) 


where ¢,, is the value of the (program) signal on the nth cycle, which is determined by the sum of the errors aris- 
ing from preceding cycles of the given series, and the preliminary program ¢» accumulated in the integrator be- 
fore the commencement of operation with the given series: 


*Ab 
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|| 


n—1 
Kind) (n>2). 


Kme ¢an be taken as unity without reducing the generality of subsequent discussions, 


Assuming that before the beginning of each series the integrator is cleared of all previously accumulated 
information, 1, e., that g) = 0, and introducing the notation K = KinKreg: we obtain 


K 2}. (3a) 


Let us expand expressions (3a): 


n—1 n—2 n—2 
aed j=1 j=1 j=1 


n—2 
— K K(1—K) = — K (1 — K) x 
j=1 j= 


n—3 n—3 
x (see + » 25) = 
j=l j=1 


= fn— Kfn—1 —K(1 — K)fn—g —K (1 — K)* 2) =... ete. 


Finally 


In= fn —Kfn — K (1—K) In—s —K(i—K)*fn—s—...— K(i—K) fy. (3b) 


Let us consider the disturbance to be a function of the argument n: fy, = f(m). On the basis of (2), f(n) is 


composed of two components — a constant, for the given series, component f, and a component described by the 
function A f 


= fot Afn. (2a) 


The value of Af, averaged over all cycles of the series must be equal to zero, insofar as in the opposite 
case it would be included as part of fy). From this it is easy to see that in the case when A fy is a stationary func- 
tion of the argument n, its mathematical expectation* A f , is equal to zero. 


When (2a) is satisfied, the deviation of the output coordinate can also be represented as a sum of two func~ 
tions of the argument n: 


Ln = = a(n) + %(n), 


where x,(n) is the regular component of x(n), dependent on the presence of the component f, in the disturbance, 
and x, (n) is a random component dependent on the presence of the component A fp in the disturbance. 


Substituting (2a) into (3b) we obtain 


x(n) = x, (n) + = — K —K (1 —K)—K (i —K)*—-- 
K (1 — + Af, — — K (1 — K) — 


* A bar above a function will denote its mathematical expectation, 


— K (1 — = fo (1 — + Afn — KAfn—1 — 


— K (1 — K) Afn—, —-+-— K (4 — (4) 


Examining (4) we see that the condition for stability in a system with an inertialess regulator and a con- 
trolled member is the satisfaction of the inequality |1—K|< 1, l.e., K< 2. 


In the case when the system 1s stable, limx,(n) = 0 and lim x(n) = x,(m), Let us find the value of the square 
n> 0 n> @ 
of x(n), averaged over the set (ensemble) of the series: 


= Tz, + = + 22(n) (n) = f(t — + 
+ K (1 — K) — K (1 — + 
+ 2(1 — fo — — K K(i — 


Assuming that there is no correlation between the random and regular components of the disturbance, i, e,, 
foaf,, = 0 for any n, we obtain 


= Ba — + 


+ Aff + + — KPAP_, + — 


— —2K(i—K) — 2K (1 — -—2K(1—K)" + 
+ 2K? (1 — K) Af, + 2K? (1 — AT, + 
+ 2K* (1 — KP (1 — + 
+ 2K* (1 — + 2K? (1 — + 
+ 2K? (1 — +++ 2K2(1 — K)" AF, + 
+ 2K* (1 — + 2K? (1 — + 
+ 2K? (1 — +++ ++ (1 — + 


the | 


+ 2K? (1 — + 2K? (1 — + 
+ 2K2(1 — 


We shall assume that A f, is a stationary, random function of the argument n. In this case 
Af, =R(0)=D, 


where R(m) is the correlation function and D is the dispersion, 
Taking (6) into account, we transform (5) and obtain 


= fo (1 — + 
+D+K°D KP + (1 4+ ---4+(1— Kp) — 
— 2KR(1){i —K (A—K) — K K (1 — - -— K(1— K) — 
—2(1— K) KR(2)(1 —K (1 — K)— K(i— K)®—- --— K(i— 


(6) 
— trolle 
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—2(1 — K)" KR (n—1). 


Let us represent (7) in the form 


n—1 
= — + B,D +S) R(m)B(m). 
m=] 


From (7) we determine 6,» and 6,): 


Bo = 1+ K®(1 + (i — K)? + (1 — (1 — = 
K?/1 } 
i—(i1—Ky 


+(1— 
B(m) = — (1i—K)™—1 + (4 — 


Making use of (8) and (9), we have 


z(n) = — 4 — (1 — 


—K Rim) («a — — (1 — }- (10) 


In this way, making use of (10), we can determine the mean-square deviation of the output coordinate of 
the member x*(n) from the given fy, R(m), and K, 


For example, let the correlation function have an exponential form: 


R(m) = D=R(0) =a. 


Substituting (1) into (10) and carrying out a summation we obtain 


1—e~® (1— K) i—K—e* 


Determination of the Optimum Amplification Factor 


The problem being solved by the control system under examination consists of the stabilization of the con- 
trolled member's output coordinate; 1, e., the reduction of the value of x(n). 


In that case when the random component of the disturbance is absent, i, e., f(m) = fy, the value of x(n) 
determined from (4), will be equal to 


x(n) = 2, (n) = f(t — 


—2(1 — K)?KR A(i — K) — K (1 — K)® —---—K (1 — — 
(7) 
B(m) = — 2K (4 — (1 — K(1 — K) (1+ (1 — KP? + (1 — 
(11) 


At the same time it is obvious that the optimum value of K from the point of view of the fastest decrease 
of x(n) will be K=1. In this case x(n) = 0, when n > 1, 


However, when a random component A f,, is present in the disturbance, x(n) is also a random quantity, the 
mean~square value of which (12) is a function of n and K and cannot be reduced to zero for any value of K, 


In this case, to determine the optimum value of K, it is necessary first of all to determine such a function 
of x(n), the value of which would be the best criterion for the quality of the system's operation, 


Let us denote the function by #(n, K), The coefficient K, ensuring the most profitable extreme value of 
¥(n, K), will, obviously, be the optimum coefficient. 


We shall proceed from the assumption that the number of cycles N during which the given disturbance acts, 
is known as a real, initially given quantity or as an averaged statistical quantity. 


Thus, in the case of the control of the thickness of a hot-rolled sheet, the quantity N, corresponding to the 
number of strips in a set being rolled, can be given before the rolling of the given set as the true number of strips 
in the particular set, or as the mean (for the particular machine under the given technological conditions) number 
of strips in the set being rolled. 


For a known N, the function ~(n, K) can be determined as 


N 
K) = 9(N, K) = >) FQ). (13) 
n=) 


Let us note that the specified requirements imposed on the quantity x(n) in a concrete case, can result in 
(13) not being the best evaluation, Thus, for example, if the value of x(n) exceeds some set tolerance xX» and 
results in the production of a completely defective product, while for | x(n)| « x», then (13) will not be the best 
evaluation of the quality of control, 


Consequently, it is expedient to introduce a more general evaluation 


N 


K) = (n) (14) 


n=1 
where a = a(N, ff, a) is an empirical coefficient chosen from technological and technoeconomic considerations, 
Evaluation (13) is obtained from (14) when a = 1. 
Substituting (12) into (14), we obtain ° 


N 


n=l 


(1 (e~* + 1—K) (4 —(1— 
2aK 2(1—K —e—®)(2— K) 


Ke~®(4 — e—®)(2— K)[1 — (1 — 
(14 —e — K — 


(15) 


For given values of N, ft, B and a, the quantity Kopt can be found from the condition for a minimum in 
expression (15), 


The determination of Kop; from (15), which represents in general, a very laborious problem, fs simplified 
to a considerable degree when the random component of the disturbance is white noise. 
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The Case of White Noise 


Let us examine the case when Af is white noise, i, e., R(m) = 0(m=1,2,...). At thesame time (15) 


+ + 2%(n)] = 


where 


8 (16) 


To investigate (16), let us prove that Kopt = 1, Let us rewrite expression (12) for the case of white noise, 
= ft (1 — + a7) 


Insofar as from the condition for stability 0 < K =< 2, both terms of (17) are always positive. At the same 
time, x2(n) grows monotonically with the growth of (1— K)*, while x%(n) grows monotonically as K changes from 
_ 0to 2, From here, it is easy to see that the minimum in (17) for any n, and the minimum in (16) for any N(N > 2) 
can lie, if stability conditions are satisfied, only in the region 0 =< K <1, 
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Fig. 3. 


In this way, when investigating (16) for a minimum, one can assume that 2—K > 1 and O=1-K <1, 
and, consequently, for a sufficiently large N , (16) can be represented with any degree of accuracy as 


N 
+ (n) + a(n)|= W (18) 
n=1 


The value of K which ensures a minimum in (18) will be 
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A graph of Kope as function of is shown in Fig, 3, 


As K -—> 0, the second term in expression (16) 


8 
V. 
tends toward an indeterminate quantity of the form 0/0, and its replacement by the expression wae be- 
comes incompetent, 


Thus, to ensure that the Kop; found for the approximate equation (17) {s also the optimum one for equation 
(16), it is necessary to show that (16) has one minimum in the range 0 <K <1 for sufficiently large N. Let us 
prove this. 


Taking into account that when 0 < K <1 the term [1 — (1— K)*§}/[2(2— K)N] tends toward zero more 
rapidly than does */,N, for a sufficiently large N and for any 0 < K <1, equation (16) can be written as 


N 
= 20 —K (20) 
n=1 
Differentiating (20) with respect to K we obtain on 
inc 

N 
1. 

inf 
in 
The first term in (21) represents, for changes in K for 0-1, a positive monotonically increasing function; - 


the second term, under the same conditions, represents a negative function with a monotciically decreasing 
modulus, 


It is obvious that under these conditions, the equation 


N 
N 
2 1—K)*"""' =0 
can have not more than one root in the interval of real yalues of K between 0 and 1, which is what was required 


We must add that the conditions for (16), to be a good approximation of (18), almost always exist. 
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ANALYSIS OF THE COMPENSATION FOR STRIP THICKNESS 
VARIATIONS BY MEANS OF AN ELECTRONIC ANALOG 


S. A. Doganovskii and A. A, Fel*dbaum 
(Moscow) 


The paper deals with the theory, calculations, and the possibilities of 
applying computer circuits for compensating strip thickness variations in a roll- 
ing mill, Test results of the strip thickness variations compensation obtained 
on an electronic analog of a rolling mill are quoted. 


The development and perfection of modern continuous cold-rolling mills is to a great extent the result of 
achievements in the design of their electric drives and the automation of their control [1-3], The change-over 
to complex automation of rolling mills requires a further study of a number of questions connected with the auto- 
matic control of the strip thickness in rolling. The solution of these problems will ensure production of strips with 
minimum tolerances considerably smaller than those specified by the existing GOST (State Standard), and will 
increase the productivity of the rolling mills, 


Various ways of using computers for the achievement of these aims are being considered at present, 


1, Development of a Compensating Computer 


Figure 1 shows the principle of operation of a continuous cold-rolling mill; below we give only as much 
information about the mill as will be required for understanding subsequent statements, Those who are interested 
in studying the principles and theory of the rolling mill technique should read the specialized literature (for in- 
stance, [4-7]). 


Figure 1 shows diagrammatically, not to seale, four stands and the strip-rolling process, 
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A strip of thickness Hj passes through the stands of the mill, {s rolled down in each stand and leaves it with 
a decreased thickness, The main reason for the inaccuracy in preserving the set dimensions of a rolled strip is 
the variation of the stress on the rolls of the mill, which leads to variations in the stand strains and hence spaces 
between the rolls, 


Thus the thickness H, of the rolled metal at the output of the mill varies due to the operating conditions 
(heating of rolls, changes in lubrication, etc.), properties of the rolled metal and especially to the changes in the 
thickness H; of the strip at the input, the so-called input thickness, Variations in the input thickness are usually 
due toinaccuracies in the operation of the hot-rolling mill which produces the strips for cold rolling. Improye- 


ment of the accuracy of hot-rolling is a separate, very complicated problem which will not be dealt with tn this 
article, 


The variations in the input thickness have a large amplitude and low frequency and can be reduced to a 
minimum in the first stand by means of a thickness regulator operating the pressure mechanism PM, of the stand 
(Fig. 1). 


A suspension micrometer My, which measures strip thickness H,, is placed at distance I, behind the first 
stand, At present contactless measuring devices are coming into use, which measure the absorption of x-ray or 
radioactive radiations passing through the rolled metal [2]. 


If the rolled metal thickness Hy, is greater than required, the electrical pressure mechanism operates, de- 
creases the spacing between the rolls and reduces the thickness to the required value, If the thickness is smaller 
than its nominal value, the spacing between the rollers is increased. 


Due to an insensitive region, inertia of the pressure mechanism, lag in measuring and slow displacement of 
the rolls, this control system operates slowly, It can only eliminate thickness variations which ‘ave a sufficiently 
low repetition frequency, The use of a hydraulic drive would considerably speed up the operation of the pressure 
mechanism, but it would require large structural alterations in the existing rolling mills, 


Strip thickress variations of higher repetition frequencies and smaller amplitudes, remaining after the pres 
sure mechanism regulation, can be compensated by controlling, for instance, the tension between the second and 
third stands, If the forward tension F, is increased at the second stand, 1, e,, at cross section AA' (Fig. 1), the 
thickness of the strip leaving that stand will decrease, This is due to a decreased pressure on the rolls resulting 
in a smaller strain in them and their associated mechanisms, which leads to a smaller spacing between the rolls, 


Tension F, can be controlled by changing the angular velocity Q, of the third stand rolls, When Qy is in- 
creased, the strip enters the third stand at a higher velocity. The strip between the two stands can be considered 
as an elastic spring whose stretching increases tension F,, With a decreasing 9, tension F; also decreases, Thus 


by appropriately changing the driving speed of the third stand rolls, it becomes possible to compensate for the 
strip thickness variations so that at the output of the mill, thickness H, remains constant [8-10]. 


Since the tensions between the stands, the torques, etc, must not exceed certain limits, it is possible to com- 
pensate by this method only small thickness variations, but not the large ones (see below). 


In order to change the third stand driving speed appropriately, it is necessary to measure the strip thickness 
by micrometer My, before the given section approaches the second stand and feed the information to the compel 
sating computer (abbreviated CC) marked CC;. This computer calculates the value of the correction y to be 
added to the set speed Qgp of the rolls, In calculating the correction, the following factors should be kept in mint 


a) the time lag of the metal cross section passing from the place where it is measured by micrometer My 
to the second stand, This lag is beneficial since it provides time for the stand drive to be “prepared” for the 
changed conditions. 


b) the dynamics of all the stand drives, Since all the stands are interconnected by the elastic strip, the 
dynamics of all the drives influence the compensation to a certain extent, 


c) relationships characteristic of the stands, for instance, between the moment of the drive loading and the 
front and back tensions, also the thickness of the strip at the input and the output, etc. 


d) elasticity of the strip between the stands. 
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e) limitations imposed on the variable quantities characterizing the operation of the drive and stands. 
These variables can only change within limits permissible for them. 


With ideal compensation, the thickness at the output of the mill should be strictly constant, In practice, 
however, owing to certain additional factors which are not taken into account such as errors of measurement, 
changes in the parameters of the stands and drives, variations in the properties of the rolled metal, etc., and also 
to limitations imposed on the values of variables, the thickness at the output cannot be strictly constant, In an 
actual compensation system, device CC, can only decrease the thickness variations at the output, The thickness 
variations can be decreased to a greater extent if the thickness is measured by micrometer Mg at the output of 
the mill and the information is supplied to a second computer marked on Fig, 1 as CC, The second computer 
CC, affects CC, in such a way as to decrease the thickness variations, Such an arrangement should be classified 
as a self-regulating system and will not be considered in this article. 


2, Development and Design of a Compensating Computer Circuit 


For designing a CC, the relationships noted above must be known, The most important of them are the equa- 
tions of the main roll drives and those of the stands, In their original form these equations were obtained by the 


TsKBMM TsNIITMASh (also see [11]). In order to reduce the stand equations to a form convenient for simulation 
and analysis, let us adopt the following notations: 


a) Let all the quantities belonging to the ith stand have subscripts i. 


b) Let all the quantities qualifying the state of metal before it enters a stand be marked by a dash; those 
occuring on leaving a stand without a dash, 


c) Let us write all the quantities with capitals (Fig. 2). For instance, Let H be the thickness, F the tension, 


V the speed of the strip motion, Mg the moment of loading of the roll drives, and N the number of roll revolutions 
per minute, 


d) Let all the quantities corresponding to a certain nominal conditions have the subsscript zero, Let us de- 
note deviations from the nominal condition by small letters, Thus: 


H=H,+h, V=Vy+v, F=Fo+f; My +m N= N,+n7. (1) 


Equations for the ith stand can be written by means of the above notation in the form: 


My, = hy ky, — hihi, — + fia 
(Vio + v1) (Hig + hi) = (Vio + 04) (Hin + hu), (3) 


(4) 
= By (fi — fir) (Neo + mi) + Dini. 


hy = — — + (6) 
hy (t) = his (t — 1), (1) | 

é (8) 


Let us explain the meaning of these equations, 


Equation (2) represents the moment of loading, It shows that the moment of loading depends on the forward 
and back tensions and on the thickness variations both in front of and behind the stand, This relationship is in gen- 
eral, nonlinear and can be represented by expression 


M a='Mq(Hi, Hi, Fi-1, (9) 
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If the deviations from the normal condition, however, are considered to be small, it is possible to use a 
formula with differentials 


aM, aM, 


and to replace the differentials by small deviations, By replacing coefficients °8M4/0Hj, 9Mg/0Hj,dMq/3Fj-s, 
and 0M, /dF; , which are constant for a given nominal condition (but vary for different conditions), by kj,, kis, 
kgs and k obtain formula (2). 


Equation (3) represents continuity. It shows that the volume of metal entering a stand in unit time is equal 
to the volume leaving it in the same time. This relationship is nonlinear. 


. Equation (4) represents the elastic strip between the stands. If the elastic strip of metal between the stands 
is considered as a spring, the tensile force Fj will be proportional to the difference in the coordinates of the two 
extreme cross sections at the beginning and end of the spring, with the nominal values for which the tension ts 
equal to zero taken as the basis of reckoning: 


Fy = ky — Xi), (10) 
where k; = const, By differentiating this expression and substituting the values of 


ete, 


we arrive at formula (4). 


Equation (5) represents the forward flow law. The speed of the strip at the output is usually somewhat higher 
than the angular speed of the roll, This relationship is also nonlinear; the nonlinear term in this expression, how- 
ever, is not important. 


Equation (6) represents the law resulting from the pressure on the rolls and their strain. It shows that the 
strip thickness variations on leaving the stand depend on the thickness variations at the input of the stand, the 
forward and back tensions in the strip and the speed of the strip, Variations in the strip tension during rolling 
change the direction of plastic strains in the flowing region of the metal, resulting in changes of the specific pres 
sure on the rolls and of the spacing between them. It has also been shown experimentally, that changes in the 
speed of rolling lead to variations in the friction coefficient and hence, the thickness of the rolled strip. 


Equation (7) accounts for the time lag rj. in the strip movement,due to its finite speed between the stands, 
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Finally, Eq, (8) determines the value of this lag. This equation only holds with Vj = const. If, however, 
V; = Vj (t), Eq. (7) must be replaced by a more complex relationship: 


- 


\ Vi (t) dt. (11) 


Quantity 1; is determined from Eq. (11) with given values of 2; and Vj (7%), Since this law is complicated, 
and with a slowly changing Vj(t)-is transformed into Eq. (8), we shail use henceforth formula (8), 


The above statements represent part of the work carried out in developing an electronic analog of a rolling 
mill, Equations of the main roll drives of a stand,* produced by the TsKB MM, were simulated with the aim of 
obtaining simplified forms of equations, They are not quoted here in view of the limited size of the article, An 
oscillogram of the transient function of the drive speed control closed system, simulated electronically, showed 
that with an accuracy sufficient for practical purposes,it is possible to substitute this closed system in a linear 
approximation by a series connection of a time delay and lagging elements (see Fig. 3a, where y is the speed 
regulator setting correction). 


If the effect of the moment of loading on the drive speed is taken 


o% og into consideration, the circuit in Fig. 3a should be replaced by that in 
ey y : 7 Fig. 3b. Let us denote the original quantities by letters without a bar 
> Pr and the simulated ones by letters with a baron top of them, Then equa- 


tions of the circuit in Fig. 3b will become 


kyy’ — kmmg = (1 + pT)n, 
yer. 


(12) 


Expression (12) represents a simplified drive equation where quan- 
Fig. 3 tities rT») and T are determined from oscillograms of the drive transient 
> function, and constants ky and ky, from the static drive equations, 


Let us first calculate the basic characteristics of the computer. 
For this purpose, let us further simplify the stand equations (1)-(8) keep- 
ing the simplified equation (12) for the drives,** In simplifying the 
stand equations let us make the following assumption (Fig. 4). 


1. Let us consider the drive of the second stand to be ideal (with- 
out inertia) and the speed at the output of the stand to be constant, We 
shall neglect the forward flow at the second stand, Then: 


ic> Drs = 0. (13) 
(p) 


Fig. 4. 


On the basis of condition (13) we obtain from equation (4): 
Phi = (14) 


2, Let us also consider that the increment in the metal thickness hj at the output of the second stand de- 
pends only on hi and fi: 


hy = high, — (15) 


3, Let us consider the forward flow time r between the instant a given strip section is measured by microm- 
eter M, and the instant it arrives at the second stand to be a constant, Then the following expression will hold: 


— 
* The circuit of this drive is described, for instance, in [1]. 


** A check by means of an electric analog (see below) showed that such an approach to the problem was per- 
missible. 
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hy = (16) 


4, Let us assume that the object of compensation {s to provide a constant thickness at the output of the 
second stand, {, e., to fulfill the condition hj = hj,,;= 0, In fact this condition is not quite accurate, since even 
if hj,s = 0, quantity hj,; may not be equal to zero because the back tension at the third stand {s also changed, 
It would be more correct to assume that hy,, = 0, but this complicates the conclusions, With hj = hy,) = 0, the 
continuity equation for the third stand assumes the form: 


H i+1,0° Vita = (17) 


5. The speed at the exit ofthe third stand is determined from formula (5). By linearizing this formula and 
assuming that ¥;,; = 0, 1. e., that the forward tension at the third stand {s constant, we find: 


= — Bip ft + (18) 


Finally, the expression for the moment of loading mg ;,, 1s simplified owing to the assumptions made 


= = Figs = 0): 


ita = ky, fi. (19) 


Equations (13)-(19) are used below for calculating the computer circuit. 
The compensating device formula can be written as 


y= Ko (P) (20) 


and the operator Ky(p) can be found from the system of equations (13)-(19), if we assume that hj =0 (condition 
of compensation). Omitting intermediate operations, let us write the compensating computer transfer function; 


Ko(p) = (agp? + ap + ap). (21) 


The coefficient here represent the following quantities: 


Tk,, ky, (Thur + 4) 


Nits 


The amplitude-phase characteristic of the compensating computer is represented by the expression 


Ko (jo) = [ag (jw)? + ajo + ao). 


The amplitude- frequency characteristic is calculated from the formula: 


A(w) = | Ko (jw) | = V + (40 — 


and the phase characteristic has the form 


@ (m) = — w (t — to) + are tg 
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Let us find the limiting amplitude of fluctuation hy., which can be ideally compensated, taking the stand 
possibilities into account, Let there be limitations imposed on amplitudes of fj, y and on dnj,,/dt; quantity 
hy. varies sinusoidally, The disturbing frequencies hj., of the input thickness is usually limited to a range of 


2 cps, mainly due to the eccentricity of the hot-rolling mill rolls, Let us express quantity hy., in the form of a 
rotating vector: 


hy_, = Heit, (26) 


From Eq. (15) we find for hy = 0 


k k 
= hy = othy_, = 
kig kis 


kis =n 


dn. 
Similarly, formulas for vectors my41 (1) and y can be found, The moduli of these vectors give ampli- 
tudes corresponding to the fluctuations: 


max 


ly| = HV (a0)? + (a — < 


k 
= <f 


Thus, the limiting conditions have the form: 


H< max 
V (ayo)* + (a9 — 


Variations of the input thickness with an amplitude larger than the permissible one, 1, e., which does not 


satisfy one of these inequalities, cannot be compensated by tensioning the strip. For this purpose the pressure me- 
chanism regulator should be used. 


K, (p) j 

Fig. 5. 


Let us now design a block schematic for a compensating computer with the help of formula (21). The ele- 
ment with operator e~P(T-To) can be obtained in various ways. It is possible to use, for instance, a controlled 
delay unit, which has rotating brushes sliding over segments to which condensers are connected, The voltage 
Value proportional to h, is first "remembered" by the condenser and then, after an interval of time 1, it is mea~ 
sured by another brush (see, for instance, [12}). 
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For correct operation of such a delay unit, the speed of rotation of the brushes (v,,,+) must be synchronized 
with that of rolls or the strip (v,9}). 


If the speed of variations h does not exceed 2 cps and the delay time r is no greater than 1 sec, it is quite 
sufficient to have 10-15 condensers, A similar method with variations h recorded on a magnetic tape can be used, 


In order to account for the effect of rp it is necessary to change the speed of rotation V,o¢ according to a 
more complicated law. Let 


(30) 


where L is the distance between the recording and measuring brushes, 


Taking into consideration that r = 1/v,.), where I is the distance between the suspension micrometer and 
the second stand, we find v,..+: 


L 
Vrot = Yrol (31) 


Thus the block schematic representing operator Ky(p), will take the form of Fig. 5, Here RDU stands for 
regulated delay unit and NT for nonlinear transducer, At the circuit output, we obtain y, the correction of the 
set drive speed of the third stand. Other methods of construction are also possible. 


Thus a theoretical design of a sufficiently simplified circuit shows the basic features of a compensating 
computer block schematic. 


3. Testing of a CC Circuit by Means of a Rolling Mill Electronic Analog 


Formula (21) was derived on the basis of a simplified circuit analysis, Further work consisted in the construc- 
tion of an analog of a three-stand rolling mill. In this connection, the full system of technological stand equa- 
tions (1)-(8) (the number of equations amounted to 22 for all the three stands) was preserved, and each main wll 
drive equation was made to correspond to expression (12). 
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When the stand equations were simulated, it was assumed that deviations f, and f, remain equal to zero. 
Further, the continuity equations were replaced by linear equations which in this case assume the form: 


Viohs + = Vig ht + Hor (32) 


A full system of equations for the three stands falls into three groups of equations similar to each other, 
This can also be seen from the block schematic of the electronic analog consisting of separate units simulating 
each stand (Fig. 6). Groups of units are interconnected by means of circuits containing constant time delay units 
CDU-1 and CDU-~2 simulating type (7) equations and circuits simulating elastic strip (4) and continuity (3) equa- 
tions. 

In the first stage, the work was carried out with a constant strip speed analog. Hence, constant delay units 


were used, based on an approximation of function e~P™ by means of a series of n terms (Fig. 7), where each term 
is a ratio of two polynomials of p. 
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The block schematic of the mai™ roll arive has been previously analyzed (Fig. 3b). Units Dr, and Dry 
(Fig. 6) have no compensating inputs and thus represent typical lagging elements whose inputs are fed with quan- 
tity mgj. Unit Drs has also an input for the correction (y) of the third stand drive regulator setting, and there-~ 

fore has a separate block schematic given in Fig. 7. The delay element CDU-3, 
h | | contained in the Dr; unit of the rolling mill electronic analog, is constructed in 
— | hy a different form from CDU-1 and CDU-2 so as to obtain a better zero stability 
ae | > of its output voltage which is important for the regulator setting correction cir- 
g cuit, The principle of operation of this element consists in representing the func- 
Fig, 8. tion e“P™ by means of a Pada series [13]. 

In order to obtain the required parameters of the compensating computer, it is advisable to examine the 
electronic analog of a rolling mill in the form of the system shown in Fig. 8, This system has two inputs h,, and 
y and one output hs. Two quantities are fed to the input of the system, hy , the quantity measured by micrometer 
M, and y, the correction for the third stand speed regulator setting. The output quantity is hg, the strip thickness 
at the output of the third stand. 
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Fig. 9. 


Let us examine the frequency characteristics of the system for each of its inputs taking the voltage at the 
other input being equal to zero. In such a case when sinusoidal oscillations hy = A;)sin wt at y = 0 are trans- 
mitted, it is possible to determine the relationship h/h,, = (jw). Similarly by feeding sinusoidal oscillation 
y = Bmsinwt to the other input with h, = 0, it ts possible to o the relationship h,/y = Ky (jw). 


If the system is linear (and experience shows that for small deviations it can be considered linear), the 
following equation must hold: 


h; = Kp, (jo) hy + K, (jo) y. (33) 


Under conditions of a full compensation of thickness variation h, = 0 


Ky, (ja) 


Expression (34) is the compensating device equation [see formula (20)] and relationship 


Ky, (jo) 


Ko (jo) = (35) 


represents the amplitude-phase characteristic of a compensating computer. 


Thus, with the knowledge of the frequency characteristic of the system for each of its inputs, it becomes 
possible to deduce from (35) the required frequency characteristic of the compensating computer, This charac- 
teristic can also be found by a somewhat different method by feeding oscillation = A,)sinwt and y = Bp sin 
— g) simultaneously to the inputs of the system in Fig, 8 and selecting values for B,, and ¢ for each value of 
so that h, = 0, 
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Fig. 14. 


Figures 9a and 9b show the amplitude-phase characteristics 
Khy (jw) and Ky(jw) obtained by means of an analog and Fig, 3.0 
shows their ratio, i. e., Kg(jw), the amplitude-phase characteristic 


of the compensating computer. The frequencies on these curves are 
in cps. 


By choosing suitable parameters in formula (21) it is possible 
to obtain an amplitude-phase characteristic very close to the one 
shown in Fig, 10 and thus obtain a synthesis of the compensating com- 
puter, This task can be simplified if the amplitude-phase charac- 
teristics shown in Fig. 10 are first analyzed, Curve Khy (Jw) can be 
approximated with sufficient accuracy to a curve of the amplitude- 
phase characteristic of a series connected first-order lagging element 
and a constant delay element with a time delay of r » 0,568 sec. 
This value is obtained in the electronic analog of the rolling mill 
by means of delay units CDU-1 and CDU-2, The lagging element 
paramters are found by determining the limiting (i. e., having a 
zero delay) characteristic Kim) of the curve Khy (jw). The time 
constant of this element also depends to a small extent on minor 
capacities introduced into the feedback circuits of some of the analog 
amplifiers in order to prevent oscillations, 


Curve Ky (jw) can be sufficiently well approximated by a . 
second-order element with a transfer function 


k 
Ky (P)~ EO OM (96) 


Now let us determine the parameters of the compensating 
computer, For this purpose, let us represent the system of the rolling 
mill electronic analog and the compensating computer by the block 
schematic of Fig. 11, taking into consideration the assumptions made 
with respect to operators K hy (P) and Ky (P). Obviously, the com- 
pensating computer circuit must contain an element with a transfer 
function 1/Ky(p) and units simulating the element with delay r and 
the lagging element with a transfer function k, /(1 + pT), which 
forms part of Kh, (p). 


The element inside the compensating computer surrounded by 
a dotted line was designed as follows, Let us examine the circuit 7 
shown in Fig. 12, If voltage ¢ is taken to be equal to zeio, it is 
easy to see that the following equality will hold: 


U2K(p) __ U 

(37) 
tp) 
(p) = = — (38) 


Hence, knowing the transfer function Ky (p) [see (36)] and the circuit which produces it, it is possible to 
obtain a circuit which will reproduce function 1/ Ky(p). Figure 14 shows a well-known circuit which simulates 
asecond order element of the form Ky(p) = + CaRyp + 
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By connecting this element into the feedback circuit of Fig. 12 instead of the element of the transfer func- 
tion K(p) and connecting an inverter with resistances p, and pz, required for obtaining a negative feedback, we 
obtain the required circuit (Fig. 14). A lagging element with a transfer function k,/(1 + pT) (Fig. 11) is connec- 


ted at the output of this circuit. The complete circuit corresponds to the part of the compensating computer sur- 
rounded by a dotted line in Fig. 11. 


Since the tests were made on an electronic analog of a rolling mill with a constant strip speed, a unit with 
a constant time lag and without an adjusting device, based on the known approximation of the function e~M(T~"%) 
by a ratio of two polynomials of p was used for an element with time lag r. 


It should be noted that the synthesis of computers for compensating input strip thickness variations can be 
carried out by means of more complex and universally applicable devices with respect to obtaining the required 
amplitude-phase characteristic [14]. Other circuit combinations are also possible, including more economical 
ones with the use of discrete and pulse techniques. 


Tests of the compensating computer circuit, carried out by means of the rolling mill electronic analog pro- 
duced the following results; this computer compensated variations in the strip thickness h on an analog ofa three- 
stand rolling mill in the frequency range of 0-2 cps with a limiting amplitude of about 25 yu [when this limit was 
exceeded some of the quantities attained their limiting values, see (29)], The thickness variations at the output 
were reduced to a fraction of their input values and were practically eliminated (Fig. 15), 


Each of the oscillograms of Fig. 15, taken for several values of f , the repetition frequency of the strip thick 
ness variations, shows on its left-hand side, values of hg, hy and y compensated by means of the computer. A 
good compensation is observable at the output of the rolling mill (hg » 0). 


Each oscillogram shows on its right-hand side, values of hg without compensation, i, e., when y = 0. In 
this case, h, differs considerably from zero. Thus the oscillograms of Fig. 15 show the possibility of compensa- 
ting for the strip thickness variations by means of a compensating computer, Small changes (up to 5%) in the 
stand and drive parameters do not affect substantially the conditions of compensation, 
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A COMPUTER FOR CONTROLLING THE POWER SUPPLIED TO 
A STEEL-SMELTING ARC-FURNACE 


Yu. M. Alyshev, L. N. Fitsmer, and L. I, Shevchenko 
(Moscow) 


Design principles of a computer for controlling the power supplied to a 
steel-smelting arc-furnace are cohsidered. Results of testing the device with 
a commercial furnace are given. 


Block Schematic of a Computer Automatic Control System 


The electrical operation of existing arc furnaces is controlled by sustaining set values of currents, or more 
precisely, relations between the phase currents and voltages [1]. 


Tests have shown that such a method of controlling the power supplied can vary considerably from smelt 
smelt, This is due mainly to supply voltage variations and the imperfections of existing regulators, 


Owing to the presence of many disturbances, the operators are unable to evaluate the actual value of the 
current and power and thus correctly set the regulator. 


t In this connection, Yu. E. Efroimovich posed the problem of control- 
“ co : ling the electrical operation of the furnace by sustaining the power supplied 

5) to it [2]. This problem was solved by introducing into the existing contol 
Y system, consisting of the controlled object CO, the regulator R and the as- 
signment unit K (Fig, 1), a computer (encircled in Fig. 1 by a dotted line), 
The latter consists of a multiplying unit MU, integrating unit IU, and a 
setting unit SU, 


The purpose of the MU element is to multiply the instantaneous value 
of current I and voltage U. The output MU voltage, proportional to the 
arc power P, is compared with the voltage proportional to the set power Po. 


The difference of these voltages: 


AP = P,—P (1) 


is fed to the input of the integrating unit IU. 


In this case, voltage a, at the output of IU is proportional to the deviation from the energy supplied to the 
furnace: 


=a APdt. (2) 


The set voltage to current ratio 
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U 
k= > (3) 


can be changed by the action of the two voltages a,» and a. Voltage a» corresponds to the set value of k and ts 
supplied directly from the setting unit SU; voltage a, corrects k when for any reason the energy actually supplied 
differs from its set value, Variations of coefficient k under the effect of voltages a, and ay is attained as it ts 
shown below by means of a saturation coil. 


The above block-schematic provides a constant set power supply in one of the furnace supply phases, When 
controlling the total power of a multi-phase furnace, the voltages proportional to the separate phase powers should 
be added; the voltage sum thus obtained should be compared with the voltage proportional to the set total power. 


The integration of the difference of these voltages can, obviously, be carried out by one integrating unit 
common to all the phases. 


Computer Schematic Circuit 


The schematic circuit of an arc furnace computer is shown in Fig, 2. The circuit calculates at every in- 
stant the total power supplied to the furnace, compares it with the set power and integrates the difference between 
them. The signal fed from the computer output operates the furnace electromechanical regulator (EMR) so as to 
eliminate the difference between the supplied and the set power. The power dissipated in each phase of the fur- 
nace is calculated by means of the multiplying unit MU [3]. For this purpose voltages proportional to the current 
and voltage of the furnace phase are fed to its input. These values are obtained by means of the instrument volt- 
age transformer VT and the current transformer CT. The products of the voltages of both polarities for the nor- 
mal operation of the unit, are obtained by grounding the midpoint of the voltage transformerVT and the potential 
divider R, and R; on the secondary winding of the second current transformerCTs. 


The multiplying unit obtains its set reference voltage U, =—190¥v, proportional to the set phase power of the 
furnace, by means of resistor Ry. As the result of this, the MU output provides a voltage proportional to the devia- 
tion of the actual phase pow=r from the set value. In the steel smelting process ,it is necessary to change the sup- 
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plied power within wide limits, If these power settings are made by means of varying Uy, or Ry, the MU will ob- 
tain,at certain stages of smelting,small input voltages as compared with the scale reading and will therefore pro- 
duce considerable errors, For this reason the reference voltage in this computer circuit is kept constant and the 
power is set by means of rheostat R, in the secondary winding of the current transformer CT . 


For this method of power setting the following equation holds: 
k,Um —U,=0. (4) 


Where Uy, is the voltage proportional to the measured power,and k, the coefficient depending on the value 
of resistor R. 


Circuits for measuring power in the remaining phases are designed in the same way, Computer circuits 
similar for all the three phases are encircled by a dotted line in Fig. 2. 


The output voltage of each phase MU is fed to the input of the operational smplifier Y. The voltage at the 
output of this amplifier is equal to: 


ou 


t 


where kg is the division ration, T the time constant of the integrating unit, P the real value of the total power 
supplied to the furnace, and P, the set value of the total power. 


Voltage Upyr, proportional to the supplied power deviation from its set value, is fed to the control winding 
CW of the saturated coil SC, The latter acts as a shunt in the current supply circuit of the electromechanical re- 
gulator, The shunting loss of the saturated coil is a function of the currents through windings CW and SW. The 
nominal furnace phase current is set by means of rheostat Rs. For each value of the power supplied to the furnace 
there is a definite relationship between the position of the sliders of rheostats Ry and Rs. When disturbances arise 
in the furnace, the signal fed from the output of amplifier Y reacts on the sa- 
turated coils in such a way that the electromechanical regulators compensate 
the deviation of the power supplied to the furnace from the set value, In addi- 
tion to the units enumerated, the computer contains some auxiliary elements 
required for adjusting the circuit and monitoring its operation, Rheostat R, serves 
to adjust the time constant of integration, which is required for adapting the 
computer to various smelting periods, The controls of rheostats Ry, Rs, and Rg 
and those of switches S4-S6 are ganged. 


The computer supplies are obtained from electronically stabilized recti- 
fiers, The multiplying unit, which forms part of the computer, includes thyrites 
and one dc operational amplifier which is widely applied in the electronic simu- 
lation technique [4]. For temperature compensation of thyrites, thermistors are 
used, 


The integrating unit IU is also designed on the basis of a dc operational 
amplifier, Only the last stage of this amplifier is not completely conventional 
and is designed for higher than normal power with a more economical utiliza- 
tion of the supply energy [5]. In addition to the elements already enumerated, 
the schematic shows switches S1/1 to S1/7 which switch the computer unit cir- 
cuits so as to be able to set their output voltages to zero when the instrument is 
connected for use. Neon lamps LS;1 LS3, and LS, with resistors Rg and Rg provide 
a warning when the linear portion of the output voltage scale is exceeded. 


Test Results 


The computer described above (Fig. 3) was tested on a 20 ton arc furnace 


by a regulator type RMD-45. The accuracy of sustaining a given power , supplieé 


of the “Elektrostal'” plant, The electrical operation of this furnace was controlled 


to the furnace, by the existing method and by means of the computer was measured. 
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The mean value of the total power supplied at various conditions of operation of the furnace was measured, 
for this purpose, the increment of energy AW was measured over a determined time interval At. The mean power 
was determined ‘rom this data: 


AW 
P= (6) 


The table contains data characterizing the accuracy with which a given total power can be maintained, 
from one smelting to another, by means of the computer and without it for two basic operating conditions of 
smelting. 


Method of | Measured 
Smelting stage 
process values ing 
control 
Fusion 
Without CC P, kw 5076 4915 5440 5833 5843 5513 6130 4933 
At, hrs 0.8 10 14 08 #40 4.2 0.7 0.7 
w cc P, kw 5358 5382 5278 5354 5295 5278 5355 5323 
oF At, hrs 0.8 1.0 09 0.7 08 O05 O.5 0.5 
Refining 
Without CC Pi kw 1456 1681 1536 1635 1566 1528 1489 41467 
At. hrs 12 12 08 43 410 410 0.7 O68 
With CC P, kw 1606 1600 1600 1606 1600 1602 1601 1602 
At, hrs 03 04 06 05 05 06 O04 0.5 


It follows from the data obtained during fusion (the most difficult stage in smelting with respect to distur- 
bances in the furnace) that the maximum error in maintaining a given total power without the computer is equal 
to +11%, and with the computer it does not exceed + 1-2%, 

SUMMARY 


Exhaustive tests of an electronic computer carried out in the laboratory and with a commercial arc furnace 
showed that with the help of the computer it is possible to control with considerably greater accuracy the electrical 
operation of the furnace as compared with the existing method, 


The computer developed by us could also be employed with ferroalloy furnaces, 


Prolonged testing under factory conditions of the computer electronic unit models (the multiplying, inte- 
grating and summating units) point to the possibility of their application in equipment for production automation. 


The development was carried out with the participation of Yu. E. Efroimovich, A, A, Fel‘dbaum, K, S, Gensht, 
V.B, Tikhmenev, and Ya. P. Stiop. 
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ELECTROPNEUMATIC TRANSDUCERS OF THE IAT AN SSSR* 


Yu. V. Krementulo 


(Kiev) 


The paper describes the layout,principles of operation and experimental 
investigation results of an electropneumatic transducer developed by the Auto- 
mation and Remote Control Institute of the Academy of Sciences USSR, The 
paper provides a technique of determining the error due to the nonlinear charac- 
teristic of the primary pneumatic relay of the “nozzle-vane” type, and discusses 
the gain of the secondary amplifying pneumatic relay. The technique of design- 
ing a temperature compensating circuit is provided. 


In production automation, it is often expedient to use an electrical measuring system, but pneumatic equip- 
ment for regulating production processes, In such instances, it is necessary to have devices interconnecting instru- 
ments working on different forms of energy. Let us call a device which hes an electrical parameter as its input 
and compressed air pressure at its output an electropneumatic transducer. 


The application of electropneumatic transducers provides monitoring and control systems working at greater 
distances, makes the substitution of costly pneumatic lines by electrical ones possible, etc. 


Originally, the development of the electropneumatic transducers was carried out by the IAT AN SSSR, in 
connection with the electrical instruments used in the unified assembly system (UAS) of automatic control and 
recording. In the course of development, other possible spheres of application of the electropneumatic transducer 


became apparent. 
Two types of transducers, the fpp-1 and fpp-2, are described below. 


1, Electropneumatic Transducer Type fpp-1 


Layout and principles of operation. The schematic layout of transducer EPP-1 is given in Fig. 1. 


The transducer slide wire Rs, is connected to the electronic instrument setting rheostat Rs,.* The position 
of the slide wire Rs, slide 1 is proportional to the measured parameter, When slides 1 and 2 of slide wires Rg, and 
Rg, are out of step an unbalance voltage is produced. It is amplified and fed to the winding of a reversible motor 
DA-1, which rotates by means of reduction gear P slide 2 in the direction which reduces the unbalance. Thus, the 
position of slide 2 is proportional to the parameter being measured, An Archimedes spiral-shaped cam K is 
mounted on the same spindle as slide 2. The rotation of the cam moves vane V of the primary pneumatic relay 
changing the opening of nozzle N. The interthrottle chamber IC is connected to the secondary pneumatic relay 
SR, . When the vane is approached to the nozzle, the pressure in the interthrottle chamber increases, This depres 
rod R, which is fixed on two membranes M, and Mg, and opens valve Vv. The air from the feeding line Px passes 


* The UAS electronic instruments type P-1, P-2, and P-2V have three slide wires, whose slides are mechanically 
ganged: the first slide wire forms part of the measuring circuit, the second provides the unbalance signal and the 
third is used as a setting rheostat, Setting rheostats are provided in other mass produced instruments, for instance, 
in EPD-227, EPP-120 and others, 


* Automation and Remote Control Institute, Academy of Sciences, USSR 
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Fig. 1. Layout of an electropneumatic transducer EPP-1, Rg, ts the setting rheostat, 
Rs, the slide wire of the electropneumatic transducer, A the electronic amplifier, PP 
the power pack, G the reduction gear, DA-1 the two-phase reversible motor, W,. and 
Wy the supply and control windings of DA-1, K the cam, V the vane, FU the feedback 
unit, FT and N the fixed throttle and nozzle of, respectively, the primary pneumatic 
relay, IC the interthrottle chamber, SRythesecondary pneumatic relay, Poy the out- 
put pressure, and P¢ the feeding pressure. 


> ks Jem? between the valve and its seating to the transducer output and on 
out’ to the feedback unit FU. This pulls the nozzle away from the vane, 
since the former is fixed to the moving bottom of the feedback bel- 
w lows, The transducer has a displacement negative feedback, 
WA Thus, at the output of the transducer a definite pressure is 
uw y, established, which depends on the position of vane V and the stiff- 
Vee ness of the feedback bellows springs, During operation, the nozzle 
P, "follows" the vane with an accuracy of 0.015 mm, I, e., the dis- 
ad placement of the bellows can be considered to be the same as that 
Yo of the vane (the working stroke of the bellows is 0.6 mm).* The 
@ LZ ' displacement of the vane and hence, the transducer output pressure 
“a is thus proportional to the measured quantity, 
Experimental investigation results, The basic transducer sta- 


L7/A tistical characteristic is given in Fig. 2, Ratio Ry/R,, where R, is the 


Fig. 2 ’ resistance of arm I, is plotted along the x axis, and the output pres- 


tion from linearity does not exceed 0.5%, This characteristic was obtained as the result of many measurements, 
The transducer error does not exceed 0.7%, the dead zone lies within the limits of 40.25%, The dead zone was 


determined as the minimum relative change in resistance of Rp, which produced a noticeable change in the output 


pressure of the transducer, 
A change in the supply pressure from 1.1 to 1.6 kg/cm? produced a mean error not exceeding 0.5%, 


*The displacement of the feedback bellows when the pressure in the bellows chamber changes from 0.2 to 1 kg/cm? 


ls considered to be the working stroke of the bellows. 


sure along the y axis, The characteristic is almost linear, the devia- 
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2. Electropneumatic Transducer EPP-2 


Layout and principles of operation, The transducer schematic layout is shown in Fig. 3. 
It receives de signals, The current flows through winding ot he the control coil which is suspended on two 


flat springs S in the airgap of a permanent magnet M. Vane V of 


+ 
‘y 


SRe 


fed 


a5 a 5 790 
f cps 
- 
Fig. 5. EPP-2 frequency characteristic, A/ A, is the 
ratio of the P,,,,, variations in amplitude to its value 
at zero frequency, ¢ is the phase difference between 


Pout and ty. 


e primary pneumatic relay is fixed directly 


2 2 
Pout” kg/em' 
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las 
122 
Fig. 4, EPp-2 characteristic. 


to the coil, To each value of current {_ in the control 
coil winding there corresponds a definite position of the 
vane, and hence, a definite output pressure Poy. In 
transducer Epp-2, similarly to transducer £pp-1, a dis- 
placement negative feedback is used, The pneumatic 
part of the layout is similar to that of Epp-1, The 
transducer is adjusted by means of the screw PB. 


Experimental investigation results, Figure 4 shows 
the statistical characteristic of the transducer, The 
transducer ‘is adjusted so that at ly = 0, the output chara 
teristic Pout = 0.6 kg/cm*, The output characteristic 
Pout = f (ly) is practically linear, the deviations from 
linearity do not exceed +0.7%, The current range ts 
+3 ma, 


Transducer tests have shown that temperature change 
of +2 to +50°C do not produce a noticeable error if the 
coil current remains constant, If, however, the trafs- 
ducer operates on voltage and not current, a tempeé- 
ture compensation circuit must be used, Its schematle 
is given in appendix 3. 


Figure 5 shows the amplitude- frequency and phase 
frequency characteristics of the transducer for the 7 
dition when output volume capacity was about 600 cm. 
The EPP-2 technical specification is given in the table, 


8. Application of Electropneumatic Transducers 


Some of the possible applications of electropneumatic transducers are shown in Figs, 6, 7, and 8, 
Figures 6 and 7 show the application of transducer type EpP-1, 
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EPP-1 


EPP-2 


Setting rheostat resistance measurements, Resistance 
100-300 ohms, supply voltage 12 v, 50 cps obtained 
from the transducer power pack. Use of inductive po- 
tential dividers possible. 


Direct current + 3 ma 


Air pressure (0.1-1)kg/cm*® 


Air pressure (0,64 0.4) 


kg/em* 


Basic error 


Additional error due to 


a 10% change in pres- 
sure 


Working temperature 


+10 to +50°C 


Ambient medium 


Neutral 


Supplies Air pressure 1.3 kg/cm*; supply voltage 220 v, 50 cps Air pressure 1.3 kg/cm* 

Weight = 5 kg m2 kg 
Remarks Slide wire, motor, electronic amplifier, power pack the | Power dissipated by the con- 


same as for electrical instruments of a UES, Used in 
the transducer as complete units 


Fig. 6. A control circuit using an EPP-1 transducer, 1) Controlled object, 2) elec- 
tronic measuring instrument, 3) electropneumatic transducer, 4) pneumatic control 
unit, 5) membrane actuating mechanism, Rs,) the setting rheostat, a) "measuring," 
b) “setting,” T) transducer, 


The circuit shown in Fig. 6 refers to the control of any process by means of an electronic measuring instru- 
ment, an electropneumatic transducer and a pneumatic actuator, Air pressure P; at the output of the transducer, 
Proportional to the value of the controlled quantity, is fed to the measuring chamber of the pneumatic control 
unit, Pressure P,, proportional to the set value of the controlled quantity, is fed to the setting chamber of the 
unit.‘ The output of the pneumatic control unit is connected to the actuating mechanism, 


Figure 7 refers to the case when the control law is worked out by the electron regulator, At the output of 
the electronic measuring instrument connected to the transducer, we obtain a voltage proportional to the difference 
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Fig. 7. Control circuit using an EPP-1 transducer (notation is the same as in Fig. 6). 3A is 
the electronic control device, P the pneumatic part of the £pP-1, 


between set and measured values of the controlled quantity, This voltage is fed to the electronic proportional 
integral regulating device, Its operation time is set by varying the ratio of the reduction gear G and the feedback 
coefficient is set by rheostat R,. 


The transducer output voltage is fed to the membrane actuating mechanism (MAM). 


Fig. 8. Control circuit using an EPp-2 transducer (notation is the same 
as in Fig. 6). 


A circuit using transducer £PP-2 is given in Fig. 8. A signal proportional to the deviation of the regulated 
quantity from its set value is fed from the setting rheostat to the input of the transducer, At the output of the 
transducer we obtain signal* 

—k(U —U,)+ Po, 


where Us is the voltage proportional to the set value of the controlled quantity, U, is the voltage proportional to 
the actual value of the controlled quantity, k is the voltage ratio, and Py, is the transducer output power at zero 
current in the control coil, 


The summation unit of the pneumatic regulator forms signal: 
A = ky (P — Ps) = [P, — Po + k(U, —U, 


* If the input resistance of the transducer is large compared with that of the setting rheostat. 
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where k, is the power ratio. 


In order to make A proportional to U, — Uj either the setting chamber should receive a pressure equal to 
Py {then A = ky {Po — Po + k(Us— Uj)} = kg(Us— Uj], or the reference point in the pneumatic regulator should be 
displaced by the amount Ps— Po, 


In conclusion, let us point out some other instances where electropneumatic transducers may be applied, 


Often electronic regulator circuits are designed in such a way that the feedback rheostat of the electrical 
actuating mechanism forms part of the electronic regulator circuit, If a pneumatic drive is required, it is possible 
to replace the electric drive by an electropneumatic transducer type Epp-1¢ connecting its output to a membrane 
actuating mechanism. 


Type EPP-2 transducer can also successfully be used for remote setting, remote control, and as a transition 
unit between pneumatic equipment and electronic analogs,** or as a summator.*** This type of transducer can 
also be used with automatic compensators which have no balancing motor or rheostat, but provide an output current 
proportional to a thermocouple emf, and with similar devices, 


Appendix 1 
Error Due to the Nonlinearity of the Primary Relay Characteristic 


Even if the feedback has a linear characteristic in the pneumatic layout examined in the text, its resulting 
characteristic Poy, = (Ax) (where Ax is the change in the position of the vane) will be nonlinear, This is caused 
by the nonlinear relationship between the interthrottle chamber pressure*and the gap between the vane and the 
nozzle, 


The variation of the gap between the vane and the nozzle is 
Adz = Ary— An, 


where Ax yis the displacement of the vane, Axy the displacement of the nozzle. 


Hence, 


Ax, == Az + AAg. 


Az, =kiPy AAr= +a, 


where kgP.,,¢ is the linearized characteristic of the primary relay, and a its deviation from linearity. 


If the maximum deviation from linearity is equal to a,,,, the error given above will become 


@m ax 


Here APour = Pout max ~ Pout min: 


Quantity k,AP,.., = Axay represents the working stroke of the bellows when the pressure changes by APoyt, 
and kKpAPo yt = AAXxp is the aan in the working gap due to the linearized characteristic, 


The transducer relative error due to the relay is 


ax 


If Axay > Adxp, the error becomes: 


* In this case the electronic amplifier of the EPP-1 transducer is no longer required, 

** Thus in one of the IAT AN SSSR laboratories, type Epp-2 transducer was used for investigating the operation 
of an experimental pneumatic regulator with an dectronic analog. 

*** In this case the coil must have several control windings, 


But 
| 
207 


a 
— 
Brgy 100%, 
Applied to the EPP-2 design 
APout = 0.8 kg/em*, ogy = 0.001 mm, Xgy = 0.5 mm, Ax, = 0.015 mm. 
For 5 Ax, 0.16%. 
Appendix 2 
Secondary Relay Gain 


Often the secondary amplifying relay is made with a pressure gain kp* equal to unity. 


In this case, adjusting difficulties arise due to the gap between the vane and the nozzle of the primary relay, 
Often it becomes advisable to make relays with a gain greater than one, 


Pours 


\ 


al 


a6 : Fig. 10. Temperature compensating circuit, 


In type EPP-1 and EPP-2 transducers the static 
characteristic of the relay is represented by the expres- 
sion: 


a 


az 


{ 2(Py —0,2) kg/cm’, 
aol Pout) 0 


4z,mm 
P.>0.2 
Fig. 9. Combined characteristic of the primary and > kg/cm : 
Py, <0.2 kg/cm 


the secondary relays. Ax is the gap between the vane 
and the nozzle of the primary relay, Poyt the output 
pressure of the secondary relay. Here the pressure gain is equal to two. 


Figure 9 shows the combined characteristics of the primary and secondary relays for kp = 1 and kp = 2, A 
increase in k, moves the characteristic to the right and makes a careful fitting of the vane to the nozzle unnece 
sary. In this respect the increase in kp has the same effect as using a primary relay with a constant pressure drop 
and either a constant or a variable throttle. 


Moreover, an increased kp requires a smaller gap and thus reduces the error discussed in appendix 1. 


Appendix 3 
A Temperature Compensating Circuit 


Let us now examine the temperature compensating circuit shown in Fig. 10. The circuit contains three re- 
sistances Ry, Rg, and Rs, which can be chosen arbitrarily, Resistors R, and R, are made of manganin and R; of cop 
per, Rj is the resistance of the control coil Wy- 


Let us assume the following requirements: 


a) with a given Uj, the current fy in the control coil must also be given; 


* Quantity kp is understood to represent the ratio between the interthrottle chamber pressure variations APy and 
output pressure variations APoyt. 


b) the input resistance of the circuit consisting of Ry, Rg; Rs and Ry, must have a maximum value; 
c) the temperature must be compensated. 


The design of a circuit satisfying these requirements can be calculated from formulas 


Vin 
(m should be > 4) kio = Vm —1, 


kio— 


Under these conditions the circuit input resistance will be 


m 2kiop—i1 
Ain ky + 1° 


SUMMARY 


i. 2 The two transducers, types EPP-1 and EPP-2, have the same accuracy, The absence in transducer type 
EPP-2 of such a sluggish element as a motor provides a better frequency characteristic, 


2, Operating experience has shown that transducers of this type can be successfully usec in automatic con- 
trol and tracking systems. 
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RELIABILITY CALCULATIONS FOR GAMMA-RAY RELAYS 


A. G. Vasil'ev and K, S. Klempner 
(Kharkov) 


The relation between reliability and source activity for gamma-ray relays 
is considered, Formulas are derived and optimum counting rates are calculated 
for sets of reliability indices and relay parameters. 


Transducers in automatic devices that use radioactive sources often have advantages over other types of 
transducer, Gamma relays that use Geiger counters are now common [1-3]. These relays contain a y -ray source, 
a detector (one or more counters) and a two-state electronic circuit, 


The relay changes its state when a moving absorber enters or leaves its line of sight. The state with no 
absorber we denote by 1; that with an absorber, by 2. 


The device will work reliably only if the source's activity A is suitably chosen, 


The total counting rate n, (counts/sec) given by all m counters when the moving absorber is absent must be 
higher than the threshold rate ny, at which the relay transfers from state 1 to state 2, 


The activity is defined by 


4nL*ny 
1 
3.7-107 abm 2 (1) 


where Mp is the counting rate for which states 1 and 2 are equally probable, a is the working length of a counter 
in cm, b is the diameter of a counter in cm, L is the distance from the source to the relay in cm (L » a), q ls 
the number of y -rays of the ith energy emitted per disintegration, « is the counter's efficiency for these y-rays 
(i. e., the ratio of the counting rate to the incident y-ray flux), dj is the path length in the fixed absorber of the 
jth sort, in cm, 4; is the absorption coefficient for these y -rays in this absorber in cm~*,-and 3,7. 10" is the num- 
ber of decays per second for A = 1 millicurie, 


We combine some of the constants to give 


A>+. 


Ng = BAg and n = BA, 
If the device is to work for a time t, we must have that 


Ae > Ao, 
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where d is the decay constant of the isotope, 


Formulas (1)-(4) have been put as inequalities because radioactive decays are random in time, and so ny 
fluctuates. The fluctuations in ny must not be such that ny =< npg [of (5)], else the relay would operate wrongly. 
Hence, the conditions ny > ny and A > Ag. 


The source should not be excessively active, though, because this would cause shielding difficulties, would 
reduce the counters’ working lives, and might reduce the counting efficiency. We shall also show that it becomes 
less probable that the relay will stay stably in state 2 as the source becomes more active, 


We have thus to find the minimum activity (Amin), the minimum rate in state 1 (Min), and the maximum 
permissible activity (Amax) that will ensure that state 2 is stable. 


We deal here with the criteria for the maximum and minimum activities. 


Every such relay contains an integrator, The pulses that reach the integrator during its integration time* 
r(which number on average Ty) show a Poisson distribution. 


The integrator's mean output is proportional to the mean counting rate, We may assume that the output 
also follows a Poisson law. This is not strictly true, though. The problem cannot be dealt with rigorously by clas~ 
sical probability theory; the theory of random processes is needed. 


The treatment that follows is adequate for many practical purposes. 


1. The Minimum Necessary Activity 


Suppose that a source, of activity A mC, causes a mean count rate n, in state 1, Any possible fluctuations 
in ny must be such that the result 


thy (5) 


isof adequately low probability. 


The formula that shows whether we shall have (5) complied with for a Poisson distribution is 


W (xn, < tM) = 1 (6) 


If we assign some small value to the probability of rng, we may use the integral probability tables for the 
Poisson distribution [6] to find My =Tnyby interpolation, Then 


M 


0.5000 where Nin is the minimum necessary count rate in state 1. Then (3) and (4) give us 


0 

0 | 0.1587 the minimum necessary activity as 
0 | 0.0228 
0 

0 

0 


0.0014 
00000317 Anin= M 
3x 10-7 in 


(8) 


But the Poisson function is tabulated in [6] only up to M = 20. It is also impossible to solve (6) for rn, in 
terms of simple functions, We can use the Gaussian distribution for M > 20, since at these M's, with z <5, the 
Gaussian distribution is very nearly the same as the Poisson one, To use the Gaussian form we put (5) as 


3) 


m, —2V (9) 
4) 


* We assume here that r = 2RC (see [4, 5] etc.). 
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where v7 is the standard deviation, and z its multiplicity for a given reliability. 
The probability of (9)'s being complied with is given as a function of z by 


W (z) = 0.5 — O(z) = 0.5 — 


where #(z) is the normalized Laplace function (Table IV of [6]) and 2@(z) is the probability integral, 


The table gives some values of W(z). The value of z that is used [as also the value of W(rny < rng) in 
(6)] is determined by the purpose of the relay (i. e., by how undesirable the consequences of false operations are), 


We transfer from inequalities to equalities in (9) and solve for ny; the minimum n, fs then given by 


tnt — + 2*) n, + = 0, (12) 


one solution of which, 


satisfies ny > Np. 
We get the minimum activity, by allowing for the count-rate fluctuations and using (3) and (4), as 


Amin= (Ao + je. (13) 


The second term in (12) is often very large (several times ng). 


For example, if r = 1 sec and the sensitivity is average (ny = 50 counts/sec), and the relay works in a con- 
" trol system in which frequent false signals are impermissible (z = 4), we have that n,,jn = 87.5 counts/sec, The 
value for ng = 15 counts/sec is, if the other conditions are unchanged, Nypin = 40.5 counts/sec. 


Figure 1 shows a set of N,,in = f (Mg) curves; these curves have been calculated for practically possible com- 
binations of thresholds, integration times and reliabilities. 


The ,,jn have been found from the Gaussian and Poisson distributions, The results are almost the same, 


The reliability, with ny, ny and r given, can be found from (6), or from the tabulated normalized Laplace 
function [6] in conjunction with the more convenient formula 


14 
(14) 


which can be derived from (9). 


Fig. 1, The minimum necessary count rates, 
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We checked that these relationships agreed with experiment, We used various values of ny, Mg, and r (rT = 
= 2RC) and recorded the state of a relay with an oscillograph, The experimental W(z) was measured as the ratio 
of the times spent in state 2 (wrong) and state 1 (right). 


The experimental W (z) agreed with those given by (14) and (10) used with the tabulated Laplace function, 


9. The Maximum Permissible Activity 


A source of the activity given by (13) will ensure that state 1 will be stable with a set probability, 


The stability in state 2 depends essentially on yd, which is a property of the variable absorber. We have 
assumed that the absorber is a homogeneous material, and that the y “rays are of one energy only, 


The relay will be stable in state 2 if the probability of 
tM, > Mo 


is sufficiently small; here ng is the mean count rate in state 2. 
if Mg = TN is small, ng cam be found from the Poisson law, here 


W (sng >> tM) = 
M, can be found by interpolation from the tables [6]. The maximum permissible count rate is 


ny = = 


The corresponding maximum activity is 


Amax= Ms oud, (18) 


We can calculate fiyjay by the method that we used for Nyjjn, if Mg is large. State 2 will be reliable if 
the chance that 


tM, + (19) 


will be complied with is small enough, 


Pax} 
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Fig. 2, Limiting permissible count rates, 
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We replace this inequality by an equality and solve for ng, with ng < np; the maximum permissible count 
rate is 


The maximum permissible rate in state 1 is 
d 
= 


The maximum permissible activity is found by putting (21) as 


2? —V 24 
Amax= (40 Je 


Figure 2 shows maximum permissible rates for z = 3 and for various values of yd andr. 


It is clear that relays with preset rT , Mo, wd and t can work reliably if the minimum necessary activity is 
less than or equal to the maximum permissible activity. In such cases the activity should be Ap jn. 


A given reliability is impossible within the limits set by these parameters if Apjn > Ayjay- The problem 
can be solved by using a less penetrating radiation, for instance. 


3. Corrections for Background 


The background caused by cosmic rays and natural radioactivity may affect the accuracy of these calcula- 
tions if the background rate m, is appreciable relative to n, or ng. 


The background depends mainly on the type of counter and on the number of counters used. The natural 
background {s stated on the counter's label; it is 0.1 to 1 count/sec. 


Then (9) and (19) become 


t(m +n 2V + + <M, 
+My) +2V E(t > M9. ptm 


If conditions (n, + my, > Mp and ng + Mm, < Ng) are suitable, the solution to these inequalities gives 


22+ V 24 + 4t (no + ny) 2? 
ny = M+ 2t — Np» 


224+V 24 + 4t(no + m,) 2° 
2t 


Ng = No + 


SUMMARY 
1, It is shown that the source activity must be chosen properly if a y-relay is to work reliably. 


2. Criteria for the minimum necessary and maximum permissible source activities are deduced and inves- 
tigated. 


3. A method is given for calculating the minimum counting rates that ensure reliable operation under mally 
of the conditions that are possible in practice (integration time and sensitivity threshold); values of the rates are 
calculated, 


4, The reliability is calculated for various values of the main parameters. Experiment confirms the cal~ 
culations, 


5. Limiting permissible count rates are calculated with respect to some parameters that determine how 
these relays work. 
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THE OPTIMAL PARAMETERS FOR NOISE STABILITY IN A 
TELEMETRY SYSTEM 


V. A. Kashtirin 
(Moscow) 


The stability of telemetry signals against weak and fairly strong fluctua- 
tion noise is calculated for pulse-time (PT) and pulse-width (PW) systems, It is 
shown that there are optimal transmission parameters, which ensure minimal 
referred root-mean-square error for a given noise level, The various types of 
modulation are compared from the point of view of noise stability, 


The various types of modulation have already been compared [1] for noise stability for the case where the 
signals are of limited dynamic range and the fluctuation noise in the channel is weak; the theory of potential 
noise stability was used, The referred root-mean-square error was anindex, This earlier treatment showed that 
the error caused by noise may be reduced by increasing the channel width, for most forms of modulation, The 
reason {s that the error is caused by the noise's altering the pulse front, if the noise is weak. The noise error ts 
increased, however, if the channel is made too wide, The cause is that false pulses are produced, or signal pulses 
‘are suppressed, Code- (C) and frequency-(F) modulation methods have been compared for stability against strong 
fluctuation noise by statistical and correlation methods [1, 2]. It has been shown that there are optimal transmis- 
sion parameters, which provide minimal noise error for both these forms of modulation, if theres a given noise 
level, It is of interest to derive the optimal parameters for PW and PT modulation systems, We have to derive 
the total error caused by the noise's distorting the pulse fronts, by false pulses, and by lost signal pulses. The 
optimal parameters are found from the minimum total error. 


1. Noise Stability of Pulse-Time Modulation 


Let us consider the noise stability of PT modulation; let us suppose that the noise can cause false pulses, 
and can suppress signal pulses, The effects of strong noise can be considered to be those of pulse noise with its 
pulses distributed randomly in time (mean number of noises per unit time S). The effects depend both on how the 
signal is received, and on how the synchronizing pulses are transmitted. 


We assume that we have a rigidly synchronized receiver which responds only to the first pulse that arrives, 
no matter whether it is a signal or a false pulse. Not every noise pulse causes an error at the output of such a re- 
ceiver, even if the pulse arrives during a modulation interval T. If the noise pulse comes after a signal pulse, the 
receiver will react only to the first. If the noise pulse comes before the signal one, the receiver reacts only to 
the noise pulse. The noise can cause the receiver either to respond to a false pulse, or to fail to respond to a 
signal pulse, 


Let us consider the error that occurs in the above method of detecting PT signals; the signals are taken © 
be time-modulated triangular video pulses of amplitude 2U;, and of length r. (We have taken the amplitude as 
2Um for convenience in comparing the various transmission methods, with the dynamic range of the signal, kept 


constant.) Suppose that some value of the measured parameter, say A, is transmitted, A noise pulse that occurs 
at an arbitrary time t;< t, causes an error 


A= ty foe ti, (1) 


where ty is the time when the signal pulse corresponding to \, is transmitted. We consider the case when the 
wobability of there being a false pulse during the time T is small, t.e., ST« 1. 


The probability of there being a false pulse in the time from t;— At; to t fs 
Py, = (2) 
We now assume that the probability that a false pulse will occur in time r {s equal to the probability that 


a pulse of leagth r will be lost, and that S (the number of false pulses per unit time) is defined as the probability 
V(a) that there will be one false pulse in time r when an ideal Kotel'nikov receiver [3] is used: 


(3) 


(4) 


Here o {s the spectfic noise voltage [in (cycles /sec)”*] and Q, is the signal pulse energy, which {s given 
by 


=2 (Sy dt = ba 


We use (4), (5) and the generalized transmission parameters p and y [1] to get 


V 


Fim is the highest frequency of the telemetered parameter, 


Af 
* 


and Af is the signal's bandwidth. 


We neglect the chances of two false 1 cee in time T, and of a signal pulse being lost at the same 
time as a false pulse occurs; we determine the error A’pr, which is the sum of the squares of the reduced root- 
mean-square error and of the reduced mean error, A%_,, and A,,, respectively, We average the error of (1) over 
all possible . and t,,and use (2) to get 


“S=+V@). = 
Here 
V'(a) = 2 dz 
is a tabulated function whose argument is, for our signal, given [4] by 
a? = ( = 
} 
Qe (5) 
where 
U 
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t 


—t ty T—t 
= S (th — t,)* dt; dt, V (a) ti dt,. (9) 
0 0 


Here the first term represents the error caused by the false pulses, and the second term that caused by lost 
signal pulses, Integration gives 


If 1/y « 1, (10) simplifies to 


ty T+ 


The reduced mean~square error is 


Now (11) and (13) show that when ST « 1 (1. e., when yV(a) « 1) we have At, py « We can 
therefore assume that 


Now (14) shows that the error increases with the time deviation, We have here assumed that the probability 
of two false pulses appearing in time T is much less than that of one, 1. e., that 


Cin V2 (a) (a) or (a) <1, (15) 


4 
Chr [y] (ly (16) 


and the [y] are integral values of y. 


It is of interest to see what is the minimum signal-to-noise ratio at the receiver input for which (15) is 
complied with. Suppose that we have 


(a) =0.1. (17) 


| 
The reduced mean error is | 
A i 
ms ~ 7a (12) 
0 
$0 
| as 
PT = ~ pr: 
where 


The input signal-to-noise ratio can be expressed 
in terms of the generalized transmission parameters: 


(18) 


We give y various values to find p and the mint- 
mum ratio for which (17) is complied with, Figure 1, 
curve A, show = f (y) for yV(a) 0.1, Then 

n 
(11) gives us the error caused by the fluctuation noise 
accurate to 5%, for signal-to-noise ratios defined by 
ww > this curve, 


Fig. 1 The reduced root-mean-square error caused by the 

tes pulse fronts’ being shifted by noise when the signal's 
peak power is limited (the signals were triangular video pulses of amplitude 2U;, and duration r) has been de- 
duced [1] as 


(19) 


Errors caused by noise pulses and by distortion to the signals’ leading edges may be considered as mutually 
exclusive events, Then the total error caused by leading-edge distortion, by false pulses, and by lost signal pulses, 
is found by averaging over these components, with due allowance for their statistical weights. 


The symbols we shall use are P¢(X,) = St, ( the probability of there being a false pulse in time t,), ty, 
Py (Ay) (the probability that a signal pulse will be lost), and 5,,5¢,5, (the quadratic errors caused by leading- 
edge distortion, by false pulses, and by lost signal pulses). 


If dy), is transmitted, the expression for the total error can be put as 


AL pr (hx) = [1 — Pr Pz Me Oe) + 
+ Py (de) + Py (x) (20) 


Now ST« 1 and, moreover, St, <« 1, 80 we can neglect Pr(A,) and P7(A,) relative to unity, We then 
have 


AS pr(dx) = +P fru) Be) + Pp Ou). (21) 


We must average over all transmitted values of ) in (21) in order to get thereduced root-mean-square error 
for all 4; we must also average over all possible positions of the false pulse in the time T: 


T—¢? 


(te — ti) dts + V (a) th dts. (22) 


Here the probability that a false pulse will occur in the interval ty < t< ty + (where t is any 


moment in the range 0 < t < t,), subject to the condition that only one false pulse occurs in time 0< t< t. We 
calculate the integrals in (22) and get 


V (a) 


As + (23) 


7 
10 a 
‘ 
| 
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ms PT = | | 
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If y is large, 1. e., y >» 4 (as is common in indus. 
trial telemetry systems) (23) simplifies to: 
/ 2 4 6810 28 @ 60000 20 


0.06 , V (a) 
A’ PT = (24) 


Figure 2 shows Ay py = f(y) for limited peak 
signal powers and for various values of p. The curves 
have been calculated from (23), It is clear from Fig, 2 
that, at the noise level assumed here, it is advantageous 
to widen the frequency band only up to some limit (Yopt) 
above which limit the error caused by the noise rises 
rapidly. There is thus some optimum time deviation 
for the PT system, and hence an optimal frequency band 
(y= Yopt)» i, e., one that gives minimal noise error, 
The optimum time deviation can be found as a func- 
tion of p from (23), 


Figure 3 shows how yopr varies with p for PT sys- 
tems; within the range in p covered,the curve can be 
approximated by 


oh? 
opt 25.2" (25) 


Figure 3 also shows y p+ V8. p for frequency modu- 
lation (FM) and code modulation (CM) [1, 2], for con- 
venience in comparing PT and PW systems with others, 
The symbols used for code modulation are CAM (code 


Fig. 3. 


modulation followed by amplitude modulation),CFM (code followed by frequency modulation), and CPM (code 
followed by phase modulation); cy, is the mark-to-space ratio of the binary code. 


/ 


The broken lines A and B in Fig, 2 show the regions where our formulas are applicable: 1) where the noise 
is weak [error mainly due to the first term in (23)], 2) where the noise is fairly strong [error mainly due to the se- 
cond term in (23)], and 3) where the noise is strong. In the last region errors caused by two false pulses appearing 
in time T should be allowed for. The noise error is then so large, though, that the system is unusable, The error 
can be reduced by increasing the signal voltage, or by using other transmission methods, 


Let us determine the noise stability of a PT system that has a limited mean signal power P,,, pr. 'm our case 


Pp 407. 
(28) 


We denote the square 'm PY p*: 


We combine (23) and (27) to give the reduced mean-square error as 


A? pr 


where 


* 


(29) 


Figure 2 (broken lines) shows A*» py as functions of y for various p*. It is clear from Fig. 2, and from (28), 
that there is an optimum time deviation (or passband) in this case also, but the optimum value is less critical, 
The physical cause is that pulse energy stays the same as y varies,when the mean power is limited, whereas the 
pulse energy falls as y increases,if the dynamic range is limited. 


2. Noise Stability of Pulse Width (PW) Modulation 


We consider the case where the noise level in the PW channel fs fairly high, We assume that the receiver 
consists of a memory device, and of a device that measures for how long during a time T the signal has exceeded 
some threshold value, The threshold level is U,, = */, Uc, where Ug is the signal's amplitude. 


If the probability P, that a noise pulse will occur in time T is small, we can neglect the ——* P, that 
two noise pulses will occur in this time, Then 


P,=M(T)==V@), (30) 


where M(T) = ST is the mathematical expectation of a false pulse in time T 


We assume that the leading edges of the PT and PW pulses last for equal times, and that the probability 
that a false pulse occurs in a time r equals that of losing a signal pulse in this time: 


Pe(r) = Py(r) = V(r). (31) 


The error caused by a false or lost pulse does not depend on the time when the error occurs, We take the 
mean length of the interfering pulse as r, whence we have 


(32) 
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where a is the reduced error caused by one false 
pulse, 


The reduced error averaged over the working range 
in ) is zero, Hence, we have 


\ 


A* ims PW = =. (33) 


It is easily shown that the absolute value of the re- 
duced error, averaged over the range in ), is much les 
than the reduced root-mean-square error in our case 
(P, « Py). The formula for the latter error, as caused 
by noise that shifts the leading edge, is [1] 


0.12 
PW = (34) 


when the peak power {s limited. 


Errors caused by noise pile-up are neglected in (34), 
Now the errors of (34) and (33) are statistically inde- 
pendent, so the total error from these two sources fs the 
simple sum of the two: 


0.12 , V(a) 


Figure 4 (full lines) shows A> pw vs. y for various values of p. There is an optimum width excursion (or 
Yopt)» i. e., one which gives minimum noise error, Figures 2 and 4 make it clear that the optimum y for PW is 
larger than is that for PT (if p is the same for both), Figure 3 shows yont pw =f(p), which, within the range in 
p covered, is represented adequately by 


pls 


Figure 1 shows the limits to the applicability of the formulas for weak noise in PT, PW and FM systems, in 


U 
and coordinates, 
Un A f 


Let us determine the noise stability of a PW system with limited mean power Pm pw. For a signal of am- 
plitude 2U,, we have 


Assuming 3y >» 1, we have 


We now use (38) and (35) to find the error when the mean power is limited 


V pw 
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(34), 


5) 


(38) 


(39) 


a" (40) 


Figure 4 (broken lines) shows A’; pw vs. y for various p*. The optimal width excursion when the mean 
power is limited fs somewhat larger than when the dynamic range is limited, Figures 3, 5,6 and 7 show Ay min 
vs. p and A's min V8. p° for PW, PT, FM and CM systems; they also show yopt vs. p and Yopt V8. p*. These 
curves enable us to choose the bandwidths that are optimal as regards noise stability for these forms of modulation, 
subject either to limited mean power or to limited dynamic range, 


3. Comparison of the Noise Stabilities of Various Methods 
of Transmitting Telemetry Signals 


We have shown above that PT and PW systems have optimum bandwidths (ones that result in least noise 
error), It has also been shown [1, 2) that frequency~ and code-modulation methods also have optimum bandwidths, 
We now compare the minimum errors (for y = yopt) these systems give, firstly for limited dynamic range, and 
secondly, for limited mean power. 


Figure 5 shows how the minimum errors depend 
onp; we include code modulation followed by various 
forms of secondary modulation, with a mark to space 
ratio a, for the code, The curves relate to equal dy- 
namic ranges in the signals, It is clear, frequency and 
code methods are the best when used with their opti- 
mum bandwidths, Each type of secondary modulation 
has some range in p for which the code system is 
more stable against noise than is the frequency-modu- 
lation one, The range in which the code system Is 
most stable against noise is that where the errors are 
small (A < 107° 10~4, When < A < to 
10~*, the frequency-modulation system gives less error 
than the code one does, ~ 


0 20 bop 


We have calculated for the code system on the as~ 
sumption that T< 1/2 F,,,, which true if the re- 
corder has a storage device. If the recorder has simply 
a low-frequency filter, it is usual to have T < 1/2 F,,, 


(Ts om The error then increases substantially. 
-9)Fm 


4; min The PT and PW systems are equally stable against noise, 
Fig. 5, but are less than are the frequency- and code-modula- 
tion ones for all valuesof p covered, Figure 3 readily 
gives the optimum bandwidths for the various forms of modulation when the notse levels and response speeds (p ) 
are known; the demand for bandwidth is thus readily evaluated, 


Consider an example of calculating the errors of PT and PW systems, Suppose we have U,, = 90 mv, o = 
=2 mv -cps-1 and T = 0.1 sec. We have to find the mean-square errors and the optimum bandwidths. We find 
the value of p from (7), using F,, = 1/2 T: 


p = 20. 
Formula (25) now gives us the optimum y: 


Yopt PT = 6.46, 
The bandwidth (8) now gives us. 


where j 
| q 
: | 
the 4 = : 
or 
6) it 2 q 
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Afpr = 64,6 cps. 


We then get the error from (23) as 


V (4.54) [4 + 6.46] 
13 = 0.57%. 


0.25 
Va x 646 
The signal-to-noise ratio at the input to a receiver with a 64.4 cps passband is 


= 8. 


We now use the same data for a PW system. Here (36) gives 


Yopt PW = 11,3. 


Af pw =113 cps. 


Then (35) gives us the root-mean-square error: 


__0.12 V (3.45) 


The signal-to-noise ratio at the input to a receiver with a 113 cps passband is 


(ni! 


n.ef Af= 113 


It is thus clear that the errors are equal if the peak powers are limited and the optimum bandwidths are 
used, The PT system demands only about half the bandwidth, though. 


Figure 6 gives the minimum errors for the various forms of modulation, subject to limited mean power, The 
results show that PT systems give less error than PW at ones at allp, and if the fluctuation noise is very weak 
(As < 10-* to 10~), PT systems are the best (except for code modulation followed by phase modulation). If 
A= 107 to 107%, then frequency modulation gives 

. the least error, Figure 7 gives the optimum bandwidth 
when the mean power is limited. A point here is that 
code systems can be made more stable at any noise 
level, if codes more complex than binary ones are used, 


We have shown above that the optimum bandwidtls 
(y = yopt) gives the best noise stability, The various 
forms of modulation have been compared on the assump 
tion that optimum bandwidths are used, It is easy 
provide the optimum bandwidth in code modulation, 
since there are no special demands on the apparatus, 
It is often difficult to ensure the optimum bandwidth 
in PT and PW systems, because the apparatus does not 
function precisely enough; values of y > yop, are used. 
It is therefore of interest to compare the vations forms 
of modulation when the bandwidths exceed the optimum 
values, and the noise levels are fairly high. 


0 20 JO 


Here the errors of PT and PW systems are deter- 
mined mainly by the second terms in (23) and (35). If 

the peak powers are limited we find that, when y > 

> Yopt* the PW system always gives less error, for a 

given p and y, than does the PT, Hence, if the band- 

widths are to be larger than the optimum ones, PW should 

be preferred to PT. For example, if p = 20, (25) and 

(36) give 


100 


We use (23) and (35) to derive the errors for y > 
> Yopt: Suppose that y = 15. Then Ay py = 4.76% . 
and A 5 pw = 0.66%. Figures 2 and 4, and the results 
of [2], imply that FM systems have optimum bandwidths 
very much larger than do PT and PW systems, Hence, 
p and y values that cause PT and PW systems to work 
in the h igh-noise range (y > yp,) will cause FM ones 
to work in the low-noise range a © Yopv) The FM 
system gives then very much smaller errors than do the PT and PW ones, If the optimum bandwidth is for some 
reason impractical in a PT, PW or FM system, all these systems become less suitable than a code modulation one, 
for which there is usually no difficulty in providing the optimum bandwidth, 


Fig. 7. 


Noise stabilities can sometimes be evaluated in terms of the probability that an error exceeding some set 
value (say, Ay = AT /T) will occur, If the noise is weak, the probability distribution is normal, and if it is strong, 
the distribution is uniform, so for PT systems we have the formula 


App = — @ tt — + 27 V 7 ). 


The first term gives the error caused by false pulses; 
the second, that caused by lost pulses; and the third, 
that caused by shift to the leading edges, Figure 8 shows 
the probability that the error will exceed A, as func- 

tions of y for several values of p, for the case where . 
/ the peak power is constant, The curves show that there 
is a well-developed minimum in all cases, 


The optimum y (optimum, that ts, as regards the 
probability of the error's exceeding A,) depends on p 
and A;. For PW systems we must know the probability 
P(a, t)dr that fluctuations of duration range r ,r + dr 
will exceed the threshold level in order to find the 
probability of an error, The finding of P(a, r)dr has 
been considered several times, e. g., in [5, 6}, These 
papers indicate that, if the correlation coefficient of 


the noise is R(r) = (where and 
if r; = 1/Af, the mean number of times the level 
a is exceeded by fluctuations of duration range r,T + 
+ dr in unit time is given by 


ul 
2%, 


N (a. t) dt = — (41) 
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Formula (41) is valid if a® < 20, 


If the noise is not too strong, so that the probability of two false pulses arising in a time T can be neglected, 
the mean number of times the level a is exceeded for a time r > 1, will equal the probability that there will 
in a time T be an error greater than A, = 1,/T: 


_ 
P,(JA|>A,) =T =0.3% 


Thus we can use probability tests to evaluate noise stabilities, in addition to mean-square and mean error 
tests; the former are especially convenient for use with systems whose working 1s disturbed by errors that exceed 
some set value, These tests are also useful as extra ones, for use in conjunction with the others. 


We murc also consider how well the optimum parameters can in fact be realized when we use these proba- 
bility tests. ‘fhe actual conditions of use may have a major influence on the choice of the system that gives the 
best noise stability. 


1 am indebted to G. A. Shastova for some valuable commentson this paper. 
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NOISE STABILITY IN PULSE-WIDTH AND PULSE-TIME 
TELEMETRY SUBJECT TO STRONG FLUCTUATION NOISE 


N. V. Pozin 
(Moscow) 


Noise stability in pulse-width and pulse-time telemetry subject to strong 
fluctuation noise is examined; the main methods of transmission are considered, 
The mean and mean-square errors are used as tests of the noise stability. Errors 

in the leading edges of pulses are neglected, but allowance is made for errors 
caused by fluctuation noise which splits up the pulses or the spaces between them. 
The signal is sampled for analysis, 


Formulas for evaluating the noise stability are deduced, The noise sta- 
bilities of certain specific telemetry systems are compared, 


The mean and mean-square telemetry errors of pulse-width (PW) and pulse-time (PT) systems are dealt 
with here, for the case where fluctuation noise is fairly strong.* Here we consider steep-fronted pulses; errors 
within the front have no appreciable effect on the telemetry. 


Other noise stability tests may be better if pulse noise, or strong fluctuation noise, is present, or when some 
types of moculation and demodulation are used. Such tests are, for example, the probability of a false pulse, 
the mean number of false pulses, the probability of errors that exceed some set value, etc.** Some of these tests 
are derived as intermediate steps in this paper. ; 


Particular attention is given to the mean and mean-square errors, because these are the most suitable for 
use with many telemetry systems, This is especially so if integrating terminal equipment is used which has 
its integrating time T much larger than the time for a telemetry cycle r (T » 1,). Such devices are often used, 
The region where these errors increase rapidly (the noise threshold) is often a useful rough test even when the 
averaging times are short, where, in general, the magnitude of the errors is not of major importance, 


We shall first deal with the low-frequency side of the receiver, and shall assume 1) that unipolar pulses with 
silent spaces are received, 2) that fluctuation noise is present in this section of the receiver, and 3) that this sec- 
tion contains a low-pass filter with an upper frequency limit f, a double-acting limiter, which selects some nar- 
tow band at the level Up (Fig. 1) in the incoming pulse, and finally a demodulator, 


In a PW system the demodulator is equivalent to an instrument that measures the mean output current from 
the limiter, For a PT system, we envisage a demodulator that is equivalent to a PT to PW converter, and that 


produces width-modulated pulses with their origins at the PT working pulses, and with their ends at strobing pulses 
produced by the receiver. 


* We envisage fairly low signal-to-noise ratios (e. g., Uc/Un > 3-5), which Iie at the noise stability threshold, 
Correlation methods of reception are outside the scope of this paper. 

** The mean-square deviation is clearly most applicable {f the errors are distributed nearly normally, The mean 
and mean-square tests are unsuitable if the incoming pulses are recorded and then processed, etc, 


ed, 
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1. Pulse-Width Modulation 


We divide the channel time rj, (Fig. 1) into parts At = 1/,f and call At the length of unit (elementary) 
time interval, Let 


om (1) 


by which we mean that n is the number of unit intervals in time r,). 
We call"unif(elementary) distortion a false space 


= oe niieh ensovaies — of length At occurring at any point in the output pulse 
uff seEegsesec sssssesaeeecs from the limiter, or a false pulse of the same length 

7 in any part of a true space, We assume that the proba 
Foyt bilities (P, and Ps, respectively) of these untt distor- 


tions are known (see section 3), The At defined above 
is the correlation interval of the noise output from the 
filter, so the unit distortions can be taken as indepen- 
dent [1, 2). 


The mean number of such distortions in a space interval (rT, — r) and in a pulse interval (r) are respectively 


Fig. 1, PW. Division of the time channel into n parts 
of duration At = 1/2 f. 


m, = (n—m)P,, = (2) 


(3) 


We define the mean absolute error Agy as the difference between the mean numbers of unit distortions to 
space and pulse (i, e., between the positive and negative errors), multiplied by At: 


Aay = At(m,; —m) —+(Ps + A). (4) 


Now (4) shows that A,,y varies linearly with the pulse time r. In telemetry,it is usually some combined 
value of the error referred to the time range Ar (i, €., to the total range of change in the parameter) that Is of 
interest. 


We define the reduced mean error 5 gy as the ratio of the absolute error, averaged over Ar, to the time 
range itself. If Ar 


5 ay = ea (5) 


We average A,, over T,}, and integrate the modulus of (4): 


where rp is found from the condition Agy = 0, 
Integrating and using (5), we have 


| 

where 

: 
| 
Te t 

4 
d 

Sev" (1) | 


If P, = P ,orif Py» Pp 


Say = 2° (8) 
Passing now to the mean~square error, we consider a time T that contains N cycles of length r,: 


(9) 


yhere 1 is the number of channels, The time T specifies the lag in the telemetry The mean-square error A* rms 
js the dispersion in the unit distortions in time T; it can be expressed in terms of the dispersion for a time At 
(P, and Pp, respectively for space and pulse), in accordance with the rule for adding dispersions, as: 


= N (n—m) (10) 


The reduced mean-square error is found by averaging over Ar for X* rms and forrsing the ratio to the square 
of the maximum change in the parameter in time T, 1. e., to nnt(2) : 


A* ims 
To average pn we integrate (11) with respect to ny, assuming that n, is continuous: 
n 
0 
We substitute the result into (11) and get 
P, + P, 
If > Pp 
P, 
(14) 


2, Pulse-Time Modulation 


The time r between the received and strobing pulses (Fig, 2) is termed “the clear space," The interval 
Tt.) ~ 7 is the "the uncleared space." To find the errors we need to know the probability P({) that erroneous recep~ 
tion will occur from a false pulse that lies iAt away from the transmitted one, We consider the probabilities for 
the above two spaces separately; they are those for deviations fy and ig, where 0 < 1; < ny and 0 > fy > —(n— ny), 


We find P(i,) as the probability that two independent events occur together; the events are that there is no 
distortion in the ny — 1, section, and that there is unit distortion at a distance 1, from the leading edge: 


P(i,) = Py (1 —P, 


(0<i1<m)- ™ 


* It is convenient to represent T as the time taken to set up a steady reading in therecetver , if the links with 
lags in the transmission system are aperiodic, N is a coefficient that depends on the method of recording the read~- 
ings, T=r, and N = 1 if the system is optimal as regards response speed, 
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P(ig) is the probability that three independent eveny 
will occur together; the events are that the true pulse 
is lost, that there {fs unit distortion exactly away from 


the leading edge, and that there no distortion in 
— Toy! * section n 


U P = R, P, (1 — 
Fig. 2, PT. Division of the time channel into n parts (16) 
of duration At = 1/2 f. [0 >i >—(n—n)]. 
The main source of errors is clearly distortion in the uncleared space, In the future we shall assume that 


P(t) = P(ty). 
The absolute mean error (averaged over a range Ar T,») is 


(17) 


If 1 and ny change continuously (and not by steps), the sums can be replaced by integrals, Then (15) gives 


n ™% 
\ (4 —B)™ i(1 —B )~didn,. (18) 


™=0 i=0 


The integrals are readily calculated, Formula (27), derived jo the appendix, can be used to find the mean 
error, but it is rather cumbersome. Expanding In(1— P) and(1~— P) inseries form (for small P) we get (28), which 
is moreconvenient, The reduced mean error is determined from (18) and (28) (see appendix) by using (5). 


The simplest expression for the mean error (but one that gives the error too large) is obtained by replacing 
the double integral in (18) by its limiting value when P—> 0, 1. e., by (29). Then (18), (29) and (5) give us 


nP, 
Say = (19) 


The absolute mean-square error eS: » averaged over Ar = T,y, is given by 


(20) 


This integral is evaluated exactly and approximately (for small P) in the appendix [ Formulas (30) and (31)} 
We use the analogy with the largest limiting value for (32) in conjunction with (20) and (11) to get a very simple 
expression for the reduced mean-square error (which gives values that are too high, though): 


2 nP, 


rms 


(21) 


3. Some Remarks on the Initial Assumptions in the Analysis 


We have neglected distortion to the leading edges. PT and PW telemetry systems can use pulses with lead- 
ing edges that rise Very slowly or very sharply (relative to the time deviation), The choice of slope for the lead 
ing edges is not always made from considerations of noise stability. 


=0 i=0 
Bins \ didn, = 
m=0 i=0 
NP 
(1 — P = 
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The reduced telemetry error can be put as 


{, e,, a8 the resultant of independent errors caused by the transmitter, 54, by noise in the channel, 53, by variable 
relation between the levels of signal and of limiter (which defines the error in measuring the duration that occurs 
within the leading edge), 53, and, finally, by instability in the recording or indicating device, 5,. 


Slowly rising leading edges can be used {f the level at the limiter is proportional to the signal levels, 1. e., 
if Uy = k-U, (variation in k, determines 53), The relations that are used to calculate the main parameters r.h, 
Ar, low-frequency passband in the video section f) are then; 


4 
= SENT Ard, 2 (22) 


Here F is the highest frequency with which the telemetered parameter changes, A, is the absolute (maximal 
or mean) instability in the modulator, Typ in is the minimum pulse duration in PW modulation (the pulse length 
in PT systems). 


(22) shows that, if the modulator is good (small A), only a very narrow bandwidth is needed, 


Steep leading edges are used if U, and U, are rather unstable, to reduce errors in duration measurements on 
demodulation, i. e., to reduce 53 = Ay/Ar (here As = ker¢, where rr is the length of the leading edge and k is 
a coefficient; k < 0,5 usually, because of instability in k,). 


The relations to be used in the initial calculations on a telemetry system are 


4 Bat 
Toh = Tch 4* f (23) 


Let us now consider the effects of noise, Errors caused by noise (53) can be represented as due to leading- 
edge distortion and to false pulses (or spaces), If narrow-band transmission is used, with slowly rising leading edges, 
the first type of error is the more probable, The error must then be calculated* by a method in which this distor- 
tion is allowed for, e. g., by using Kotel 'nikov's potential noise stability theory [3]. Wider channels (1. e., steeper 
leading edges) reduce this type of error. The noise level is thereby increased, however, and so the probability 
of false pulses (spaces) rises, It may necessary to allow for both sorts of error under some circumstances, If the 
leading edges rise sharply, however, distortion in them can usually be neglected, 


Many existing telemetry systems use fairly steep leading edges, of very short relative durations, The for- 
mulas of sections 1 and 2 are then best used, with 5, = Say and 5,=4r,* in accordance with which test is the 
more convenient. 


We have so far considered the low-frequency side of the receiver, We must know the noise distribution in 
this,if we are to find P, and Pp, The Laplace integral gives us them,{f the distribution is normal, This usually im- 
plies direct transmission without secondary modulation, though; and this is rare, Secondary modulation is more 


common, The low-frequency section is then preceded by a filter of passband W and by the primary demodulator, 
The relation 


WwW 
(24) 
enables us to pass from one circuit to the other in our calculations, It is not always easy to find P, and Pp by cal- 
culation, since the noise and signal distributions after the demodulator are not always known, An approximate 
method for P, and Pp is available,if secondary amplitude modulation ts used with linear detection, and the noise 
in the high-frequency section shows a normal distribution [1, 2]. This special case is a common one, 


We assume that we know the effective noise voltage Up, the passband W, the signal level U,, and the level 
in the limiter, If this latter is Us for a space, it is reckoned as U, — Up, when a pulse is transmitted. Then P, for 
the space can be taken from the tables for the integral Rayleigh distribution 
* The root-mean-square error plays the main part here, 1. €., 53 = 5 rms. 
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Now Pp can be found from the generalized integral Rayleigh distribution 


U, 
it’ U*+U? UU, 
Vex (— ) av: (26) 


Here |, is a zeroth-order Bessel function. 

The integral is evaluated from the tables (or graphically, see, for instance, [1, 2]), so Pp is found as a func- 
tion (Up — U,,)/Unfor sets of values of U./Un. 
4. Noise Stability Comparison for PT, PW and FM Telemetry Systems 


The simplified formulas for PW, PT, and FM systems are given in the table; the values for the last group 
have been derived in an analogous fashion [4)].° 


of Here M = TAf, where T is the integration time 
Say 5rms and Af is the frequency range [4]. 
FM Ps Ps It is of interest to compare the noise stabilities of 
2 2M these systems. 

ne Ps Ps The noise stability depends not only on the type 
2 2Nn of modulation, but also on the methods of transmission 
nP, nP, and reception (i. e., of limiting, recording, etc.) if the 

PT oF Ton noise level is high. We must therefore allow for the 


actual ways in which telemetry systems work, We need 
% to know the numerical values of the main generalized 
Fi ' parameters M, N, and n, at least. We also need to know 
the ratio of passbands, Wi/W. We clearly need to com- 
Ka pare the systems in relation to devices that have the 
fa) same main errors, peak signals, etc,, as well as the 
| same methods of secondary modulation, 


We now list some features of the systems we have 
chosen to compare. 


\ 1, Comparisons are made againstFM,. We assume 
4 that the indicator system is solely responsible for the 
lag. The frequency sweep has a range Af = f..ay~fmin” 
; 3 and the input filter a passband W. 
\ \ 2. We assume that the PW and PT systems have 
indicators that are free from lag, and that can follow 
a 2 4 6 a 0 the change in the parameter during a single cycle, Then 
& T= Tchy and Ny=1. Hence (1) and (9) give ny= 
% = TW,/N,. The input filter's passband is the same as 
Fig. 3, Relations of 5,,, (curves a) and rms (curves for FM, i. e., W1=W. 
PW, ond 3. We consider also PW and PT systems which 


eT en, pemmene. differ from the above only in that the readings are taken 
from instruments that have the same lag as the FM system. 


* The simplified expressions have been derived for Py >» Pp. 


** We assume that f,,jn and T are appropriately matched, i. e., that at f,jp the reception is stable and the 
Output reading free from flicker. 


~ 
0 6 8 ” 4 


Fig. 4, Relations of 5, (curves a) andé,,,, (curves b) for: 1) FM, 2) PW, 
and 3) PT to signal-to-noise ratio referred to the frequency band for FM, with 
identical methods of indication, 


The cycle repetition frequency must be clearly not less than f,,jp in the pulse-time system, if reception is 
to be stable. Now, Tchyy = 1/fmin» Nu = T/Tchyy, taken with (23) and (24), give us 


Wry = ny. 


Hence, in the third case Wy is Nj)/N times larger than in the other two cases, Now Un = ow, where 
o is the noise intensity, so, if the conditions are comparable (i. e., the o are equal), we have the U,/Uy, ratio 


in this case a factor Nit / larger than in the first and second cases, 
I 


Be: 
Figure 3 shows day and 6,.,, vs. signal-to-noise ratio for FM with M = 15, and for PW and PT with n, = 
= 50 and Nj = 1, Figure 4 shows the same curves for FM, but PW and PT have here Ny = 5, 


_ We have found Ps for secondary amplitude modulation from (25), with Uy = 7,U,. The broken lines in 
Fig. 4 for PT with Nyy = 5 are for Ug = */,U,; they show how the noise stability depends on the method of recep- 
tion, 


The above values of M, N and n relate to specially designed telemetry systems which satisfy roughly the 
same requirements, 


SUMMARY 


The’ following conclusions are drawn from analyses of pulse telemetry systems in which steep~fronted pulses 
are used with limited peak signals in the presence of strong fluctuation noise, 


1, PW-AM systems are much more stable against noise than are PT-AM ones, 


2. FM-AM systems are usually more stable against noise than are PW-AM ones, if the mean current output 
from the demodulator is recorded, Faster recording instruments, e. g.,thoserecording the instantaneous current 
(the information within one measurement cycle) can be used to give devices of optimal response speed with PW 
systems (N = 1), This means that PW-AM systems can be made to be about as stable against noise as FM-AM ones 
are, The results are derived by considering the region where the errors increase rapidly, i, e,, the noise-stability 
threshold. 


£% | 


3, The mean and mean-square errors are of about the same order when the fluctuation noise is strong. 


mean error becomes the better test as the integration time increases, 


The 


4, It becomes particularly clear how the noise stability depends on the method of reception if the noise 


is of pulse or fluctuation type. 


APPEN DIX 
1, The double integral in (18) is evaluated by integrating by parts: 
n 


™% 
0 


n? i—(i1— 


If —In(i— P)=P, (1— 


h~ip—pt 
The integral becomes largest when P—> 0: 


n® 
The double integral in (20) is evaluated by integrating by parts: 


™ 


n 
n= — Py didn = 


n? n? 2n 


If P« 


nt 


gp — pr + — 
The integral becomes largest when P—— 0: 


n 
n* 
—{ 
lim didn, = 
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THE KST-1 REMOTE CONTROL SYSTEM FOR DISTRIBUTED OBJECTS 


V. A. Il'in, K. P. Kurdyukov, and V, I, Stepanov 
(Moscow) 


A remote control system for use with distributed objects is described; it 
uses purely electronic switching units. Frequency selectors (resonant electric 
circuits) are used with transistors in a new short-period pulse-time telemetry 


system. 


Local automatic devices can be used to perform remote control operations on objects any distance apart, 
in the same way as such operations are done within a small factory or works. Combined systems are then needed 
to signal the states of the remote objects, for telemetry, for remote control, and for telephony. 


A system in which all these functions are combined is the best for controlling the processes from a central 
panel without the need for staffing the remote objects; it is a step towards fully automatic systems, The KST-1 
bia 1 system largely fulfills these requirements; it was de- 
veloped at the Institute of Automation and Remote 
2 Control, Academy of Sciences, USSR, for remote con- 
2g trol of distributed objects, 


Studies on sound design principles for such systems 
[1, 2] have been used to base the KST~1 system on fre- 
quency-combination methods of selecting the objects; 
oscillatory circuits are used as the frequency selectors, 
The lines used to connect the objects may be run in 
various ways, Radio links can also be used; only slight 


changes are required. 
Fig. 1. Frequency selector in the KST-1 system, No local supplies are needed in the KST~1 system, 
PB) Local relay; TP) matching transformer, K, and This makes the system more reliable, extends its uses, 


Kg) series-resonant circuits; T) telemetry transmitter. and makes it available to objects that are not electri- 
fled (gas lines, gas works, irrigation systems, etc.), 


The reliability of the system as a whole is determined mainly by that of the units at the object sites, since 


these are unattended, These units have therefore been designed round highly reliable magnetic and semiconduc- 
tor devices in the KST-1 system. 


A frequency selector in the KST~1 has two resonant LC circuits which are tuned to the fixed frequencies 
assigned to their particular object. The circuits are coupled to the line via a transformer, Transistors receive 
the output voltages from the circuits, These transistors are connected in series with a relay that switches the ef- 
fector circuits (Fig. 1). Magnetic amplifiers can replace this relay. 


The object is activated by pressing its button on the control panel. The two fixed frequencies assigned to 
that object are generated and fed to the line. Frequencies from 70 to 300 cps are used in order to ensure maxi- 
mum range, The resonant circuits have passbands of 5-7 cps. Audiofrequencies can also be used if necessary 
without changing the design. 
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Fig. 2, Block diagram of the KST~-1 system, 1) Signaling unit; 2) auto- 
matic test unit; 3) automatic control unit; 4) control unit; 5) unit con- 
taining generators; 6) amplifier unit; 17) telemetry reception unit; 8) mea- 
suring unit; 9) comparator; 10) selector unit; 11) telemetry transmitter; 
12) operation selector; 13) object position transducer; 14) object number 
indicator; 15) effector; 16) measured parameter transducer, 


The signals cause the two transistors in series to conduct; these transistors then allow dc current to flow to 
the relay. The dc operates the relay, This dc method is a typical feature of these frequency selectors; it greatly 
increases the range that can be covered with a low-power generator, 


Telemetry, remote control (switching) or adjustment (smooth variation of the object's final position can also 
be established. Control and adjustment operations are effected by frequencies common to all the objects. The 
frequency here determines the type of operation. The signals act on frequency selectors which each consist of one 
resonant circuit only; these circuits cause the appropriate effector relays to operate. 


The VST-1 short period pulse-time telemetry system [3], which was developed at the Institute of Automation 
and Remote Control, Academy of Sciences, USSR, is used to measure and signal the object's parameters, The 
transmitter at the object consists only of highly reliable components, and contains no moving parts, contactors, 
or tubes, The extreme positions of the object are signaled on demand by connecting to the transmitter one or 
other of two fixed resistors, equal in value to the maximum and minimum resistances of the potentiometer trans- 
ducer, 


The telemetry system can also be used to confirm that the object has responded (checking), When the object 
responds, it connects a fixed resistance to the telemetry transmitter; every object has its own different resistance, 


The object's number is shown on the measuring instrument scale, the needle of which instrument shows the 
resistance indicated, 


Figure 2 shows the block diagram of the system. 


Extra standard units (shown by broken line) can be fitted to do routine telemetry, adjustment, etc, auto~ 
matically, The other units are also built up on separate chassis, whose number may vary in accordance with the 
size of the system and the with the number of operations to be performed. 


The KST~1 system has the following specifications; 
Range (bimetallic two-wire line suspended in air, wires 4 mm in diameter without intermediate units) 500 
Frequency range for calling and specifying operations, 70-300 cps; 


Time to call one object 0,2-0.3 sec; 
Capacity (number of objects connected in one direction to a two-wire line) up to 100 *; 


* The capacity can be raised by changing the frequency combinations; it {s limited by the loading conditions 
at the control desk. 
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Fig. 4. General view of the local equipment, 1) Calling unit; 2) control unit; 
3) telemetry transmitter. 


Number of control or adjustment operations, 2* ; 

Number of signals available per object, 2°; 

Number of parameters that can be measured per object 2°; 
Time for one measurement cycle, 2 sec; 

Reduced telemetry error (neglecting transducer error), 1%; 
Repeaters and teeing equipment can be fitted; 

Supply to control desk, 220 v; 

No supplies needed at objects. 


* The number of telemetry, control, etc. signals can be increased by introducing appropriate extra frequencies 
for the operations required. 
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Figures 3 and 4 show photographs of the desk and of a local unit in the KST-1 system, 


A prototype installation was tested in the laboratory, and then on the Tula~Moscow gas main, The range 
was found to be 230 km with an fron circuit, equivalent to more than 500 km with a bimetallic one, 


The Union Office of Orgenergoneft', in collaboration with the Institute of Automation and Remote Control, 
developed the apparatus, made the apparatus for the Kazan-Gorky and other gas mains, installed the equipment, 


and tested it, 
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A NEW ELECTRICAL PROPORTIONAL-ACTION EFFECTOR 


V. 1, Litvak and A. Ya, Shapiro 
(Kiev) 


Proportional control, in which to every value of the controlled parameter (density, temperature, pressure, 
humidity, level, etc.) there corresponds a definite position of the controller, has sometimes to be used with auto- 
matic control systems, The Central Design Office of the Ukrglavpromdmashdetal' has designed a new electrical 

_ proportional-action effector (type IM) which is intended to work with various control heads and {s meant for use 
with manual remote control, 


The Main Characteristics of the IM Effector 


The torque at the output shaft is 3.5 kg-m; one complete movement of the output shaft (through 90°) takes 
150 sec; maximum movement of the output arm 42 mm; resistance of feedback potentiometer 50 ohms; power 
consumed about 35 w; size 348 x 247x 150 mm; weight 5.9 kg. 


The unit is meant for repeated short-term operation. The supply {is 20 v at 50 cpa. 


Principle and Circuit 


The IM effector consists of a reversible electric motor R, a reduction gear G, a balanced relay B, a feedback 
potentiometer P, limit switches K and an output unit O (Fig. 1). The motor R is asynchronous (capacitor). It {s 


of the short-circuited rotor type, with. a shorting ring and four wound poles on the stator, The sense of rotation is 
determined by the pair of windings to which the matching transformer T and phase~shift capacitor C are connected. 
The sense of rotation is altered by switching the supply from terminals I and II to terminals I and III. 


The motor drives the effector via the gear G; the output unit consists of a crank mechanism and shaft, and 
of the feedback potentiometer P. 


The balanced relay B alters the sense of rotation of the motor. 


The relay has two coils K; and Ky, each double-wound (main and auxiliary windings) and a core Cr, which 
is a horse~shoe of low-carbon steel and is coupled mechanically to the contact plate CP. Core Cr rotates as Ky 


| 
Fig, 1. 


View with rear cover removed 
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View on A 
(front cover removed) 


Fig. 3. 
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and Ky; dictate, The supply voltage is fed via the con- 
tact plate either to terminals I and II or to I and III, 
The auxiliary windings make the operation more stable; 
they are wired in series with the stator coils, They in- 
crease the attraction exerted on the core, which is im- 
portant to stable operation if Ky exceeds by a little 

the force of Kyy, or vice versa, 


The potentiometer slider is fitted with two cams, 
which press on the limit switches K when the output 
unit O reaches its extreme positions. This unit O is de- 
signed to work with controllers that operate by recipro~ 
cation, 


The IM works in conjunction with automatic devices 
that have error transducers, The transducer may be a 
potentiometer, an induction device, or some others de- 
vice that gives an output voltage which {s linearly re- 
lated to the error, The slider on the transducer poten- 
tlometer, Pz, is displaced if the controlled parameter 
is in error; the circuit composed of this potentiometer 
and of the feedback one, P, then becomes unbalanced, 
Coils Ky and Kj, then receive unequal currents, so the 
relay operates, and the supply voltage is fed via the 
contact plate either to terminals I and II or to I and Il, 
The motor runs until the balance condition has been 
restored, 


* The operation is analogous if the automatic regu- 
lator uses an induction transducer. 


Design of the Effector 


The asynchronous motor consists of the stator 20 (Fig. 2), the rotor 19, the capacitor 13 and the transformer4, 


The stator core is laminated, The four stator coils are fitted to the poles, The rotor spins in two radial ball- 
bearings, which are fitted in the body 3 and in the end-plate 21. The shaft bears a coupling bush that transmits 
the motion to the reduction gear 18, 


A standard capacitor type KMBG is fitted, of capacitance 1 yf, designed for a maximum working voltage 
of 600 v. It is fitted in a special recess in the body 3, and is held by a spring. Transformer 12 is a step-up one 
(ratio 12). 


The reduction gear 18 is fitted inside the effector; it has six pairs of straight-cut gears, The gear that en- 
gages with the motor's shaft is made of textolite, to silence the reduction gear. The output gear isin the form of 
a sectar, fixed to the shaft 2 of the crank mechanism. 


The balanced relay is made up of two coils 15 and a core 10, with the contact plate 11 fixed on a plastic 
strip 1, which is fixed to the motor's end-plate 21. Strip 1 also carries the fixed contacts 16 of the relay, the con- 
tacts 8 of the limit switches, and the feedback potentiometer 14 (a textolite body with a constantan winding), The 
slider 5 is fixed to the shaft 2 of the crank mechanism by a grub screw. 


The limit switches consist of the contact plates 7, of springs, and of fixed contacts 8, These switches open 
when the cam on the slider 5 presses on the contact plate, and thereby break the motor's circuit, The plate ts re- 
stored to its initial position by the spring when the cam withdraws. 


The body 3 carries the motor and the reduction gear 18; the limit switches, and the strip which holds the 
balanced relay and the feedback potentiometer are enclosed in the section 34 with end-plate 24 and 26 (Fig. 3). 


Fig. 5. 
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A peg holds the crank 28 onto the shaft 2. A special bolt 30, 
fitted with a spacer 29 and a nut, holds the crank onto the rod 31 that 
drives the effector shaft 39, The distance at which this rod is fixed 
from the center of the crank can be adjusted, to alter the travel of the 
effector shaft. The crank is graduated, The setting of the required 
travel is made easy by bringing the shaft out through a housing 36 that 
contains a spring shock-absorber 9; this shock-absorber prevents the 
effector's fracturing should the moving parts seize, or be locked by, some 
foreign body. 


The length of the output shaft 39 that projects from housing 36 
can be adjusted, as is necessary when theunit is being fitted to various 
control units, The length is adjusted with the nut 17, which is fitted 
in the hexagonal projection on the housing 36, The shaft is held by a 
grub screw, which can be reached via a hole in 36 which-is normally 
closed by a plug 12, 


The cutput shaft 39 has a wider section at the end, in which there 
is a hole that connects the shaft to the throttle, etc. The peg 38 and 
ring clip 37 hold the parts together. The crank mechanism and output 
shaft are housed in the case 27, which has a bush for the output shaft, 
The shaft is oiled via the oiler 33, The lower part of body 27 has a 
flange, which is used to fasten the effector to the regulator, Figure 4 
shows the effector coupled to a valve. 


The front cover 22 has a hole through which a shaft is brought out to the manual controlhandle 32. This 


shaft carries a ball clutch and spring 6, which protects the mechanism from breakage. This clutch is so designed 
that the force needed to causeit to slip can be adjusted, 


The side of the body 34 has a recess that is fitted with a terminal block, through twhich the supply voltage 
and transducer are connected, The leads are taken in through thebushes 35, The terminal block is covered by 
the cover 40, 


Figure 5 shows a photograph of the mechanism. 


This effector is simple in design, is light, gives a fairly large operating force while of small size, and has 
a large travel (up to 42 mm) on the output shaft, It is cheap to use and has a long life, The first batch of these 
effector has been made by the Kiev Control and Measuring Instrument Factory. 
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DISCUSSIONS 


THE PRESENT AND THE FUTURE OF AUTOMATIC CONTROL 
THEORY AS A SCIENCE 


(Discussion in the Institute of Automation and Remote Control of the 
Academy of Sciences, USSR) 


From the editorial board, A discussion on the subject "The Present and the Future of Automatic Con- 
trol Theory as a Science” was held on March 12, 1958 at the seminar on nonlinear problems of the automatic 
control theory. The discussion attracted some 250 people from various scientific research institutes, development 
laboratories, and enterprises of Moscow. Fifteen people took part in the discussion expressing various points of 
view on the present state and lines of development of automatic control theory as a science. In publishing 
some of the contributions to the discussion, the editorial board requests the readers of this journal to take part in 
the discussion and send their remarks on the subject to the journal, 


M. A. Aizerman's Contribution 


Whenever we reflect on the destiny of our science, we inevitably raise, side by side with purely utilitarian 
problems, such important questions of principle as: “What is the theory of control? What does it study? What is 
its aim? Where is its place among other scientific disciplines? * 


When examining monographs on the automatic control theory,one is compelled to think that the publication 
of these monographs and textbooks is due to the necessity of training and retraining in the art of control, engineers 
who are not acquainted with a number of required scientific disciplines, 


Let us try to exclude from these textbooks everything which refers not only to automatic control systems, 
but also to a wide range of dynamic circuits (for instance, circuits used in the transmission theory, radiotechnical 
devices, etc.). Then, the only thing which will remain of our automatic control theory courses, i, e., what is 
specific to automation and control, will easily fill a few dozen pages, One naturally comes to the conclusion that 
the term “automatic control theory” in the sense it is now understood is due to the manner of studying, rather than 
to the development of science, and is at present justified only by the necessity of training and retraining specialists, 
who do not possess the required grounding in bordering disciplines. 


It is not easy to answer the question: “What is the object of investigation in the automatic control theory?" 
To say it is “automatic regulators" or “automatic control systems" is equivalent to saying nothing, What is the 
difference between an automatic control system and any other dynamic system? Why is the knowledge of general 
design methods of dynamic systems insufficient for designing automatic control systems or for understanding their 
operation, and why is a specific special automatic theory required? I think these questions arise precisely be~ 
cause by “a theory of automatic control” we do not understand the sum of scientific problems specific to control, 
and distinct from problems which arise in other technical spheres, but a study course which should be more accu- 
rately called “development and design of automatic regulators and control systems." Is it possible that automa- 
tion, i. e., automatic tracking and control,have no specific, scientific problems of their own, distinct from those 
in the remaining technical spheres? I do not think this is the case, The objective of automation does not lie in 
the sphere of dynamic calculations of automatic control systems, since these calculations differ but little from 
dynamic calculations in other technical spheres, Its objective lies in studying the specific laws and properties 
of control processes, Controlling production, a process or an object involves finding the best method of placing 
the control organs for the achievement of this aim. It is not often that a simple connection of a regulator will 
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solve the problem, The problem of finding the correct method of placing the control elements consists,in the first 
instance, in choosing an optimum method of “conducting” the process, of controlling it. 


During many decades the basic automation problem was put as follows: there is given an object which is 
fully specified,including the element which could serve to control it,and the method of automation is obvious, 
i, e., the method of connecting of a regulator or regulators, 


The design problem under such conditions became a problem of dynamics; how to find regulator parameters 
which would make the system “regulator-object" satisfy certain predetermined requirements, Such a posing of 
the problem is only specific for the automation of relatively simple objects by means of individual regulators which 
are specially made for these objects, At present, the problem is often posed in a similar manner which makes the 
dynamic calculations predominant, For anyone, however, who observes the development of regulator manufacture, 
it is clear that a different posing of the problem is beginning to predominate, 


A new branch of industry, that of building regulators, has developed, which supplies customers with regula- 
tors of increasingly universal application. The people designing regulators do not and cannot know where, under 
what conditions, with what objects in view and in conjunction with what type of other regulators each regulation 
must operate, Every regulator is made with a very large range of adjustment. A basic assumption is made that 
the regulator and the object will not and cannot be calculated, and that a suitable adjustment can easily be found. 
Who, it may be asked, can design a control system under such conditions? The designer of the regulator? He 
cannot do it, since he does not know the intended object of regulation or the conditions under which the regulator 
is going to be used, The purchaser of the regulator? He need not be, however, a control specialist and need not 
know the details of the regulator design and less still the equations controlling its operation, To him a regulator 
isa “black box” with definite factory guaranteed possibilities and with more or less convenient adjusting elements. 
The choice of an optimum adjustment can under such conditions require certain calculations, but of the most sim- 
ple kind, Moreoever, a self-adjusting regulator is not a dream but a thing of the near future, 


This tendency is developing, The industry is changing over from the production of universal regulators to 
the manufacture and placing on the market of “regulator units," separate self-contained elements, from which a 
customer, unknown to the designer, can build up most diverse regulators and regulating systems, 


The present development organization for production automation hardly ever carries out dynamic calcula- 
tions,not because it is incapable of making them, but because they are seldom required, The dynamics of the sys- 
tem seldom worry the designer, It is another range of questions which worries him, How should a technological 
process be controlled? Which parameters should be measured and which should be controlled? How should a cir- 
cuit be designed? The design of the automation schemes is as yet an art, Nevertheless, there certainly exist and 
there can be established scientific foundations of automation proper, scientific methods of approach to control 
problems, to circuit design. The history of the relay-contact circuits, whose design remained for a long time an 
art, but can now be entirely based on scientific foundations, inspires confidence in the attainment of this aim. In 
my opinion, it is precisely this range of questions, specific to automatic control and completely different from the 
problems studied in other disciplines, which will in the future comprise the main content of the theoretical foun- 
dations of automation as a science, Whereas the problems which at present form the basis of the automatic control 
theory will gradually stop being considered as problems of the theory of control and will become subjects ofinves~ 


tigation of other disciplines: the general theory of dynamic circuits, the theory of oscillations and the theory 
of stability. 


Then the theory of automatic regulation will gradually develop from a study course into an independent 
science, a science of control and the means and methods of automation, 


All this, of course, does not in any way decrease the importance of study courses whose object is to provide 
the automation specialist with a combined knowledge of boundary disciplines (theories of oscillations, stability, 
differential equations, etc.) which are required by a modern designer of regulating devices, 


Neither does this decrease the importance of dynamic calculations, They provide the understanding of 
many laws controlling the regulating processes, i. e., facilitate the accumulation of facts concerning the control 
processes. They are indispensable for designing universal regulators and regulating unis unconnected with the ob- 
ject of control, since these regulators and units contain internally closed circuits, Finally there will always re- 
main individually constructed regulators and,hence,the classical problem of avestigating the closed system of a 


“given object-regulator.” Therefore, it will always be useful and important for regulator construction to develop 
the theories of dynamic circuits, oscillations, stability, and several branches of mathematics, etc, 


We cannot ignore, however, the direction in which the regulation technology is developing, since it deter~ 
mines the changes in the meaning we ascribe to the term of the “theory of automatic control." 


N. N, Moiseev's Contribution 


Nearly all the previous speakers have, in my opinion, correctly characterized the present condition of the 
theory of automatic control. In the “heroic” period of its development, the automatic control theory adapted 
various methods of bordering disciplines for the solution of concrete problems connected in one way or another 
with production technique. The basic problems of the automatic control theory are so closely connected with 
problems of bordering disciplines (mathematics, mechanics, the theory of circuits, etc.) that they can be rightly 
considered as belonging to these disciplines, Therefore, as it was pointed out, it is very difficult to separate ques- 
tions which are peculiar only to the theory of automatic control, Recently,completely new lines of development 
arose, peculiar to the control theory only, which developed as the result of natural practical requirements, 


I should like to make two remarks with respect to these new problems. 


1, In the literature on the automatic control theory,there appeared in the last decade an enormous number 
of papers which in fact belong to the theory of differential equations. Many concrete linear and nonlinear pro- 
blems were studied and often refined quantitative methods developed. This traditional development {s very fruit- 
ful and has not exhausted its possibilities by a long way, A complete solution of many important problems was 
obtained, thanks to this method, Many dynamic properties of nonlinear systems are as yet completely unknown 
to us. Nevertheless, this qualitative aspect of the theory by itself is obviously insufficient. The present state of 
the theory of differential equations is such that it is necessary for the solution of concrete problems to simplify 
the system to the limit. It is necessary, therefore, to resort to methods of numerical analysis. The availability 
of rapid electronic computers brings this problem to the foreground, since their use opens up completely fantastic 
prospects, It is possible to say, of course, that these problems are the concern of mathematics and not the theory 
of automatic control. Yet the level of development of the modern calculating technique does not as yet provide 
any standard methods of calculation. There are even quite simple problems whose solution encounters almost 
insurmountable difficulties. As an example, it is sufficient to quote the border problems of differential operators 
of the elliptic type. Among the problems of the theory of automatic control systems with distributed parameters 
(for instance, a system with hydraulic elements), there arises the necessity of solving Neumann's problem in La~ 
place equations, The numerical solution of these problems will require,in certain instances, the development of 
special methods, 


Obviously, all these questions must become the subject matter of a special science, that of numerical cal- 
culations, This branch of mathematics, however, is in its infancy and the elaboration of concrete numerical cal- 
culation methods for the solution of various practical problems is being carried out as yet only in bordering theoret- 
ical disciplines: in gas hydrodynamics, mechanics, etc. 


It is, by the way, sometimes simpler for us than for the mathematicians to invent such a numerical method, 
The point is that it is often possible to take into account the physical peculiarities of the problem. In this we are 
greatly assisted by intuition and experience. Therefore, such difficult problems as finding periodic conditions in 
complex nonlinear systems with a large (or infinite) number of degrees of freedom by means of rapid computers 
must be kept in mind by the specialists in the theory of automatic conttol. 


2. It has already been pointed out in this discussion that one of the new problems is arising in optimum 
control, I should like to note that such a condition is not always required and the very posing of the problem of 
finding an optimum method of application isan academicone, There exist, however, problems where the ques- 
tion of optimum control is reasonable. In this connection, I should like to draw your attention to the remarkable 
investigations of L. S. Pontryagin and his pupils who are beginning to publish their work, To me it seems especially 
important that their method reduces the finding of optimum control to the solution of ordinary differential equa- 
tions by Cauchy's method. The possibility of *mechanizing™ this problem provides a basis for arriving ,in many 
cases, at most valuable practical conclusions, It should be noted, of course, that in order to attain this aim, 
specialists in the automatic control theory will have yet to overcome many difficulties. 
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At present there is a number of mathematical problems arising in the midst of the automatic control theory. 
A successful solution of these problems, it would appear, is only possible if the physical pecularities of the auto- 
matic control systems are made use of, Therefore, such problems, in my opinion, should be considered as pertain- 
ing to the automatic control theory, The participation of mathematical specialists in the solution of these pro~- 
blems is necessary, but the main burden, in any case at first, must probably fall on spectalists in the theory of 
automatic control, 


v. V. Nemytskii's Contribution 


Not being a specialist in the sphere of automatic control, I cannot, of course, point to any-lines of develop- 
ment of this science, I can only speak as a mathematician and specialist in the sphere of differential equations, 
I agree with the opinion of some of the previous speakers that the part played by the theory of differential equa- 
tions in the theory of automatic control has been considerably exaggerated and that the time has now arrived to 
apply, for the solution of automatic control problems,other mathematical devices. 


On the other hand, however, the possibilities of applying the qualitative theory of differential equations are 
far from being exhausted, This applies,in the first place,to the study of the behavior of multidimensional dynamic 
systems, It should be remembered that many of the problems of automatic control have been historically connec- 
ted with problems of radiophysics, which deals in the main, with dynamic systems with one degree of freedom, 
Concepts and types of movement, for instance, the representation of autooscillations as limiting values of Poincaré 
cycles, are closely connected precisely with the behavior of systems with one degree of freedom. 


In multidimensional systems, especially for forced oscillations, completely different types of movement can 
arise,as shown by the qualitative theory of differential equations, It appears to me, however, that the automatic 
control theory specialists are not taking into consideration the possibilities this circumstance presents and are not 
investigating the part played by these complex movements in systems encountered in practice. In the solution of 
all these problems, the theory of differential equations will yet render valuable service to the theory of automatic 
control, 


D. I. Mar'yanovskii's Contribution 


The theory of automatic control is an engineering science whose objects of investigation are actual regula- 
tors and systems of automatic control, used in industry and scientific research, 


The task of the automatic control theory is to study the general laws governing the behavior of these objects. 


It is reasonable to include in the automatic control theory only matter which provides an understanding or 
deeper insight into the processes taking place in regulators or control systems,or assists the engineer in designing 
them. The criterion for such a classification is convincing,practical applicability of this matter. 


In order to make above general indications of the tasks of the automatic control theory more concrete, a 
few examples from the range of problems which should be the subject of this science will be quoted, 


1. Examination of elements which perform the same functions in different systems. 


2. Comparison of a wide range of different types of systems,with the objectiveof finding common features 
and differences inherent in them, Evolving general concepts for such types of systems and requirements for the 
transient and steady-state processes in them. 


3. Evolving common features and peculiarities of control systems in certain technical spheres, Evolving 
for these spheres of common automatic control system parameters, having definite effects on the controlling pro- 
cess, 


Is it advisable to limit the automatic control theory by the range of problems outlined above, i, e., by the 
study of existing regulators and systems,or is it better to extend it to the investigation and development of a theory 
for any regulators and systems satisfying definite technical requirements encountered in practice? 


There is insufficient material available on this subject and it is difficult as yet to answer it categorically, 


It is possible to cite,as an example,an interesting idea which, as it appeared at first, could be widely applied 
and which attracted,there fore,considerable attention, 
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We are referring to the so-called extremes regulator, whose principle of construction in a general form and 
for controlling thermal process was proposed in 1940, 


Despite the fact that successful designs embodying that principle were proposed which would provide its 
effective application in practice, users of thermal regulators did not display any interest toward it, and, as far 
as I know, with the exception of two or three experimental models, not a single thermal regulator of this type was 
ever put into use. 


It would appear that the principle of the extremes regulator has not been applied in other branches of in- 
dustry either, 


Despite the fact that extremes regulators are not as yet used in a single branch of industry, the theory of 
extremes control, judging by the number of publications, is being intensively developed, It is obvious that the 
expenditure of efforts on establishing a theory of such a "phantom regulator” cannot be as yet justified. 


In sciences which use extensively mathematical devices, as for instance, in astronomy and physics, there 
are instances when important scientific discoveries are made on a purely theoretical basis. 


Hence, there can obviously arise a question: Can the automatic control theory as a branch of science which 
uses mathematical devices extensively, provide new types of regulators and principles of control? 


It should be noted that the appearance of regulators and control systems based on new principles is usually 
due to the necessity of solving problems which arise in the practical activity of men, taken in the narrow sense of 
“practical.” Despite the attractiveness of the ideas evolved by the automatic control theory, one should follow the 
path of speculative creation of new systems and regulators very cautiously in order not to fall into visionarism, 
Without excluding in principle the possibility of discovering, by purely theoretical means, new regulators and con- 
trol systems which can be used in industry, it should be considered that such a possibility is remote. 


Let us now examine to what extent the existing monographs and courses of automatic control theory reflect 
the applied character of this science and to what extent they can satisfy the requirements of those designing, pro- 


ducing and employing new regulators and control systems, 


The majority of monographs and textbooks on the theory of automatic control usually begin with a chapter, 
‘comprising not more than 5-10% of the book, which contains descriptions of control systems, The remaining part 
of the book is usually filled with mathematical calculations, which are elementary or otherwise, depending on the 
author's level of mathematical training and the content of his dissertations, The greater part of these calculations 
is not and cannot be illustrated by concrete examples, The reader waits with impatience for the explanation of 
the relation between the mathematical calculations and the promising beginning of the book, He waits for it in 
vain to the very end of the book which he slams with disappointment, convinced that he has wasted the greater 
part of his time. 


The useless material consisting of theories and methods which have not justified themselves in the sense 
that they provide no concrete instructions for designing control devices, and do not help to understand the pro- 
cesses occurring in them. 


The most characteristic example of such a theory {s the so-called *method of small parameters” which was 
proved to be completely barren. Not a single practical control problem has as yet been solved by this method, It 
seems to me that the inclusion of such methods in books dealing with the automatic control theory considerably 
lowers the authority of that science, 


The inclusion of various graphic, analytical and analyticographic methods of approximate calculations and 
the plotting of curves of transient processes in automatic control systems seems equally pointless, 


The proposed marketing of portable electronic rapid computers will soon make them as commonplace as 
equipment in development laboratories and scientific research institutes as, for instance, the photostating equip- 
ment is at present, 


The advantage of such computers as compared with the calculating methods of transient processes described 
in the book is obvious. 


Current practice requires the solution of such problems as; How and in what direction should certain param~- 
eters of a system be changed or, how should the structure of one of its part be changed in order to obtain decreased 
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overshooting or an increased controlled-quantity speed of variations without overshooting after a disturbing effect, 
or finally, how to obtain an aperiodic approach of the controlled quantity to its stable~state value after one or two 
overshots, etc, The solutions of such problems must be qualitative, and must not be obtained by calculating the~ 
transient process for several numerical values of a certain parameter,i, e,, they must be essentially experimental, 


Side by side with the overloading of books on the automatic control theory with useless material, some of 
the design problems important for the engineer, are not dealt with, among which should be included the stable- 


state system calculations, for instance, the calculation of the gain of system elements containing electrical de~ 
vices. 


What then are the future lines of development of the automatic control theory? These lines of development 
should be determined by certain measures whose execution should be assumed by the Institute of Automation and 
Remote Control of the Academy of Sciences, USSR, 


1. In cooperation with the workers in industry, development laboratories and industrial scientific research 
institutes, a plan of theoretical development should be evolved. The objectives of the plan in the form of con- 
crete problems should be distributed among post-graduate students and scientists with the proviso that the solution 
of these problems will entitle them to corresponding scientific degrees and monetary awards, 


The approval of the solutions of these problems should be entrusted to the persons or agencies which raised 
them in the first place. 


Problems solved in the form of papers should be printed in the journal “Automation and Remote Control" 
immediately in preference to other material, 


2. A plan forissuing a series of monographs of 4-10 printer's sheets should be established. Each monograph 
should be written by at least two authors (a theoretician and a practical engineer), The monographs should con- 
tain an explanation of the general principles governing the operation of automatic control systems used in definite 
spheres of industry, or establish general conceptions, properties and phenomena common to a type of control system 
classified by some common property. 


According to available information, Energoizdat (State Power Press) will shortly publish a series of small- 
size monographs on automation, This edition should be made use of,in order to provide the monographs with a 
content corresponding to the aims set above, 


3. No further editions of new basic textbooks or monographs on the automatic control theory should be issued 
or old ones reissued until new material is accumulated which,in conjunction with the valuable material contained 


in the existing books,would provide an automatic control theory worthy of our epoch of great scientific achieve- 
ments, 


P, I, Stakhovskii's Contribution 


At present a very promising type of control system is being developed — the self-adjusting system, One 
such type of system is that of automatic optimization, which is different in principle from the normal automatic 
control systems. 


In fact, the aim of an automatic control system is the elimination of the unbalance between the actual and 
the set values of coordinates, Thus the actual value of a coordinate provides the required change in the input sig- 
nal for the restoration of the unbalance, 


An automatic optimization system must provide a combination of input signals Xy(i= 1, . . ., m) which would 
give the required quantity Q an extreme value, and that obtains (if the function of Q has a first derivative at all 
the points of the required region of the function) when 
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Moreoever, neither the type of function Q(X;, . . ., X,) nor the value of Q.., are known inadvance, Only 
the conditions which function Q must satisfy in a small region around the required point are known, In view of 
these circumstances, in order to predict the movement of the system it is necessary to know the behavior of the 
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function of Q in a small region around the point, determined by the current values of X;,..., Xp, and rationally 
to extrapolate the movement of the system on the basis of previously obtained data by investigating a small area 
of the current values of coordinates, Therefore, the peculiarities of an automatic optimization system consist in 
small fractional movements, the analysis of the results of tentative movements and the prediction of the most ad- 
vantageous further movement of the system. 


A very important case in automatic searching is that of searching in restricted conditions, 1, e., the case 
when the region of possible combinations X;, . . ., Xp is limited by the inequality Hy (Xy,.. Xn) SO (J =1, 
. + «,m™). In other words, a condition is made that in searching the extreme value of Q,other functions of the same 
variables should not exceed in their numerical value certain given quantities. Moreover, the optimum value of 
Q can lie on the border of the region of permissible variations of X;, ..., Xp and the behavior of Q at the point 
of the nominal extreme (optimum) must satisfy the relationship [grad Qx grad Ho] = 0, where Hy = DHj, providing 
the sum includes only the values of Hj, for which at the optimum point equation Hj (Xj, . . ., Xp) = 0 holds. 


The objective of the automatic optimization theory is the finding of the most advantageous movement 
algorithms which lead to the determination of the optimum value in the quickest and most accurate manner, For 
this purpose it is necessary to investigate the movement of the system in the process of finding Qopr (transient 
searching process) and the stable~state movement of the system in the region around the determined optimum value, 
It is necessary to compare the searching process indexes (the time of searching, accuracy of determining Qopt, and 
other indexes) obtained by means of different methods and with different types of function Q(X,, .. ., Xn) with 
the aim of determining the optimum algorithm, The latter can differ fo- different types of functions of Q, It is 
important to determine the searching process conditions of stability. 


The effect of random interference on both the transient and stable-state searching processes follows from 
the principle of operation of automatic optimization systems, since the prediction of the movement of the system 
is determinedby small changes in the input signal, The analysis of the movement of systems subjected to random 
interference and the synthesis of noise~proof systems is, therefore, a very important task of the automatic optimiza- 
tion theory. 


Thus, the theory of automatic optimization and that of self-adjusting systems includes in general a wide 
field of investigation. 


I should like to support the opinion expressed by several speakers before me that the further development 
of automatic control will consist in the development of the theory of self-adjusting systems, 
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CHRONICLE 


GRIGORII MITROFANOVICH ZHDANOV 
(ON HIS 60TH BIRTHDAY AND 30 YEARS OF SCIENTIFIC AND TEACHING ACTIVITY) 
The Moscow Order of Lenin Power Institute has celebrated the 60th birthday and 30 years of scientific and 


teaching activity of the prominent scientist Prof, Grigorii Mitrofanovich Zhdanov, responsible for training specialists 
in the sphere of automation, remote control, and calculating tedhniques, 


A hard life preceded G, M. Zhdanov's scientific activity. By 
combining study with work, he obtained his secondary and part of 
his higher education, In 1919 he was called up to the Red Army 
where he served as a Red Army man, political worker, and teacher 
of social sciences, From 1928 to 1932, G. M. Zhdanov studied at 
the Moscow University of National Economy and the Moscow Power 
Institute, He remained there for post-graduate studies and as early 
as autumn of 1932 he began to lecture on the theoretical foundations 
of electrical technology. 


When his candidate's dissertation was accepted in 1937, G. M. 
Zhdanov was appointed to the recently created chair of Automation 
and Remote Control, From then on, G. M. Zhdamnoy led a life of a 
tireless organizer, teacher, scientific and social worker, devoting 
all his energy to the task of creating and developing new branches 
of technology. He developed the study of and lectured on remote 
control, remote measuring, automatic control, and remote control 
equipment, foundations of automatics, foundations of remote con- 
trol, automation and remote control of power systems, computer technique, and mathematical analog computers, 


In theyears 1951-1953, G. M. Zhdanov established in the Moscow Power Institute a faculty and a laboratory 
of mathematical computers which he heads at present. In 1954 G, M, Zhdanoy obtained his doctor's degree and 
in 1955 he was awarded a professorship. 


Many young specialists, who are now successfully working in various branches of automation, remote con- 
trol, computer technique and electrical instrument construction have been trained under the guidance of G, M. 
Zhdanov. He has also trained a number of teachers, scientific workers and aspirants. 


G. M. Zhdanov is responsible for 44 scientific works and inventions, His two volume work “Remote Mea- 
surements” (1952-1953) remains the foremost systematic textbook of its kind, Many of his works (1953-1956) con- 
sist of basic textbooks on mathematical analog computers, G. M. Zhdanov created in the Moscow Power Institute 
a calculating center and put in hand the organization of an automation problem laboratory. 


In 1955, G. M. Zhdanov, organized and became head of the section of automation, ermote control, calcula- 
ting and measuring technoldgy of the journal *Elektrotekhnika“ of the All-Union Institute of Scientific and Tech- 
nological Information. 


G. M, Zhdanov is a member of the Scientific and Technical Council of the Ministry of Higher Education 
and chairman of the Instrument Section of the Council, he is a member and deputy chairman of the Scientific 
Methodological Commission of the Central Administration of Polytechnical and Engineering Higher Education 
Establishments, memebr of the Moscow Power Institute Scientific Council, member of the editorial board of the 
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journal “Higher School Papers" (series on electromechanics and automation), and member of the Scientific Tech- 
nical Council of the NIlschetmash (Computers Scientific Research Institute). 


G. M, Zhdanov's extensive and fruitful work was rewarded by a number of diplomas, three medals and the 
Order of the “Badge of Honor." 


Numerous colleagues and pupils of G. M. Zhdanov have expressed to him on his birthday their heartfelt 
congratulations and wishes of good health and continued successful activity. 


G, Chilikin and R. G, Romanoy 
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REVIEWS AND BIBLIOGRAPHY 


REVIEW OF YA. Z. TSYPKIN‘'S BOOK “""THE THEORY OF 
PULSE SYSTEMS" 


Ya. Z. Tsypkin's monograph "The Theory of Pulse Systems" (Fizmatgiz, 1958) contains an exhaustive de- 
scription of the theory of discrete acting dynamic systems — pulse systems. This theory has been developed by 
inany investigators in various countries, among whom one of the leading positions is occupied by the author of 
the monograph himself. It should be noted that there are no monographs on the subject with the exception of the 
two written by Ya. Z. Tsypkin (the one under review and the other entitled "Transient and Stable-State Processes 
in Pulse Circuits," Energoizdat, 1951). Yet pulse systems are increasingly widely used in most diverse branches 
of industry ; hence, a synthesis and full consistent description of the theory of systems of this class is very important 
for practical purposes. From this point of view the appearance of the monograph "The Theory of Pulse Systems," 
written in a rigorous yet understandable style for a wide circle of scientific workers, engineers and students, is 
most opportune. 


The monograph deals in the main with the theory of line pulse systems with constant parameters, which is 
exhaustively described on the basis of the author's work. Attention is also paid in the book to variable parameter 
pulse systerns and nonlinear pulse systems. The basic principles of the theory of these systems are given and some 
particular designs valuable from the practical point of view are quoted. 


The monograph describes both pulse systems without feedback, which include many radiotechnical circuits, 
and those with feedback, which include various systems of interrupted control and tracking of different objects. 


Before describing the theory of pulse systems, the author starts with a short introduction and then in Chapter I 
outlines various types of pulse systems. In particular, he introduces the classification into three types of systerns 
(systems with AIM, ShIM and VIM). He deals with the possible nonlinear characteristics of pulse system eleme nts 
and outlines closed and open pulse systems. Examples are given mainly from electrical systems, but some refer 
to the mechanical systems as well. It should be noted that most examples have been very well selected. The 
author describes very modern systems which engineers of various specialties constantly encounter in practice. 
Such systems include synchronous filters, circuits with a periodic commutation, numerical computers, etc. It 
should, however, be noted that examples referring to radio location have not been described very satisfactorily. 
They serve more as illustrations instead of being examples of modern circuits usually employed in radio location. 


The many examples quoted in Chapter I makes the description more lively and the application of the theory 
more graphic. 


Chapter II deals with the mathematical basis of the pulse system theory in particular, with difference equa- 
tions and discrete Laplace transformations. The author introduces the concept of a lattice function, in our opinion, 
a happier one than his previous notion of a step function. The properties of lattice functions and the equations 
which they satisfy (difference equations) are investigated. For the solution of difference equations, discrete La- 
place transforms are used, and the basic theorems of the operational calculus on the basis of Laplace and Fourier 
transforms are given. The relations between discrete Laplace transforms and the FourierandLament series and 
between the representations of a continuous function and the corresponding lattice function are given. The find- 
ing of these relations is very important since it facilitates the understanding of the essence of the theory. 


The mathematical basis of the theory related to discrete transforms is described sufficiently fully and convinc- 


ingly. It should be, however, noted that in the literature cited by the author at the end of the book there are 
foreign authors who use somewhat different methods of developing the pulse system theory, which are not directly 
telated to the use of discrete transformations. Considering the exhaustive character of the monograph, it would 
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be. very desirable to include a review of the various approaches to the analysis of pulse systems and point out the 
relation between all these methods. This would greatly facilitate the choice of a method for solving problems 
which would produce the required answer in the speediest manner. In particular, the works of Linville, Ragaz~- 
zini, and Zade. deserve attention. Chapter II also contains information on random lattice functions and their 
characteristics important for the understanding of the theory. 


Chapter III deals with the foundations of the theory of pulse systems without a feedback (open loop systems), 
It should be noted that many of the propositions developed in this chapter apply equally well to closed systems, 
In the case of linear pulse systems with constant parameters, however, the division of the theory into two parts 
which refer, respectively, to open and closed loop systems is justified, since the properties of the closed loop sys- 
tems are subsequently expressed in terms of the open loop systems. Such a separation of the theory loses its pur- 
pose to a considerable extent ,for variable parameter or nonlinear pulse systems. The propositions of the variable 
parameter pulse systems quoted in Chapter III also hold, with the exception of a few particular formulas, both in 
the case of closed and open loop systems. 


The theory of open loop systems, starting with a mathematical description, is developed into a unified form 
convenient for further analysis. By introducing further the concept of pulse reaction, the author makes use of the 
superposition formula and with the aid of Laplace discrete transformations introduces the concept of the transfer 
function of a pulse system. Formulas for the transfer function valid for various conditions are derived and its pro- 
perties described. The important case of pulse systems with several pulse elements is described and a transfer 
function for such systems derived. 


Basic relationships connected with the concept of a variable transfer function are derived for systems with 
variable parameters and certain particular cases of parameter variations are examined. 


In systems with constant parameters, in addition to transfer functions, frequency and time characteristics 
are also analyzed, and practical methods of plotting them are described. The concepts of pulse system charac- 
teristics permitted the author to formulate means of finding the system's reaction to an arbitrary input signal. The 
examination of random processes in pulse systems, given in the book, is of great practical value, since in many 
modern problems, it is either necessary to examine systems in the presence of random signals,or to deal with ran- 
dom variations of the system parameters. The chapter concludes with a comparison of open loop pulse systems 
and continuous action systems. Thus, Chapter III deals very fully with the foundations of the theory of pulse sys- 
tem which have no pulse operated feedback. 


In Chapter IV, the author gives many examples in which he applies the theory he has developed for analyzing 
various aspects of the processes in mechanical, electrical, and electromechanical systems. The examples cited 
are up to date as in Chapter 1. The solutions of the problems given in the book serve to illustrate the open loop 
pulse system theory and at the same time are valuable in themselves. 


Chapter V deals with the theory of pulse systems which have a pulse feedback. First, the equations which de- 
scribe these systems are reduced to a unified convenient form, then by means of the superposition method the ex- 
pressions of a transfer function for closed loop linear pulse systems with constant parameters are derived. At the 
same time relations between transfer functions of closed loop and corresponding open loop systems are clarified 
and the properties of transfer functions of closed loop systems are discussed. The same questions are examined 
for systems with several pulse elements, which considerably widens the sphere of application of the theory de- 
veloped in the book. 


Great attention is paid to the important question of the stability of pulse systems. In the case of pulse systems 
with constant parameters,this question is fully dealt with in the book. Having discussed the definition of stability, 
the author relates stability with the properties of the transfer function and formulates Nyquist's criterion for pulse 
systems. Practical methods of investigating stability are examined on the basis of this criterion. The application 
of Mikhailov's criterion to pulse systems is discussed and of algebraic criteria, in particular, that of Hurwitz-Rouse. 


For the practical design of stable systems, it is important, in addition to finding the conditions of stability, 
also to take steps for providing that stability. This important question has also been dealt with in Chapter V. 
Time and frequency characteristics of closed loop systems are also examined in detail, practical methods of plot- 
ting these characteristics and evaluating by means of them the transient and stable-state processes at the output 
of the system are given. Thus, Chapter V deals in detail with all the propositions characterizing the normal pro- 
cesses in closed loop pulse systems. : 
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The analysis of random processes is just as important for pulse systems with feedback as it is for open loop 
systems. Chapter V provides such an analysis with respect to an important particular case, thatof random station- 
ary signals in systems with constant parameters. Moreover, considerable attention is paid to the synthesis of opti- 
mum linear pulse systems. This important question is considered within the framework of Wiener's theory, i. e., 
the minimum mean-square error serves as the criterion for optimizing the linear system. It should be noted that 
this is,as yet,the only criterion used in automatic control systems,and it is not always sufficient. Often optimum 
operations are not linear and there arises a necessity for more general criteria related to the theory of evaluating 
distribution parameters. Moreover, many practical problems lead to the investigation of pulse systems in the pre- 
sence of nonstationary random signals. These remarks, however, are meant less as a rebuke to the author and 
more as a wish for a further development of the theory, since in view of the incomplete development of these 
problems, the author apparently could not include them in the monograph in a sufficiently finalized form. 


Attention is also paid to closed loop nonlinear pulse systems and a particular case of inertialess nonlinearity 
is examined. It should be noted that this case covers many problems and,hence, it is of great practical interest. 


The book is concluded by Chapter VI,which provides many closed pulse system design problems whose value 


consists, as in the case of previous problems, not only in illustrating the proposed theory, but in providing practical 
design examples as well. 


Ya. Z. Tsypkin's monograph "The Theory of Pulse Systems” is a valuable contribution to the literature on 
the theory of pulse systems. It deals fully and consistently with the theory which has important practical applica- 
tions and it examines a large number of technical problems in various spheres which facilitates the assimilation 
of the basic principles of the theory in addition to the intrinsic value of the examples themselves. 


G.P. Tartakovskii 
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